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I . INTRODUCTION 

The  purpose  of  this  document  Is  two-fold.  First,  It  serves  as  a de- 
tailed introduction  to  the  General  Electromagnetic  Model  for  the  Analysis  of 
Complex  Systems  (GEMACS).  Secondly,  Its  format  and  presentation  of 
figures  are  such  that  it  can  serve  as  a teaching  guide  and  source  of 
Vugraph  slides.  With  this  In  mind  a detailed  sample  problem  Is  presented 
and  continually  referred  to  throughout  the  report.  The  geometry  of  the 
problem  Is  shown  In  Figure  1. 

It  is  assumed  that  this  geometry  represents  the  roof  of  a cormmlcatlons 
hut.  Inside  the  hut  are  located  a transmitter  and  a receiver  operating  at 
two  different  frequencies.  The  concern  Is  centered  around  the  Interference 
in  the  receiver  caused  by  radiated  energy  from  the  transmitter.  Two  anten- 
nas are  collocated  on  the  top  of  the  hut  and  are  shown  as  dots  In  Figure  1. 

It  Is  not  the  Intention  of  this  report  to  present,  derive  or  justify 
the  basic  theory  on  which  GEMACS  Is  based,  nor  is  the  structure  or  make-up 
of  the  code  Itself  discussed.  The  Intention  Is  to  show  how  the  user  Inter- 
faces with  GEMACS.  What  are  the  kinds  and  quantity  of  input  data;  what  are 
the  cormands  and  the  language  of  GEMACS?  What  are  the  data  available  from 
a GEMACS  analysis;  what  form  do  they  take;  and,  what  do  the  data  tell  the 
user? 

Section  II  gives  a very  brief  theoretical  background  of  the  method  of 
moments  (MOM)  pointing  out  the  limitations  of  many  existing  codes  Implementing 
this  formulation.  Section  III  shows  how  GEMACS  reduces  the  impact  of 
these  limitations.  This  will  be  followed  by  a discussion  of  the  language 
of  GEMACS,  the  use  of  the  language  being  Illustrated  by  examples  taken  from 


the  analysis  of  the  sample  problem.  Section  VI  describes  the  output 
provided  by  GEMACS.  Section  VII  describes  an  alternate  solution  capability 
built  Into  GEMACS  - a conventional  full  matrix  solution  formulation.  Sec- 
tions VIII  and  IX  show  how  GEMACS  is  used  to  analyze  5catter1ng  problems 
and  how  GEMACS  may  be  used  to  perform  parametric  studies  of  a system  being 
designed.  The  Appendix  presents  the  total  computer  output  from  an  analysis 
of  the  sample  problem  diagramed  In  Figure  1. 
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II.  BRIEF  THEORETICAL  BACKGROUND 


GEMACS  Is  based  on  the  Method  of  Moments  (MOM)  formulation  for  the 
solution  of  Maxwell's  equations.  This  formulation  has  been  extensively 
written  in  readily  available  open  literature.  Most  of  this  deals  with 
the  more  specialized  topics.  References  which  give  the  basic  derivation 
of  the  MOM  algorithms  are  listed  In  Section  X fl  -3], 

GEMACS  Is  a thln-wlre  code  which  means  that  the  geometry  Is  re- 
presented by  a set  of  wires  whose  radii  should  generally  be  less  than 
one-tenth  of  the  subsection  length.  As  a rule,  subsection  lengths 
should  be  one-tenth  of  a wavelength  or  less  at  the  desired  frequency. 
Somewhat  lonqer  subsections  may  be  used  on  long  wires  with  no  abrupt 
changes  while  shorter  lengths,  one-twentieth  of  a wavelength  or  less,  may 
be  required  In  modeling  the  critical  regions  of  an  antenna  (e.g.,  the  feed 
region) . 

However,  extremely  small  subsection  lengths  and  wire  radii  should  be 
avoided.  Approximate  minima  are  0.001  and  0.00001  wavelength,  respectively. 
These  are  dependent  on  the  accuracy  and  word  size  of  the  computer. 

In  addition,  the  lengths  of  adjacent  subsections  should  not  vary  by 
a ratio  greater  than  two.  This  Is  especially  Important  at  the  junctions  of 
several  wires.  This  Is  necessary  since  the  code  use$  an  Interpolation 
scheme  to  represent  the  current  variation  over  each  segment,  and  rapidly 
changlnq  lengths  for  adjacent  subsections  may  resqlt  In  unrealistic  distri- 
butions of  current. 

The  radii  of  adjacent  subsections  should  also  be  of  comparable  size. 

This  is  In  line  with  the  physical  geometry  as  well  as  allowing  for  a 


smoother  Interpolation  of  the  boundary  conditions,  which  are  applied  at 
the  surfaces  of  the  conductors. 

Finally,  an  actual  wire  in  the  system  should  be  modeled  with  its  actual 
radius.  This  is  consistent  with  the  concept  of  modeling  the  system  in  a 
manner  as  true  to-llfe  as  possible. 

If  the  geometry  consists  of  a surface  or  a solid  sheet  of  conducting 
material , then  the  surface  is  represented  by  a wire  grid  of  thin  wires. 

This  can  be  pictured  as  modeling  the  surface  by  window  screen.  Grid  models 
should  use  a wire  radius  of  one-fiftieth  of  the  segment  length.  In  addition, 
grid  mesh  circumferences  should  not  greatly  exceed  one-half  of  a wavelength. 
Larger  circumferences  could  generate  loop  resonances  which  do  not  exist  in 
the  physical  system. 

In  all  cases  segments  with  lengths  differing  by  more  than  a factor  of 
two  should  not  be  joined.  Angles  smaller  than  about  twenty  degrees  between 
joined  segments  should  also  be  avoided.  Unjoined  segments  should  be 
separated  by  at  least  the  maximum  diameter  of  any  of  the  subsections. 

These  guidelines  are  summarized  in  Figure  2. 

Usinn  the  guidelines  that  have  been  specified  on  the  problem  that  has 
been  postulated  results  in  the  wire  model  shown  In  Figure  3.  Each  subsection 
is  either  a quarter  or  an  eighth  of  a wavelenoth  in  length.  In  the  region 
around  the  antennas  the  length  around  each  loop  of  the  mesh  Is  one-half  of 
a wavelength.  The  surface  is  thus  made  up  of  70f  subsections.  The  dots 
represent  the  locations  of  the  two  antennas.  These  are  subsectioned  as  shown 
in  Figure  A.  It  is  assumed  that  one  antenna  is  a ouarter-wave  monopole  at 
the  freouency  of  interest,  and  the  second  is  a monopole  three-quarters  of 
a wavelength  long. 
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• GENERAL 

• JOIN  SIMILAR  SUBSECTIONS  WITH  LENGTHS  AND  RADII 

• AVOID  SMALL  ANGLES 

• SEPARATE  UNJOINED  SEGMENTS  BY  SUBSECTION  DIAMETER 

• WIRES 

t LENGTH  LESS  THAN  0.1  WAVELENGTH 

• RADIUS  LESS  THAN  0.1  LENGTH 

• USE  RADIUS  OF  ACTUAL  WIRES 

• AVOID  VERY  SHORT,  VERY  THIN  SUBSECTIONS 

t SURFACES 

• WIRE  RADIUS  LESS  THAN  0.02  LENGTH 

t MESH  CIRCUMFERENCE  LESS  THAN  0.5  WAVELENGTH 


SUBSECTIONING  GUIDELINES 


Figure  2 
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«ONE  EIGHTH  WAVELENGTH 
«»ONE  QUARTER  WAVELENGTH 


The  first  antenna  is  modeled  with  two  subsections,  and  the  second  Is 
modeled  with  six  segments.  One  volt  Is  applied  to  the  bottom  subsection  of 
the  first  antenna,  and  a fifty-ohm  load  Is  attached  to  the  bottom  of  the 
second  antenna.  Note  that  each  segment  has  a length  of  one-eighth  of  a 
wavelength.  It  Is  further  postulated  that  they  both  have  an  Identical 
radius  of  0 . 001 94  wavelength.  Thus  each  subsection  has  a length  around 
one-tenth  of  a wavelength,  and  the  radius  of  a subsection  Is  less  than 
one-tenth  the  length. 

One  final  point  to  note  In  Figure  3 Is  that  the  antennas  are  located 
so  that  the  bottom  segment  of  the  antenna  is  connected  to  an  intersection 
of  the  wire  grid  representing  the  surface.  If  it  were  not  done  this  way, 

. 

the  currents  on  the  antennas  would  have  no  way  of  coupling  to  the 
ground  plane  In  the  model. 

The  MOM  formulation  results  In  a matrix  equation  of  the  form  shown  In 
. Figure  5.  Here  the  Z-matrlx  represents  the  Interaction  between  subsections, 

and  Its  elements  are  a function  of  the  electrical  distance  between  pairs 
of  subsections.  There  are  N by  N such  elements,  where  N Is  the  number  of 
subsections  used  to  model  the  geometry.  In  this  case  N Is  714.  The 
I -matrix  Is  a column  vector,  each  element  belno  the  unknown  amplitude 
of  current  on  each  subsection.  This  Is  an  N by  one  matrix.  The  E-matrlx 
Is  also  a column  vector,  each  element  of  which  Is  the  amplitude  of  the 
electric  field  across  each  subsection.  This  Is  also  an  N by  one  matrix. 

' . In  this  equation  the  elements  of  the  I-matrlx  are  to  be  calculated 

knowing  the  geometry  and  having  specified  a particular  excitation.  This 
can  be  done  using  either  a matrix  Inversion  and  multiplication  process 
or  some  form  of  elimination  and  Iteration  scheme.  In  either  case,  the 
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ALL  DIMENSIONS  IN  WAVELENGTHS 

FIGURE  4.  ANTENNA  SUBSECTIONING 
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= ELECTRIC  FIELD  ON  SUBSECTION  i DUE  TO  A CURRENT 
OF  UNIT  AMPLITUDE  ON  SUBSECTION  j 

= CURRENT  ON  SUBSECTION  i 
= FIELD  ON  SUBSECTION  j 


MOM  MATRIX  EQUATION 
Figure  5 
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process  of  filling  the  Interaction  matrix  and  solving  for  the  currents 
Is  an  expensive  one.  Computer  storage  Is  also  a problem,  the  example  re- 
quiring approximately  one  million  computer  words. 

It  will  be  well  at  this  point  to  list  the  limitations  and  assumptions 
that  will  be  built  Into  a GEMACS  analysis.  These  are  shown  In  Figure  6. 

As  GEMACS  now  exists  the  geometry  Is  modeled  as  a wire  screen  approximation 
to  the  actual  surface.  The  currents  In  the  model  exist  only  on  the  axis 
of  the  wires  In  the  grid,  whereas  physically  they  are  spread  over  the  entire 
area  of  the  surface. 

It  follows  then  that  since  the  current  exists  only  on  the  axis  of  the  wire, 
there  Is  no  azimuthal  variation  of  the  current  around  the  circumference  of 
the  wire,  as  would  exist  physically  on  the  antennas.  As  an  example  of 
this  point,  consider  two  wires  located  near  each  other  such  that  the  separation 
between  them  Is  less  than  a wavelength.  Depending  on  whether  the  currents 
In  the  wires  flow  In  the  same  direction  or  In  opposite  directions,  there 
will  be  a spreading  or  a concentration  of  longitudinal  current  on  the  surface 
of  the  wires  between  the  wires,  with  the  opposite  phenomenon  occurlng  on  the 
opposite  side  of  the  respective  wires.  Thus,  there  would  exist  an  azimuthal 
variation  of  the  longitudinal  current  around  the  wires. 

Also  because  of  the  assumed  concentration  of  the  current  on  the  axis 
of  the  wires  there  Is  no  radial  component  of  current  flow  within  the  wire. 

This  has  significance  In  the  generation  of  near-field  phenomena  at  the  ends 
of  the  wire. 

Antenna  sources  are  generally  modeled  as  a delta  voltage  source  placed 
across  a subsection.  This  may  have  no  counterpart  whatsoever  In  the  physical 
situation.  Moreover,  the  size  of  the  gap  In  the  model  usually  does  not 
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ACTUAL  GEOMETRY  REPRESENTED  BY  THIN  WIRES 
CURRENT  EXISTS  ONLY  ON  WIRE  AXIS 
NO  AZIMUTHAL  VARIATION  OF  CURRENT 
NO  RADIAL  COMPONENT  OF  CURRENT 
ANTENNA  SOURCES  ARE  GAP  MODELS 
CURRENTS  DETERMINED  ONLY  AT  POINTS 


SYSTEM  LIMITED  IN  SIZE 
PERFECT  CONDUCTORS 
EXTERNAL  PROBLEMS 


MOM  ASSUMPTIONS  & LIMITATIONS 
Figure  6 
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bear  any  relation  to  the  size  of  the  gap  In  the  physical  antenna.  The 
gap  In  the  model  Is  usually  the  same  size  as  the  length  of  the  adjacent 
subsections,  since  one  of  the  modeling  rules  of  thumb  Is  to  avoid  large  ratios 
In  the  relative  lengths  of  adjacent  subsections. 

The  current  on  a subsection  Is  computed  at  the  center  of  the  segment, 
and  the  variation  over  the  wire  Is  determined  by  Interpolation  between 
adjacent  centers,  as  discussed  previously. 

Even  though  all  these  assumptions  are  built  Into  a GEMACS  analysis, 
or  a MOM  analysis  In  general,  good  correlation  exists  between  measured 
data  and  predicted  data,  and  between  other  analytical  results  and  the  data 
obtained  by  using  GEMACS  [4  - 8]. 

MOM  analyses  have  been  limited  to  small  systems,  those  that  can  be 
represented  by  300  subsections  or  less.  Electrically,  this  roughly 
corresponds  to  a size  of  30  wavelengths  of  wires  or  a surface  with  an  area 
of  one  square  wavelength.  This  has  not  been  due  to  a limitation  of  the 
theory  or  the  technique,  but  has  been  brought  about  by  the  computer  re- 
sources needed  to  perform  a MOM  analysis.  The  computer  core  storage  goes  up 
as  N ; and  the  solution  time  Increases  as  N , where  N Is  again  the  number  of 
subsections. 

All  the  commonly  used  codes  assume  that  the  material  of  which 
the  system  Is  composed  Is  perfectly  conducting.  Brute-force  techniques 
can  be  used  to  get  around  this  limitation,  but  they  put  a heavy  burden  on 
the  user.  Alternate  solutions  are  being  pursued  at  RADC  and  will  be 
Included  In  future  versions  of  GEMACS. 

finally,  the  MOM  wire  grid  model  can  be  used  only  to  solve  the 
external  problem.  Problems  that  cannot  be  treated  with  confidence 
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Include  coupling  through  apertures  In  the  skin  of  the  structure  and  coupling 
between  antennas  located  on  opposite  sides  of  the  structure.  Since  the 
structure  Is  modeled  by  a wire  grid,  electromagnetic  energy  will  "leak 
through"  the  mesh  In  the  model,  thus  resulting  In  a form  of  aperture 
coupling.  For  the  antenna  coupling  problem,  energy  will  go  directly  through 
the  body  In  addition  to  going  around  on  the  surface.  Thus,  the  coupling 
will  be  greater  than  If  the  surface  were  modeled  as  a solid.  A solid 
surface  model  will  be  shortly  Implemented  Into  the  GEMACS  code  to  eliminate 
these  present  limitations. 

With  this  brief  background  It  Is  now  time  to  examine,  again  briefly,  how 
GEflACS  attempts  to  relieve  the  problem  of  large  computer  resources  being 
reguired  to  perform  an  analysis  on  a large  problem. 


I 
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III.  GEMACS 


As  previously  mentioned,  the  MOM  formulation  results  In  a matrix 
equatlpn  which  Is  solved  to  obtain  the  current  distribution  on  the  struc- 
ture. The  solution  is  usually  obtained  by  matrix  Inversion  or  some 
elimination  technique.  The  drawback  Is  the  time  required  to  accomplish 
the  solution. 

GEMACS  gets  around  this  problem  to  a certain  extent  by  employing 
what  Is  called  the  Banded  Matrix  Iteration  (BMI)  technique  [9].  The 
equations  are  shown  in  Figure  7. 

To  start  with,  recall  that  the  Z-matrlx  elements  correspond  to  Inter- 
actions between  wire  subsections.  The  Interactions  decrease  with  Increas- 
ing distance  between  the  subsections.  As  will  be  shown  later,  the 
segments  can  be  numbered  such  that  the  difference  between  subsection  num- 
bers for  all  close-neighboring  segment  pairs  Is  small  compared  to  N.  The 
largest  matrix  elements  can  then  be  kept  close  to  the  principal  diagonal 
of  the  Interaction  matrix,  shown  dotted  In  Figure  7. 

The  first  equation  on  Figure  7 Is  the  basic  matrix  equation  to  be 
solved,  relating  the  field  on  the  subsections  to  the  currents  on  the  sub- 
sections by  the  Interaction  matrix.  The  Interaction  matrix  Is  then 
separated  Into  three  matrices  as  shown  In  the  second  equation  In  the 
figure.  Here  B Is  called  a banded  matrix  with  lower  and  upper  bandwldths 
M.  That  Is,  M Is  the  number  of  minor  diagonals  on  either  side  of  the 
main  diagonal.  L and  U contain  the  remaining  elements  of  the  Interaction 
matrix  and  are  called  the  lower  and  upper  matrices,  respectively.  They 
are  both  triangular  matrices. 
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BI  = E-(L+U)I 


BMI  EQUATION  DERIVATION 


FIGURE  7. 


These  two  equations  can  be  combined,  after  some  minor  manipulation,  to 
result  In  the  third  equation  shown  In  Figure  7.  An  Iterative  scheme 
Is  shown  In  the  first  equation  In  Figure  8,  where  Ij  denotes  an  approximation 
to  the  solution  at  the  jth  Iteration;  and  Ij+1  Is  an  approximation  to 
the  solution  at  the  following  Iteration.  Some  starting  value  IQ  (usually 
zero)  Is  chosen,  and  Ij  Is  computed  from  this  equation.  Then  Ij  Is  entered 
on  the  right-hand  side,  and  ^ is  computed.  If  the  sequence  converges,  then 
an  approximate  solution  to  the  original  matrix  equation  Is  obtained. 

This  equation  must  be  solved  at  each  Iteration.  The  cost  Is  minimized 
by  decomposing  B Into  a product  of  lower  and  upper  triangular  matrices. 

Then  the  first  equation  In  Figure  8 Is  solved  by  forward  elimination,  as 
shown  In  the  second  equation  In  the  figure.  Finally,  this  Is  followed 
by  a back  substitution,  as  shown  In  the  third  equation  In  Figure  8. 

Note  that  the  decomposition  of  the  banded  matrix  Is  similar  to  the 
full  matrix  decomposition,  except  that  the  cost  Is  much  less  since  this  Is 
a matrix  In  which  many  of  the  elements  are  zero.  The  cost  for  the 
Iterative  solution  process  Is  shown  In  the  last  equation  In  Figure  8,  where 
K Is  the  number  of  Iterations,  N Is  the  rank  of  the  Interaction  matrix 
and  M Is  the  bandwidth  chosen.  Note  that  the  best  general  method  for  the 
solution  of  the  matrix  equation  requires  n3/3  operations.  In  GEMACS  the 
user  must  make  a trade-off  between  the  number  of  Iterations  and  the  bandwidth. 
The  larger  the  bandwidth,  the  smaller  the  number  of  Iterations  and  vice 
versa.  No  rules  of  thumb  exist  for  the  general  case  to  optimize  the  cost 
of  solving  the  matrix  equation.  However,  It  Is  apparent  that  one  should 
choose  a large  bandwidth.  This  Is  so  slnce’a  larger  bandwidth  provides 
lower  Initial  errors  and  faster  convergence  rates.  If  the  user  picks  a 
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bandwidth  too  small,  slow  convergence  and  low  efficiency  with  respect  to 
the  Gauss-Oordan  elimination  technique  will  result.  Restarting  with  a 
larger  bandwidth  might  provide  a higher  overall  efficiency,  even  with  the 

. 1 

cost  of  the  aborted  run  Included.  This  topic  Is  discussed  In  the  BOM 
Interim  reports  [4-8],  especially  RADC-TR-75-272  [6]. 

So  much  for  what  the  problem  Is,  and  how  and  why  GEMACS  solves  the 
problem.  What  Is  GEMACS,  at  least  In  Its  overall  concepts?  Basically, 

GEMACS  Is  a mainframe  data  handler,  whose  main  purpose  Is  to  store  data 
and  present  data  In  a place  where  they  are  needed  and  In  a format  required  by 
some  operational  routine.  In  this  respect  It  can  be  used  In  many  problems 
that  result  In  large  amounts  of  data,  such  data  being  used  many  times  over. 

A functional  breakdown  of  the  GEMACS  structure  Is  shown  In  Figure  9. 

The  GEMACS  executive  routines  control  the  Interface  of  the  code  with  the 
host  computer  and  perform  three  basic  functions:  Input/output  to  peripheral 
files,  taking  checkpoints  and  restarting  from  these  checkpoints;  and  the 
compilation  of  statistical  Information,  which  can  be  used  to  pinpoint 
areas  for  further  code  refinement. 

The  Input  language,  task  execution  and  run  termination  processors 
simply  read  the  user's  data  deck,  call  appropriate  subprocessors  based 
on  the  user's  commands  and  terminate  the  analysis,  respectively.  Within 
these  two  unper  levels  In  Figure  9 are  all  of  the  file  handling  capabilities 
built  Into  GEMACS.  With  a proper  Interface  for  the  new  subprocessors 
under  the  task  execution  processor,  these  new  subprocessors  can  change 
the  field  analysis  technique  or  even  apply  the  mainframe  to  a different 
type  of  problem  completely,  such  as  the  dynamic  load  analysis  of  some 
structure. 
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The  geometry  processor  generates  the  geometry  to  be  analyzed 


by  Interpreting  the  user  geometry  Input  commands.  The  Interaction  matrix 
processor  generates  the  elements  of  the  Interaction  matrix  for  the  fre- 
quency and  geometry  specified  by  the  user.  The  excitation  processor  gen- 
erates the  elements  of  the  vector  on  the  right-hand  side  of  the  MOM  ma- 
trix equation.  The  load  processor  modifies  the  Interaction  matrix  to  take 
Into  account  the  presence  of  loads  on  the  wires  or  the  presence  of 
Imperfectly  conducting  materials.  The  matrix  solution  processor  solves 
the  MOM  matrix  equation  for  the  currents  on  the  structure.  The  output 
processor  calculates  such  quantities  as  the  near-  and  far-fleld  patterns 
and  terminal  Impedance  for  antennas  or  the  backscatterlng  from  the  structure. 
The  direct  manipulation  processor  sets  such  variables  as  the  maximum 
CPU  time  allowed  for  the  analysis,  the  number  of  files  In  the  system 
available  to  the  code,  the  frequency  of  the  analysis  and  the  electrical 
character! si  tics  of  the  ground  (If  present).  It  also  performs  arithmetic 
operations,  such  as  modifying  the  frequency  by  some  factor,  a feature  which 
Is  useful  when  a loop  Is  being  built  Into  the  command  stream.  This  Is 
Illustrated  by  examples  In  Section  IX. 

What  Is  the  Impact  of  all  this  to  the  user?  This  modular  con- 
struction has  several  Important  advantages  that  he  should  be  aware  of 
at  the  start.  First,  and  most  obvious,  is  that  one  can  plug  In  any  tech- 
nique which  has  a proper  Interface  or  driver  to  transfer  data  between  the 
implementing  subroutine  and  the  mainframe.  This  Is  much  like  the  use  of 
magnetic  cards  for  programs  on  personal  calculators.  It  Is  therefore  possible 
to  have  a complete  set  of  techniques  stored  In  separate  files.  But  the 
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big  difference  Is  that  communication  with  all  these  techniques  Is  In  one 
common  language.  There  Is  no  need  to  remember  several  different  sets  of  In- 
put formats,  or  limit  oneself  to  one  specialized  code. 

A wealth  of  debug  features  are  built  Into  the  code  to  aid  In  the  exe- 
cution of  an  analysis  or  to  aid  In  the  Integration  of  new  analysis  tools. 

An  explanation  of  the  debug  messages  accompanies  the  code  documentation 
supplied  with  GEMACS. 

GEMACS  has  a checkpoint/restart  capability  within  the  code.  This  means 
that  long  analyses  can  be  stopped  and  started  In  case  the  computer  should 
go  down  or  In  case  there  Is  a need  to  change  the  conmand  stream,  as,  for 
example.  In  the  case  In  which  the  bandwidth  of  the  banded  matrix  needs  to  be 
changed  to  decrease  solution  time  or  the  number  of  Iterations.  This 
feature  of  restarting  the  same  analysis  several  times  can  be  used  to  schedule 
the  analysis  such  that  It  Is  accomplished  at  hours  In  which  the  computer  load 
Is  a minimum. 

Finally,  and  highly  significant.  Is  the  fact  that  GEMACS  Is  tied 
Into  the  Air  Force  Intrasystem  Analysis  Program  (IAP).  This  means  that  It 
will  have  the  full  support  of  the  Air  Force,  to  provide  aid  In  the  loading, 
use  and  maintenance  of  the  code,  additions  to  the  capability  of  the  code, 
updates  to  eliminate  any  bugs  that  may  be  In  the  code;  and  a corrmon  language 
among  all  users  of  the  code. 

This  then  has  been  a very  brief  Introduction  to  the  problem  and  a 
technique  to  solve  this  problem.  A quick  run-down  on  the  structure  of 
GEMACS  and  the  advantages  of  Its  use  to  the  user  have  also  been  given.  It 
Is  now  time  to  get  acquainted  with  the  language  of  GEMACS  and  to 
see  how  the  example  problem  would  be  solved  using  GEMACS. 
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IV  SEHACS  GEOMETRY  LANGUAGE 

The  function  of  the  geometry  processor  shown  In  the  block  diagram 
of  Figure  9 Is  to  translate  the  user's  Input  regarding  the  geometry  Into  a 
data  set  which  may  then  be  operated  on  by  the  Interaction  matrix  generator 
to  generate  the  Interaction  matrix  for  the  structure  under  analysis. 

The  geometry  processor  Is  entered  when  the  GMDATA  command  Is  encountered 
In  the  command  language  stream  (see  Section  V.C). 

The  basic  elements  for  GEMACS  are  points  and  line  segments.  These 
In  turn  may  belong  to  larger  groups  with  a given  name,  such  as  WING  or 
PLATE.  Any  reference  to  this  given  name  will  also  reference  all  the  points 
and  segments  within  that  group.  In  addition,  line  segments  may  also  be 
Identified  as  a group  by  having  the  same  tag  number.  Thus,  all  the 
segments  belonging  to  the  wing  may  have  a tag  number  "six"  and  may  be 
referenced  as  a group  by  specifying  a tag  of  six  Instead  of  using  the  label 
WING.  Another  use,  as  shown  later,  would  be  to  give  all  segments  loaded 
by  the  same  Impedance  an  Identical  tag  number.  This  will  obviate  the  need 
for  multiple  load  comnands.  The  same  Is  true  for  excited  segments. 

The  subsection  may  thus  be  Identified  by  either  a segment  number  or 
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a tag  number.  The  difference  Is  that  the  first  Is  unique  In  the  model 


while  the  second  may  be  shared  by  any  number  (or  all)  of  the  subsections 
within  the  geometry  model. 

The  following  basic  conventions  must  be  observed  when  the  geometry 
deck  Is  being  built.  The  first  nonblank  Item  (not  necessarily  the  first 
two  columns)  of  the  card  must  be  a two-letter  alpha  code.  Then  the  data 
following  this  type  code  must  be  separated  Into  Items  by  either  a comma 
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or  a blank.  Blanks  within  Items  are  not  allowed,  as  the  data  would 
then  be  Interpreted  as  two  separate  pieces  of  Information  when  In  fact 
only  one  was  meant. 


If  a dollar  sign  ($)  Is  encountered  on  a card,  the  data  following  are 
treated  as  a comment  and  are  ignored  for  further  processing.  Also,  keep 
In  mind  that  should  a comment  be  used  on  a card,  a blank  must  separate 
the  dollar  sign  from  the  last  piece  of  geometric  Input.  A dollar  sign  as 
the  first  nonblank  character  on  a card  will  cause  the  entire  card  to  be 
treated  as  a comment  card,  and  the  data  on  that  card  will  be  Ignored 
for  further  processing.  Comments  are  useful  when  a second  party  Is 
using  the  deck.  This  allows  easy  Interpretation  of  what  Is  being  done 
with  each  command. 

Data  for  the  same  geometry  command  may  appear  on  more  than  one  card. 
This  Is  accomplished  through  the  use  of  continuation  cards.  A continua- 
tion card  is  Indicated  by  the  presence  of  an  asterisk  (*)  In  column  one 
of  each  succeeding  card.  There  Is  no  limit  to  the  number  of  continuation 
cards  that  may  be  used,  providing  that  no  more  than  256  pieces  of  data  are 
associated  with  any  single  command.  About  the  only  comnand  that  this  will 
occur  In,  as  will  be  seen.  Is  the  renumber  (RN)  command  (IV. E). 

The  general  form  of  the  conmands  Is  shown  In  Figure  10.  In  this  figure 
TYPE  is  one  of  the  mnemonic  codes  In  Figure  11;  and  PI,  P2,  etc.  are  the 
ordered  parameter  Items  required  for  the  processing  of  the  command.  As 
shown  In  the  figure  they  are  separated  by  commas  or  by  blanks.  Note  that 
the  order  of  the  appearance  of  the  data  Items  Is  fixed,  and  any  variation 
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1.  GENERAL  COMMANDS 

CS  RA 

END  SC 

2.  SUBSECTION  GENERATION 

CP  PT 

MP  WR 

3.  ELEMENT  GENERATION 

AT  DE 

CE  DF 

4.  MACROGENERATORS 

RF 

RT 

XL 

5.  RENUMBER 

RN 


Figure  11.  GEOMETRY  GENERATION  GROUPING 
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will  cause  an  abort  of  the  analysis  or  an  Incorrect  representation  of  the 
geometry. 

Finally,  some  of  the  conmands  have  Items  which  may  be  defaulted. 

These  are  pointed  out  In  the  discussion,  and  their  existence  Is  Indicated 
by  the  use  of  brackets  [ ] In  the  figures.  In  general,  for  a command 
of  the  form  shown  In  Figure  10  In  which  P3,  P4  and  P5  may  be  defaulted. 

It  must  be  kept  In  mind  that  the  defaults  can  be  achieved  only  from  right 
to  left.  For  example,  to  default  P5,  both  P3  and  P4  must  be  specified. 

To  default  P4,  only  P3  needs  to  be  specified,  while  P5  must  assume  Its 
default  value.  When  one  or  more  parameters  are  to  be  defaulted,  then 
all  Items  to  their  right  must  also  be  defaulted.  As  a final  example 
assume  that  It  Is  desired  to  default  P3.  Then  both  P4  and  P5  must  also 
be  defaulted.  If  P4  must  be  specified,  then  a value  for  P3  must  be 
also  Input.  If  a value  for  P5  must  be  given,  then  both  P3  and  P4 
must  be  Input. 

The  geometry  commands  may  be  arbitrarily  grouped  Into  five  categories 
as  shown  In  Figure  11.  The  general  commands  essentially  set  what  may 
be  called  the  global  parameters  of  the  problem.  They  are  defined  In  the 
beginning  of  the  geometry  deck  and  referenced  thereafter.  The  only 
exception  Is  the  END  card,  which  Is  the  last  card  of  the  deck  and  Indicates 
to  the  geometry  processor  that  It  now  has  processed  all  of  the  Input  cards 
relating  to  the  geometry. 

The  basic  units  of  geometry  are  points  and  subsections  and  these  are 
defined  by  the  comnands  of  group  two.  Large  units  can  be  built  from  these 
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points  and  lines  by  defining  and  combining  elements  using  the  commands  of 
group  three.  These  units  may  be  generated  quite  easily  by  the  use  of  the 
reflect,  rotate  and  translate  commands  of  group  four. 

The  renumber  command  can  be  used  to  resequence  the  numbering  of  the 
subsections.  Thus,  the  geometry  may  be  entered  In  any  convenient  manner 
and  then  very  simply  be  resequenced  In  order  to  rearrange  the  Interaction 
matrix  In  preparation  for  use  of  the  BMI  solution  technique.  Examples  of 
this  will  be  shown  during  the  discussion  of  the  renumber  comnand. 

A.  General  Commands 

The  four  commands  within  this  group  perform  the  following  functions: 

1.  Define  all  coordinate  systems  (CS) 

2.  Set  up  a radius  table  for  all  subsections  (RA) 

3.  Scale  the  geometry  data  to  meters  (SC) 

4.  Define  the  end  of  the  data  deck  (END). 

The  first  two  commands  will  always  be  located  at  the  beginning  of  the  geo- 
metry data  deck.  The  SCALE  command  appears  before  a deck  of  cards  with 
consistent  units  which  are  not  meters.  If  a deck  of  cards  then  follows 
with  a different  set  of  units,  this  latter  deck  must  be  preceded  by  an 
appropriate  SCALE  command.  The  END  command  Is  the  last  card  of  the  geometry 
deck  at  all  times. 

Figure  12  again  shows  the  wire  grid  model  of  the  top  of  the  hut. 

Also  shown  outlined  is  an  area  of  finer  mesh  which  will  be  defined  In  a 
separate  coordinate  system.  Figure  13  shows  the  general  format  for  the 
Coordinate  System  (CS)  specification  comnand,  as  well  as  specific  examples 
suitable  for  the  coupling  problem.  The  mnemonic  CS  identifies  this 
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as  a coordinate  system  specification  command.  The  parameter  NCS  Is  a 
unique  Integer  Identifier  of  the  coordinate  system  being  defined.  In  this 
example  the  value  "two"  has  been  chosen  for  the  fine  mesh  area.  The 
parameters  XC,  YC  and  ZC  are  the  (x,  y,  z)  coordinates  of  the  global 
coordinate  system  at  which  the  origin  of  this  local  coordinate  system  Is  to 
be  located.  As  seen  from  figure  12  XC-0.0,  YC-0.3  and  ZOO. 2,  where  the 
units  are  In  meters.  The  final  set  of  parameters  RX,  RY  and  RZ  are  the 
rotation  angles  In  degrees  about  the  global  coordinate  system  X,  Y,  and  Z 
axes.  An  example  of  this  can  be  seen  on  page  61  of  the  user's  manual  [10] 
for  GEMACS.  Note  that  for  this  command  none  of  these  parameters  may  be 
defaulted. 

It  should  be  noted  that  when  more  than  one  coordinate  system  Is  present, 
even  the  global  system  should  be  given  an  Integer  Identifier.  This  Is  due 
to  the  fact  that,  as  will  be  seen,  the  coordinate  system  parameter  on 
all  other  commands  can  be  defaulted,  but  the  default  will  not.be  to  the 
global  system.  It  will  be  to  the  last  previous  coordinate  system  refer- 
enced. Therefore,  the  first  example  In  Figure  13  Identifies  the  global 
coordinate  system  as  system  number  one.  All  the  parameters  with  the  excep- 
tion of  NCS  are  zero. 

The  only  time  the  global  coordinate  system  Is  not  labeled  Is  when 
all  of  the  geometry  will  be  defined  In  a single  coordinate  system. 

The  Radius  (RA)  specification  command  sets  up  an  Internal  table  which 
contains  the  radius  In  meters  of  each  of  the  subsections  used  In  the  analy- 
sis. This  means  that  the  user  needs  to  specify  each  different  radius 
only  once  and  then  references  these  radii  by  their  position  In  the  radius 
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command  card.  Up  to  ten  different  radii  may  be  specified. 


, 


i 


Figure  14  shows  the  general  form  of  the  command.  The  letters  RA 
Identify  this  as  the  redlus  specification  command.  The  symbols  R1  through 
Rn  represent  the  different  radii  that  will  be  used  In  defining  the 
geometry.  For  the  example  only  one  radius  Is  being  used  and  that  Is  equal 
to  0.000776  meters.  Examples  of  references  to  the  radius  table  will  be 
shown  during  the  discussion  of  other  commands.  Note  that  there  are  no 
default  parameters  for  this  command.  , 

f 

The  Scale  (SC)  specification  command  allows  the  user  to  Input  the 
geometry  data  In  any  convenient  set  of  units  and  then  modify  these  data 
by  a scale  factor  to  convert  the  units  of  length  to  meters.  If  all  data 
are  Input  in  meters,  then  this  command  Is  not  needed  at  all.  Conversely, 

If  this  command  Is  not  present,  GEMACS  will  assume  that  the  Input  data 
are  In  meters. 

The  general  format  for  this  command  and  a specific  example  are  shown 
In  Figure  15.  Note  that  the  example  Is  fictitious  and  not  related  to 
the  problem  being  used  as  an  example.  The  mnemonic  SC  Identifies  the 
command  type.  The  parameter  may  be  defaulted  (with  a default  of  one). 

This  would  indicate  that  the  following  geometry  data  are  In  meters. 

When  this  command  Is  used,  a double  option  Is  available.  First,  a 
conversion  factor,  that  Is,  a scale  value,  may  be  specified  which  would 
convert  the  Input  units  to  meters.  Second,  for  the  convenience  of  the  user 
conversion  factors  for  units  of  feet  (FT),  Inches  (IN)  and  centimeters  (CM) 
are  provided  within  GEMACS.  Then  It  Is  necessary  to  Input  only  the  proper 
two  letter  code  applicable  to  the  units  being  used. 
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RA  R1  R2  R3  •••  RN 
RA  0.000776 


Figure  14.  RADIUS  SPECIFICATION  COMMAND 


The  example  shown  In  Figure  15  demonstrates  the  use  of  this  coimand.  v 
Assume  that  the  fine  mesh  of  Figure  12  Is  dimensioned  In  centimeters,  while 
the  coarse  mesh  has  units  of  feet.  (Keep  In  mind  that  this  Is  not  the 
case  for  this  problem.  All  dimensions  are  In  meters  In  the  problem.) 

Then,  the  card  "SC  CM"  would  be  placed  before  the  geometry  Input  relating 
to  the  fine  mesh  part  of  the  structure.  Later,  when  the  coarse  mesh  model 
Is  being  Input,  the  card  "SC  FT"  would  appear  just  before  the  data  generating 
this  coarse-mesh  wire  grid. 

As  an  example  of  using  the  default  option  suppose  that  the  fine  mesh 
data  were  In  units  of  centimeters  and  the  coarse  mesh  data  were  In  units 
of  meters.  The  first  scale  card  would  read  "SC  CM"  while  the  second 
would  just  read  "SC".  The  second  card  would  be  needed,  since.  If  It  were 
not  present,  then  GEMACS  would  consider  the  coarse  mesh  data  to  be  In 
centimeters  also. 

The  final  general  card  Is  the  last  card  to  appear  In  the  geometry 
data  Input  deck.  This  Is  the  END  command.  It  signals  the  geometry 
processor  to  stop  reading  Input  since  no  more  geometry  commands  exist. 

There  Is  one  further  point  to  note  regarding  this  card.  The  user 
does  not  need  the  dollar  sign  to  delineate  the  presence  of  a comment  In 
the  command.  Figure  16  shows  the  general  format  of  the  coimand  as  well 
as  a specific  example  showing  how  the  end  of  the  geometry  Input  data  Is 
typically  Identified. 

B.  Subsection  Generation 

One  of  the  key  points  In  the  MOM  technique  is  the  Interaction  between 
subsections  or  segments.  In  GEMACS  the  physical  structure  Is  divided 
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END 


END  OF  GEOMETRY  HUTTOP 


Figure  16.  END  OF  GEOMETRY  COMMAND 


4-14 


up  Into  a number  of  these  subsections.  There  are  two  basic  ways  of 
generating  these  subsections.  One  way  Is  by  a combination  of  point 
generation  (PT),  connect  point  (CP)  and  multiple  point  connection  (MP) 
commands.  The  second  way  Is  to  generate  subsections  Individually  or  as 
wires  with  a specified  number  of  equal -length  subsections  using  the  wire 
Input  (WR)  command.  Each  of  these  will  be  discussed  In  turn,  using  examples 
that  could  appear  In  the  Input  deck  for  the  example  problem. 

The  Point  (PT)  specification  command  Is  probably  the  most  funda- 
mental geometry  card,  although  It  Is  seldom  used  since  the  other  commands 
are  more  efficient  for  generating  the  geometry.  The  general  form  Is 
shown  In  Figure  17.  The  code  PT  Identifies  the  type  of  the  conmand.  The 
parameter  NPT  Is  an  Integer  Identifier  of  the  point.  Note  that  this 
must  be  an  unique  Identifier  In  the  global  sense.  That  Is,  when  a point  Is 
given  a ntmber,  no  other  point,  even  If  It  Is  In  a different  coordinate 
system,  can  have  that  same  Identifier.  The  reason  Is  that  GEMACS  eventually 
translates  all  points  Into  the  global  coordinate  system.  The  parameters 
X,Y,Z  locate  the  point  within  the  coordinate  system  specified  by  the 
parameters  NCS. 

Care  must  be  exercised  when  using  the  default  option  If  more  than  one 
coordinate  system  Is  specified.  As  pointed  out  previously,  GEMACS  will 
default  to  the  last  coordinate  system  specified  on  any  command  and  not  to 
the  global  coordinate  system.  This  further  Implies  that  the  first  command 
In  a sequence  that  references  a particular  coordinate  system  must  not 
have  a default  for  the  NCS  field.  All  subsequent  commands  referring  to 


t 
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PT  NPT  X Y Z [NCS] 

PT  1 0.0  0,0  0.0  1 

PT  2 0.0  0.0  0.0  2 


(fictitious) 

(fictitious) 


Figure  17.  POINT  SPECIFICATION  COMMAND 


this  same  coordinate  system  may  however  exercise  the  default  option  until  a 
point  In  some  other  system  Is  generated.  This  new  coordinate  system  must  be 
Identified  In  this  latter  command. 

Two  examples  are  shown  In  Figure  17,  each  of  which  specifies  a point 
at  the  origin  of  each  coordinate  system  defined  In  Figure  13.  Again, 
these  examples  are  labeled  as  fictitious  since  In  the  sample  Input  deck 
the  wire  specification  command  Is  used  rather  than  the  point  connand. 

All  the  Intersections  of  the  wire  grid  model  shown  In  Figure  12  could 
be  specified  by  points  using  the  point  specification  command.  The  sub- 
sections joining  these  points  would  then  be  generated  using  either  a 
Connect  Points  (CP)  command  or  a multiple  point  connection  (MP)  command. 

The  connect  points  command  is  shown  In  Figure  18  along  with  a fictitious 
example.  It  Is  used  to  generate  a specified  number  of  subsections  with 
a specified  tag  Identifier  and  radius  between  two  previously  defined  points. 
Note  that  these  points  need  not  be  located  In  the  same  coordinate  system. 

The  mnemonic  CP  Identifies  this  as  the  connect  points  command.  The 
parameters  N1  and  N2  are  the  Integer  Identifiers  of  the  first  and  second 
points,  respectively.  Note  that  the  wire  connecting  these  points  Is  a 
directed  line  segment  with  a reference  direction  going  from  point  1 to 
point  2.  The  positive  and  negative  signs  on  the  computed  currents  (pages  A-42 
to  A-47  ) Indicate  whether  the  current  Is  flowing  In  the  direction  of  the 
line  segment  or  against  It,  respectively. 

The  last  three  parameters  may  be  defaulted.  Here  again  they  will 
each  default  to  the  last  corresponding  value  specified  on  any  previous 
command.  The  value  of  the  parameter  NSEG  Is  the  number  of  subsections 
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CP 

N1 

N2 

[nseg] 

RA 

0.001 

0.05 

0.015 

PT 

1 

0.0 

0.0 

PT 

2 

0.0 

0.0 

CP 

1 2 

2 1 

3 

© 


o.o 

(FICTITIOUS) 


(0,  0,  1) 


(0,  0,  0) 


Figure  18.  CONNECT  POINTS  COMMAND 
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or  segments  Into  which  the  wire  connecting  N1  and  N2  Is  to  be  divided. 

NTA6  would  be  the  tag  Identifier  of  each  of  the  subsections  on  the  wire. 

The  value  of  NRAD  Is  the  location  number  of  the  appropriate  radius  for 
this  wire  In  the  list  of  radii  given  In  the  radius  (RA)  comnand. 

The  example  shown  In  Figure  18  Illustrates  what  could  be  the  first 
four  cards  of  a geometry  Input  deck.  This  Is  again  a fictitious  example 
with  respect  to  the  sample  problem.  The  radius  table  contains  four  en- 
tries. The  second  and  third  cards  would  define  the  two  points  that  are  to 
be  connected.  The  fourth  card  Is  the  connect  points  command.  It 
specifies  that  points  one  and  two  are  to  be  connected  by  a wire  that  Is 
to  be  subdivided  Into  two  subsections,  each  subsection  having  a tag 
Identifier  of  one.  The  radius  of  the  wire  Is  found  as  the  third  entry  In 
the  list  given  In  the  radius  (RA)  comnand,  or  fifteen-thousandths  of  a meter. 
The  result  of  this  operation  Is  shown  In  the  diagram  In  Figure  18. 

The  figure  also  Illustrates  how  GEMACS  will  number  the  subsections 
that  are  generated  by  the  various  geometry  commands.  As  pointed  out 
previously,  the  line  Is  directed  from  point  N1  to  point  N2  by  the  order 
of  the  points  given  In  the  connect  points  command.  Therefore,  the  two 
subsections  that  are  generated  will  be  numbered  In  sequence  along  the 
wire  from  one  to  two.  That  Is,  segment  one  will  extend  from  point  one 
to  the  midpoint  of  the  wire.  Then  segment  two  will  extend  from  the  mid- 
point of  the  wire  to  point  two.  These  segment  numbers  are  shown  In  the 
circles  In  the  diagram  In  Figure  18. 
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A generalization  of  the  connect  points  operation  Is  to  connect  a 


series  of  points.  This  Is  accomplished  by  the  Multiple  Point  (HP) 
connection  command  shown  In  Figure  19.  In  this  case  the  points  to  be 
connected  are  joined  by  a single  unbroken  line,  but  the  line  may  have 
any  number  of  bends.  The  comnand  requires,  however,  that  each  wire  be- 
tween each  pair  of  points  be  subdivided  Into  an  equal  number  of  subsections. 
That  Is,  If  three  points  are  to  be  connected  by  two  wires,  the  wire 
between  pair  one  to  have  three  segments  and  the  wire  between  pair  two  to 
have  four  segments,  then  two  separate  commands  would  be  needed  since  one  of 
the  parameters  In  this  command  Is  a single  entry  for  the  number  of  sub- 
sections Into  which  each  wire  Is  divided. 

Looking  at  the  command  In  detail.  It  Is  seen  that  MP  Is  the  two-letter 
mnemonic  for  this  operation.  The  value  for  NPTS  Is  the  number  of  points 
that  are  to  be  connected.  In  the  command  there  must  be  present  this  number 
of  point  Integer  Identifiers.  For  example.  If  points  one  through  four 
are  to  be  connected,  then  NPTS  would  be  four;  and  this  would  be  followed 
by  the  Integer  Identifiers  of  each  of  those  four  points,  given  In  the  order 
In  which  the  points  are  to  be  connected  together.  As  explained  previously 
and  will  be  shown  In  an  example,  the  order  In  which  the  pairs  of  points 
are  connected  will  determine  the  numbering  of  the  subsections. 

The  values  for  NP1  through  NPNPTS  would  be  the  Integer  identifiers 
of  the  points  that  are  to  be  connected  together  during  the  operation 
Initiated  by  this  command. 
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The  final  three  parameters  In  this  command  are  the  number  of  subsections 
between  each  pair  of  points,  the  tag  identifier  associated  with  each  sub- 
section and  the  radius  of  the  wires  connecting  the  points.  Note,  again, 
that  the  wires  between  each  pair  of  points  will  be  divided  Into  Identical 
numbers  of  subsections;  and  each  subsection  will  have  the  same  tag  Iden- 
tifier and  the  same  radius. 

As  In  the  connect  points  command  each  of  the  last  three  parameters 
may  be  defaulted.  The  default  value  would  be  the  last  previous  entry 
of  the  respective  parameter  on  any  previous  command. 

Figure  19  shows  an  example  of  the  use  of  this  command.  Assume  that 
six  points  have  already  been  defined  using  the  point  specification  command 
and  that  the  radius  table  has  been  established.  The  example  conmand  states 
that  five  points  are  to  be  connected  together  In  the  order  one  to  two,  two 
to  three,  three  to  four  and  four  to  six.  Each  connecting  wire  will  have 
one  subsection,  a tag  of  one  and  a radius  equal  to  the  value  of  the  first 
entry  In  the  radius  table.  The  resulting  subsections  and  their  subsection 
numbers  (In  circles)  are  shown  in  the  diagram  In  Figure  19. 

Note  that  points  four  and  five  cannot  be  connected  together  using  this 
conmand  as  It  stands.  A separate  conmand  would  be  needed  to  generate  a 
wire  between  these  points.  An  alternate  approach  would  be  to  connect  points 
one  to  five  together  using  a multiple  points  connect  command,  followed 
by  a separate  command  to  connect  points  four  and  six. 


This  then  completes  the  discussion  of  one  way  of  generating  wires 
and  subsections  within  GEMACS.  It  Is  a two-step  process,  first  defining 


the  points  and  then  connecting  the  points.  It  Is  not  the  recommended  proce- 
dure, and  It  should  be  used  only  sparingly.  What  Is  perhaps  the  most  funda- 
mental command  contained  In  the  GEMACS  geometry  language  will  now  be 
described. 

The  format  for  the  Wire  (WR)  Input  command  is  shown  In  Figure  20. 

The  use  of  this  command  Is  the  preferred  method  to  Input  wire  segment  data 
Into  GEMACS.  Study  of  the  Input  data  deck  for  the  solution  of  the 
sample  problem  (pages  A-12  to  A- 1 3 ) will  show  that  It  Is  the  only  method 
used  to  Initially  define  wire  subsections. 

The  mnemonic  WR  signifies  that  this  Is  the  wire  Input  command.  The 
values  of  XI,  Y1  and  Z1  are  the  coordinates  of  the  first  end  of  the 
wire.  This  Is  what  has  been  called  point  1 in  the  connect  points  (CP) 
command.  The  values  of  X2,  Y2  and  Z2  are  the  coordinates  of  the  second 
end  of  the  wire,  or  point  2. 

The  next  three  parameters  again  specify  the  number  of  equal  length 
subsections  Into  which  the  wire  will  be  divided,  the  tag  Identifier 
of  each  of  those  subsections  and  the  location  In  the  radius  table  of  the 
radius  of  each  of  the  subsections.  As  before,  each  of  these  parameters 
may  be  defaulted,  the  default  being  to  the  last  previous  entry  of  the 
respective  parameter  on  any  earlier  command. 

The  final  parameter  In  this  command  is  the  specification  of  the 
coordinate  system  In  which  the  points  and  the  wire  exist.  In  this  case 
both  of  the  points  must  lie  In  the  same  coordinate  system.  It  should  be 
noted  that  this  Is  not  necessary  when  the  connect  points  (CP)  command 
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Is  used.  The  difference  Is  that  the  latter  command  references  point 
numbers  (which  are  globally  unique).  In  the  wire  Input  command  only  the 
(x,y,z)  coordinates  of  the  endpoints  of  the  wire  are  specified. 


The  coordinate  system  specification  may  be  defaulted.  Here  again, 
as  previously  pointed  out,  care  must  be  exercised  when  more  than  one 
coordinate  system  Is  used  In  an  analysis.  As  usual,  the  default  will  be 
to  the  last  coordinate  system  specified  on  any  previous  command. 

Two  examples  are  also  shown  In  Figure  20.  Each  example  consists 
of  two  cards.  Each  pair  of  cards  will  eventually  be  used  to  develop 
the  fine-mesh  grid  and  the  coarse-mesh  grid,  respectively.  The  results 
are  shown  as  the  dark  lines  In  Figure  21. 

A note  of  explanation  regarding  Figure  21  Is  appropriate  at  this  point.’ 
Recall  from  Figure  12  that  the  coarse-mesh  model  Is  defined  In  coordinate 
system  number  one.  The  coordinate  designations  at  the  left  and  at  the 
bottom  of  figure  21  refer  to  this  coordinate  system.  The  fine-mesh  model 
Is  defined  In  coordinate  system  number  two.  The  coordinates  on  the  top 
and  right  In  figure  21  refer  to  this  coordinate  system. 

Study  of  the  first  example  In  Figure  20  shows  that  the  first  card 
generates  a wire  extending  In  the  Z-dlrectlon  from  0.05  to  0.0  meters. 

This  has  one  segment  and  a tag  of  two.  The  radius  Is  located  In  the  first 
entry  of  the  radius  table,  and  the  wire  Is  located  In  the  second  coor- 
dinate system. 

The  second  card  generates  a wire  extending  along  the  local  Y-axis 
from  0.0  to  0.05.  The  rest  of  the  values  In  this  command  are  defaulted 
to  those  specified  on  the  previous  card.  Note  that  all  subsections  lying 
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In  coordinate  system  number  two  will  In  the  sample  problem  be  given  a tag 
of  two.  This  will  aid  the  user  later  on  when  he  reviews  the  geometry  stored 
In  GEMACS. 


The  second  example  also  consists  of  two  cards.  The  first  generates 
a wire  along  the  global  Z-axIs  extending  from  0.1  to  0.0.  It  has  one 
subsection,  a tag  of  one  and  a radius  located  In  the  first  entry  of  the 
radius  table.  It  Is  located  In  coordinate  system  number  one.  Again,  the 
tag  number  Is  being  used  to  Indicate  In  what  coordinate  system  the  wire  lies. 

The  second  card  of  this  example  generates  a wire  along  the  global 
Y-axis  extending  from  0.0  to  0.1.  Its  remaining  parameters  are  Identical 
to  those  of  the  first  segment  In  this  example. 

This  method  of  subsection  generation  Is  a one-step  process.  Four 
commands  were  used  to  generate  the  four  subsections.  Specifying  points  and 
then  connecting  them  would  have  required  eight  commands. 

An  entire  geometry  could  be  generated  by  specifying  each  and  every 
subsection.  However,  for  the  example  this  would  require  more  than  700  such 
commands.  The  geometry  language  provides  other  commands  which  allow  one 
to  use  these  basic  subsections  as  building  blocks  to  generate  larger  struc- 
tures using  many  fewer  commands  than  700.  In  the  discussion  of  these  com- 
mands It  will  be  seen  that  some  are  used  quite  extensively  In  the  generation 
of  a wire  grid  model  for  the  top  of  the  hut. 

C.  Element  Generation 

It  was  shown  that  one  of  the  ways  of  grouping  subsections  having  some 
particular  characteristic  In  common  Is  by  giving  them  all  the  same  tag 
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number.  A second  technique  of  grouping  Is  by  Identifying  all  of  the 
subsections  as  part  of  some  "defined  element."  That  Is,  a specified 
group  of  subsections  Is  Identified  as  "WING"  or  "FINE"  or  COARSE"  or 
any  other  legal  name  chosen  to  Identify  that  set.  A legal  name  Is  an 
Identifier  composed  of  six  characters  or  less,  a character  being  one  of  the 
letters  of  the  alphabet  (A-Z)  or  one  of  the  numerals  (0-9).  However,  the 
first  character  must  be  one  of  the  letters  of  the  alphabet  (A-Z). 

The  generation  of  an  element  begins  with  a Define  Element  (OF)  command, 
shown  In  Figure  22.  As  can  be  seen,  It  Is  a simple  command  consisting  of 

the  mnemonic  DF  and  the  name  of  the  element.  Figure  22  also  shows  two 

examples  which  are  used  In  the  solution  of  the  sample  problem. 

The  define  element  command  Is  then  followed  by  a deck  of  geometry  com- 
mands which  may  Include  point  specification  and  connection,  wire  Input 
and  any  of  the  macrogenerators  to  be  discussed  very  shortly.  All  subsections 
and  points  defined  by  these  commands  will  become  associated  with  this 
element  and  as  such  may  be  referenced  by  calling  out  the  name  of  the  ele- 
ment. That  Is,  If  any  operation  Is  performed  on  the  element,  then  that 
operation  Is  performed  on  all  the  points  and  wires  specified  as  being  a 
part  of  that  element. 

When  all  of  the  comnands  that  are  needed  to  define  the  geometry  of  the 

element  have  been  Input  Into  GEMACS,  the  set  Is  then  closed  by  placing 

a Define  End  (DE)  command  after  the  last  geometry  command  of  the  element. 

This  command  Is  shown  In  Figure  23.  As  can  be  seen  It  simply  consists 
of  the  mnemonic  DE. 
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DE 

RA  0.000776 
DF  GRID1 

WR  0.0  0.0  0.05  0.0  0.0  0.0  1 2 1 2 

WR  0.0  0.0  0.0  0.0  0.05  0.0 

DE  $GRID1 


Figure  23.  DEFINE  END  COMMAND 


Note  however,  that  a comment  should  be  Included  to  Indicate  what  element 
Is  being  closed.  In  this  case  the  dollar  sign  ($)  Is  needed  as  well  as  the 
blank  after  the  mnemonic.  Figure  23  shows  an  example  of  the  definition 
of  an  element.  The  result  Is  the  pair  of  lines  shown  darkened  In  the  fine 
mesh  of  Figure  21.  The  first  card  sets  up  the  radius  table.  The  second 
card  defines  the  element  and  gives  It  the  name  GRID1.  The  third  and  fourth 
cards  are  the  same  as  those  used  In  the  first  example  In  Figure  20.  The  last 
card  In  Figure  23  Is  the  define  end  command  with  the  comment  that  this 
closes  the  group  to  be  Identified  as  GRID1 . 

Two  points  should  be  kept  In  mind  when  using  the  define  operation. 

First,  the  define  element  operation  may  be  nested.  This  Is  easiest  to  ex- 
plain by  using  an  example.  The  fictitious  example  In  Figure  24  starts  off 
with  a define  element  command,  labeling  the  element  as  "PLANE".  This  Is 
followed  by  a series  of  wire  Input  cards,  which  may  define  some  part  of  the 
aircraft.  These  cards  are  followed  by  another  define  element  card  with  a 
label  of  "WING".  The  following  wire  Input  cards  would  then  define  the 
wing.  Note  that  no  define  end  command  for  the  plane  has  as  yet  appeared. 

The  geometry  for  "WING"  thus  also  becomes  part  of  the  geometry  for  "PLANE". 
Therefore,  any  reference  to  the  label  "PLANE"  will  affect  all  the  sub- 
sections belonging  to  "WING"  also.  Any  reference  to  "WING"  will  only  affect 
the  second  set  of  wire  input  cards,  but  not  the  first  set,  which  are  a part 
of  "PLANE". 

The  first  define  end  card  applies  to  "WING",  and  It  closes  that  group 
only.  It  does  not  close  the  group  belonging  to  "PLANE".  Thus,  the  group 
"WING"  Is  a nested  defined  element.  Nesting  can  be  carried  to  ten  deep. 
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DE  $ PLANE 


Figure  24.  NESTED  DEFINED  ELEMENTS 
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This  Is  an  example  of  nesting  two-deep. 


Other  elements  could  be  defined  describing  the  tall,  the  vertical  stabi- 
lizer and  so  forth  In  the  same  manner.  These  would  again  be  nested  two- 
deep  with  the  group  "PLANE".  A reference  to  "TAIL"  would  not  In  any  way 
affect  the  group  "WING".  Once  a group  Is  closed,  no  new  elements  may 
be  added  to  that  group. 

The  final  define  end  command  would  close  the  group  "PLANE".  Any  refer- 
ence to  "PLANE"  would  therefore  affect  the  groups  "WING",  "TAIL",  and  any 
other  subelements  that  were  defined.  The  point  Is  that  closing  nested 
defined  elements  works  from  the  Innermost  set  outward.  Each  subelement 
Is  closed  before  the  element  to  which  It  belongs  Is  closed. 

The  second  point  Is  Illustrated  by  the  fictitious  example  shown  In 
Figure  25.  The  example  again  starts  out  by  defining  an  element  called 
"PLANE",  described  by  WR  commands.  Then,  an  operation  called  translation 
Is  performed  on  these  wires.  This  command  Is  explained  In  Section  IV. D. 

Briefly,  what  It  does  Is  generate  a duplicate  set  of  wires  using  the  first  set 
of  wire  Input  cards  as  a model.  Notice,  though,  that  up  to  this  point  the 
group  "PLANE"  has  not  been  closed  by  a define  end  command.  This  example 
shows  that  an  element  need  not  be  completely  defined  before  It  can  be 

! 

referenced  on  one  of  the  macrogenerator  commands.  Also,  keep  In  mind  that  the 

/ 

group  of  subsections  generated  by  the  translate  command  will  be  affected 
by  any  reference  to  the  element  "PLANE". 

Once  elements  have  been  defined,  various  operations  may  be  performed 
on  the  subsections  comprising  those  elements  by  simply  referring  to  the  name 
of  the  element.  Two  of  these  cownands  belong  to  the  class  of  element 


\ 
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generation  commands.  The  first  Is  the  Attach  (AT)  command.  Figure  26. 

This  command  results  In  a group  of  points  and/or  subsections  being  relo- 
cated at  the  origin  of  the  coordinate  system  specified  In  this  command. 

This  conmand  Is  useful  when  an  element  Is  defined  In  some  local  coordinate 
system;  and  It  will  be  used  on  several  different  elements,  each  of  which  Is 
defined  In  Its  own  Individual  coordinate  system.  Suppose,  for  example, 
that  the  shape  of  an  antenna  has  been  defined  In  coordinate  system  one. 

It  Is  to  be  located  on  the  left  and  right  wings  and  left  and  right  horizontal 
stabilizers  of  an  aircraft.  Suppose  further  that  each  of  the  four  aircraft 
sections  was  defined  In  Its  own  coordinate  system  labeled  two  to  five, 
respectively.  The  Implementing  commands  are  shown  In  the  example  In 
Figure  26. 

From  the  figure  It  can  be  seen  that  AT  Is  the  mnemonic  for  the  attach 
command.  Either  point  Identifiers  (PT),  tag  Identifiers  (TG)  or  defined 
elements  (DF)  may  be  referenced  by  specifying  the  number  of  the  point  or  tag 
or  the  name  of  the  element.  The  final  piece  of  data  Is  the  number  of  the 
coordinate  system  Into  which  the  group  of  points  or  wires  will  be  trans- 
ferred. It  cannot  be  defaulted  In  this  command. 

However,  It  must  be  remembered  that  the  origin  of  the  coordinate  system 
In  which  "ANTENA"  Is  defined  will  be  located  at  the  origin  of  the  coor- 
dinate system  specified  In  the  attach  command.  Depending  on  the  relative 
positions  of  the  respective  coordinate  axes,  there  may  result  both  trans- 
lation and  rotation.  Therefore,  the  origins  of  coordinate  systems  two 
through  five  must  be  at  the  locations  on  the  wings  and  stabilizers  at  which 
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Figure  26.  ATTACH  COMMAND 


the  antenna  Is  to  be  located. 


The  final  command  In  this  group  Is  the  Combine  Elements  (CE)  command 
In  Figure  27.  The  sole  purpose  of  this  command  Is  to  make  referencing 
elements  easier  for  the  user.  When  several  elements  such  as  "BODY",  "WING" 
and  "TAIL"  have  been  defined,  then  this  command  may  be  used  to  lump  their 
groups  under  the  element  name  "PLANE".  This  would  be  an  alternative  to 
what  was  done  In  the  example  In  Figure  24.  The  advantage  of  this  technique 
over  that  shown  In  the  example  In  Figure  24  Is  that  In  the  latter  technique 
all  the  elements  and  subelements  must  be  defined  In  the  same  coordinate 
system,  which  may  not  be  convenient  due  to  the  complexity  of  the  geometry. 

The  elements  combined  In  the  example  In  Figure  27  may  be  defined  In  their 
own  Individual  coordinate  systems. 

However,  once  the  elements  have  been  combined  under  one  name,  then 
the  Individual  elements  can  no  longer  be  accessed  under  their  Individual 
names.  In  the  example  just  cited,  once  the  element  "BODY"  has  been  Incor- 
porated under  the  name  "PLANE",  then  one  cannot  any  longer  refer  to  the 
Individual  element  "BODY".  However,  using  the  method  shown  In  Figure  24, 
the  Individual  elements  can  still  be  Individually  accessed. 

Looking  at  the  command  Itself  In  Figure  27,  It  Is  seen  that  CE  Is  the 
symbol  for  this  command.  Then,  NAME1  Is  the  name  by  which  the  resultant 
group  will  be  known.  'The  remaining  names  given  In  the  command  are  the  names 
of  the  Individual  defined  elements  which  are  to  be  combined.  In  the  example, 
"PLANE"  will  be  the  name  of  the  resultant  structure.  "BODY",  "WING"  and 
"TAIL"  will  be  the  Individually  defined  elements  making  up  this  structure. 

The  definition  of  elements  finds  Its  greatest  utilization  when  the 
geometry  of  the  structure  Is  being  generated  by  the  so-called  macrogener- 
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CE  PUNE  BODY  WING  TAIL  (FICTITIOUS) 


Figure  27,  COMBINE  ELEMENTS  COMMAND 


ators,  the  operations  of  reflection,  rotation  and  translation.  This  last 
comnand  Is  extremely  useful  when  a wire  grid  model  of  a structure  is  being 
generated. 


D.  Macroqenerators 

A macrogenerator  Is  a command  which  uses  as  Input  a specified  group  of 
already  defined  subsections  and  generates  one  or  more  Identical  groups  of 
subsections  by  either  reflection,  rotation  or  translation. 

Before  beginning  the  discussion  of  these  three  commands  the  numbering 
sequence  of  the  points  and  lines  should  be  described.  Remember  that  when 
points  are  defined  by  the  point  specification  (PT)  command,  the  user  speci- 
fies what  the  point  number  Is  to  be.  When  one  of  the  macrogenerator  commands 
Is  used,  the  generation  of  additional  points  occurs.  These  new  points  will 
be  ordered  in  the  same  sequence  In  which  the  original  points  were  specified. 
That  Is,  suppose  points  one  through  five  were  defined  and  then  reflected 
through  some  plane.  Point  six  will  be  the  Image  of  point  one,  point  seven 
will  result  from  point  two,  and  so  on;  point  ten  being  the  image  of  point 
five.  If  further  Individual  points  are  going  to  be  defined,  then  the  user 
must  remember  to  start  the  numbering  at  eleven,  since  point  numbers  must 
be  globally  unique. 

The  same  thing  will  happen  with  the  numbering  of  subsections.  Sub- 
sections are  automatically  ordered  in  the  sequence  In  which  they  are  gener- 
ated. This  ordering  was  Illustrated  In  the  example  of  Figure  18.  if 
the  points  In  the  previous  example  were  subsections,  then  subsection  one 
would  Image  as  subsection  six  with  the  same  relative  positions  of  negative 
and  positive  ends. 
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Subsections  and  points  In  the  plane  of  reflection  will  also  generate 
colocated  subsections  and  points  with  their  Individual  Identifiers.  Thus,  In 
the  example  If  point  or  subsection  three  were  In  the  plane  of  reflection. 

It  would  generate  point  or  subsection  eight,  respectively.  Since  GEMACS 
Is  generating  these  Identically  numbered  points,  It  will  flag  these  points 
and  Ignore  their  existence.  The  Identical  segments  will  be  zeroed  out 
In  all  calculations  Involving  the  Interaction  matrix. 

The  first  macrogenerator  command  to  be  discussed  Is  the  Reflect  (RF) 
command,  which  Is  shown  In  Figure  28.  This  command  executes  the  symmetry 
operation  of  reflecting  through  a plane  normal  to  the  axis  specified.  For 
example.  If  the  X-axis  Is  specified,  then  the  plane  of  reflection  Is  the  Y-Z 
plane;  and  all  (x,y,z)  coordinates  become  (-x,y,z)  coordinates.  Note  that 
since  segments  may  be  located  In  the  plane  of  reflection,  when  a wire  Is 
located  In  the  Y-Z  plane,  then  a wire  would  be  generated  In  the  Y-Z  plane 
and  be  superimposed  on  the  original  wire.  However,  Its  existence  would  be 
flagged  as  a duplicate  wire;  and  It  would  be  Ignored  In  all  future  calcula- 
tions. The  same  Is  of  course  true  for  points  located  In  the  plane  of 
reflection. 

As  can  be  seen  from  Figure  28,  the  mnemonic  for  this  command  Is  RF. 

The  points,  subsections  or  defined  elements  to  be  reflected  are  referenced 
by  their  Identifier  for  points  or  by  the  tag  for  subsections  or  by  their 
name  for  defined  elements.  The  symbol  A1  represents  either  the  X-axis, 

Y-axis  or  Z-axIs  depending  on  which  plane  the  reflection  Is  to  occur 
through  first.  If  the  X-axis  Is  specified,  then  reflection  will  be  through 
the  Y-Z  plane.  If  the  Y-axis  Is  specified,  reflection  will  occur  through 
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Figure  28.  REFLECT  COMMAND 
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the  X-Z  plane.  If  the  Z-axIs  Is  specified,  reflection  will  occur  through 
the  X-Y  plane. 

! 

The  symbol  A2  represents  the  axis  normal  to  the  plane  through  which 
the  second  reflection  will  occur.  Similarly,  A3  represents  the  axis 
normal  to  the  plane  through  which  the  third  reflection  will  occur.  For 
any  one  of  these  symbols  the  Input  would  be  X,  Y or  Z.  No  axis  designation 
shall  occur  more  than  once.  The  order  In  which  the  axis  designations  appear 
will  determine  the  order  of  the  planes  through  which  reflection  will  occur. 
Finally,  If  the  reflection  operation  Is  being  performed  on  an  element  being 
defined  (that  Is,  a reference  to  an  Incompletely  defined  element  - 
Figure  25),  all  subsections  generated  by  the  command  will  be  associated  with 
the  element  being  defined. 

It  should  be  noted  that  although  the  number  of  planes  of  reflection 
may  vary  from  one  to  three,  none  of  the  symbols  Al,  A2  or  A3  Is  shown  as 
one  that  can  be  defaulted.  The  code  Is  set  up  such  that  only  one  or  two 
axes  can  be  Input  without  having  to  worry  about  Including  two  or  one  default 

: 

values,  respectively.  This  Is  shown  In  the  example  given  In  Figure  28. 

The  symbol  INCTAG  represents  the  tag  Increment  parameter.  The  Integer 
input  here  will  be  added  to  the  tag  numbers  specified  on  the  subsections 
that  are  to  be  reflected.  Subsections  generated  by  the  first  reflection 
will  have  tags  Incremented  by  the  value  of  INCTAG.  Subsections  generated 
by  the  second  reflection  will  have  tags  Incremented  by  twice  the  value  of 
INCTAG  over  the  tag  numbers  of  the  original  subsections.  Subsections  gener- 
ated by  the  third  reflection  will  have  tag  Identifiers  Incremented  by 
four  times  the  value  of  INCTAG  over  the  tags  of  the  original  subsections. 
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Note  that  the  parameter  INCTAG  may  be  defaulted.  The  default  value 
Is  always  zero.  This  means  that  the  tag  Identifiers  of  subsections  gener- 
ated by  any  order  of  reflection  will  be  equal  to  the  tags  of  the  original 
subsections  being  reflected.  Also,  If  only  a point  Is  being  reflected, 
then  INCTAG  will  always  be  defaulted  since  a tag  Identifier  only  applies  to 
subsections  or  to  subsections  contained  within  a defined  element. 

The  final  parameter  In  this  comnand  Is  the  number  of  the  coordinate 
system  In  which  the  reflection  Is  to  take  place.  Note  that  It  can  be 
defaulted;  and  that  It  has  the  usual  default  value,  which  Is  the  last  Identi- 
fier for  a coordinate  system  specified  on  any  previous  command. 

An  example  of  this  command  Is  shown  In  Figure  28.  Figure  29  (a)  shows 
the  original  configuration.  Figure  29(b)  shows  the  results  of  the  first 
reflection,  and  Figure  29(c)  shows  the  final  configuration  that  results 
from  the  command. 

In  Figure  29  the  Z-axIs  Is  coming  out  of  the  plane  of  the  diagram. 

There  Is  a single  wire  In  the  X-Y  plane  with  a tag  of  two,  and  It  Is  assumed 
that  It  Is  the  only  wire  with  this  tag  number.  This  wire  Is  subdivided 
Into  two  subsections,  which  are  sequenced  as  shown.  The  comnand  In  Figure 
28  tells  GEMACS  to  first  reflect  the  wire  with  a tag  of  two  through  the 
Y-Z  plane  (normal  to  the  X-axis)  and  to  Increase  the  tag  by  two. 

Figure  29(b)  shows  the  results  of  this  operation.  The  original  wire 
Is  left  unaffected.  A second  wire  has  been  generated  In  the  X-Y  plane. 

The  wire  has  a tag  of  four  and  Is  subsectioned  Into  two  segments,  subsection 
three  being  the  Image  of  subsection  one  and  segment  four  being  generated 
by  segment  two. 


The  example  In  Figure  28  also  specifies  that  a second  reflection  Is  to 
be  performed.  This  Is  to  be  through  the  X-Z  plane;  that  Is,  normal  to  the 
Y-axis.  The  tag  In  this  case  Is  to  be  Incremented  by  two  times  two,  resul- 
ting from  the  fact  that  this  Is  the  second  reflection  Incrementing  the  tag 
Identifier  by  two. 

Figure  29(c)  shows  the  results  of  this  operation.  First  of  all,  note 
that  not  only  were  the  original  subsections  reflected  but  that  the  subsec- 
tions generated  as  a result  of  the  first  reflection  were  also  reflected. 

The  wire  with  the  tag  of  two  generated  a wire  that  has  a tag  of  six  (two 
times  two  plus  two),  which  has  two  segments.  Subsections  five  and  six 
were  generated  by  segments  one  and  two,  respectively.  The  wire  with 
a tag  of  eight  (two  times  two  plus  four)  Is  the  Image  of  the  wire  with 
the  tag  of  four.  Correspondingly,  segments  seven  and  eight  are  the  Images 
of  subsections  three  and  four,  respectively.  This  shows  how  the  subsection 
numbers  are  sequenced  as  new  subsections  are  generated  and  how  tag  Identifiers 
are  Incremented. 

The  second  macrogenerator  command  Is  the  Rotation  (RT)  command. 

The  general  form  of  the  conmand  Is  shown  In  Figure  30.  This  contnand  Is 
extremely  useful  to  generate  geometries  that  have  one  or  more  axes  of 
revolution,  such  as  a cylinder  or  a sphere.  The  mnemonic  for  this  operation 
Is  RT.  Points,  line  segments  or  defined  elements  may  be  rotated.  Identi- 
fying them  by  their  number,  tag  or  name,  respectively.  The  value  of  IADD 
Indicates  the  number  of  additional  sets  of  points  or  subsections  that  are 
to  be  generated. 
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Figure  30.  ROTATE  COMMAND 


The  parameters  RX,  RY,  RZ  are  the  total  angles  through  which  the  original 
points  and  wires  are  to  be  rotated.  Keep  In  mind  that  these  do  not  repre- 
sent Incremental  changes  around  the  axes.  For  example.  If  IADD  were  equal 
to  three  and  RX  equal  to  fifteen,  then  each  original  point  or  subsection 
would  generate  a new  point  or  subsection  at  angles  of  five,  ten  and  fifteen 
degrees  of  rotation.  The  result  would  not  be  points  and  subsections  located 
at  angles  of  15,  30  and  45  degrees. 

One  other  point  to  keep  in  mind  Is  that  none  of  these  parameters 
may  be  defaulted  or  left  out  as  was  the  case  with  the  reflect  command. 

Also,  the  order  In  which  the  parameters  appear  Is  Important.  The  rotation 
around  the  X-axis  (RX)  must  be  first,  followed  by  the  rotation  about  the 
Y-axis  (RY),  which  In  turn  Is  followed  by  the  rotation  about  the  Z-axIs 
(RZ).  Therefore,  If  rotation  about  the  Y-axis  only  Is  desired,  the  value 
for  RY  must  be  preceded  by  a zero  for  RX  and  followed  by  a zero  for  RZ. 

One  final  note  on  the  direction  of  positive  rotation  for  RX,  RY  and  RZ. 
Starting  at  the  original  point  or  subsection,  positive  rotation  Is  consi- 
dered to  be  counterclockwise  around  the  X-axis  for  RX,  around  the  Y-axis  for 
RY  and  around  the  Z-axIs  for  RZ.  For  example.  If  a point  exists  on  the  posi- 
tive X-axis  and  Is  to  be  rotated  around  the  Z-axIs,  then  the  direction 
of  positive  rotation  Is  from  the  positive  X-axis  to  the  positive  Y-axis. 

The  final  two  parameters  are  the  same  as  those  that  appear  In  the  re- 
flect coimtand,  namely  the  tag  Increment  and  the  number  of  the  coordinate 
system  In  which  the  rotation  Is  to  take  place.  The  default  values  are  also 
the  same;  zero  for  the  tag  Increment  and  the  last  coordinate  system  speci- 
fied on  ar\y  previous  command  for  NCS. 
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Figure  30  also  Illustrates  the  use  of  this  command.  There  Is  originally 
a point  located  on  the  Y-axis.  It  Is  denoted  by  the  "X"  and  has  a 
number  one.  The  sample  command  says  that  this  point  Is  to  be  rotated  three 
hundred  degrees  around  the  Z-axIs,  which  Is  directed  out  of  the  figure. 

During  this  rotation  five  new  points  are  to  be  generated.  The  diagram 
shows  the  result.  The  result  Is  six  points  spaced  around  a circle,  each 
separated  from  its  neighbors  by  sixty  degrees.  The  direction  of  positive 
rotation  Is  Indicated  by  the  arrow  and  determines  the  order  In  which 
the  points  are  finally  numbered. 

As  Is  the  case  with  the  reflect  command,  If  an  element  Is  being 
defined  when  the  rotate  command  Is  encountered,  all  the  new  points  and 
subsections  will  be  associated  with  that  particular  element  being  defined. 

The  last  macrogenerator  Is  the  Translation  (XL)  command,  whose  general 
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format  Is  shown  in  Figure  31.  In  form  It  Is  almost  Identical  to  the  rotation 
command  that  has  just  been  discussed.  The  result  of  this  operation  Is  the 
generation  of  additional  points  and  subsections  that  are  Identical  to  those 
specified  but  are  displaced  by  a specified  distance. 

The  symbol  for  this  command  Is  XL.  Again  the  operation  can  be  performed 

l 

on  points,  subsections  and  defined  elements  by  specifying  the  number,  the  tag 
or  the  name,  respectively.  The  number  of  new  sets  of  points  or  segments 
to  be  generated  Is  given  by  the  value  of  IADD.  If  IADD  Is  zero,  then  the 
original  set  Is  translated  to  the  new  location. 

The  parameters  DX,  DY  and  DZ  are  the  Incremental  changes  In  the  X,  Y 
and  Z directions,  respectively,  that  will  be  added  to  each  succeeding  new 
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set  of  points  or  subsections  that  have  been  generated.  For  example.  If 
point  one  Is  located  at  X equals  five  and  If  three  new  points  are  to  be 
generated  at  Incremental  spaclngs  of  three  along  the  X-axis,  then  point 
two  will  be  located  at  X equals  eight,  point  three  will  be  located  at 
X equals  eleven;  and  point  four  will  be  located  at  X equals  fourteen. 

Thus,  (DX.DY.DZ)  forms  an  Incremental  translation  vector. 

The  last  two  parameters,  INCTAG  and  NCS,  are  the  same  ones  that  have 
been  discussed  for  the  reflect  and  rotate  commands.  They  have  the  same  de- 
fault values  of  zero  and  the  last  specified  coordinate  system. 

As  Is  the  case  with  the  first  two  macrogenerator  commands.  If  an  ele- 
ment Is  in  the  process  of  being  defined  when  the  translate  comnand  Is  en- 
countered, all  the  new  points  and  wires  will  be  associated  with  the  element 
being  defined. 

This  command  finds  a tremendous  amount  of  use  In  the  generation  of 
a constant-mesh  wire  grid.  The  example  shown  In  Figure  31  is  taken  from 
the  input  deck  for  the  analysis  of  the  sample  problem.  The  results  of  the 
translation  are  shown  in  Figure  32  as  the  darkened  lines  between  Y*  equals 
0.05  and  0.7  meters  in  the  second  coordinate  system.  The  first  five  cards 
generate  subsections  one  and  two  and  are  the  same  as  those  in  Figure  23.  The 
defined  element  GRID1  is  shown  in  Figure  21  as  the  darkened  pair  of  lines 
In  the  fine-mesh  wire  grid.  The  translate  command  then  says  to  take  the 
defined  element  GRID1  and  generate  13  additional  elements,  each  translated 
from  the  preceding  one  by  0.05  meters  in  the  Y' -direction.  The  value  for 
INCTAG  has  been  defaulted  so  that  all  of  the  new  subsections  will  have  a 
tag  of  two,  which  is  the  tag  of  the  original  subsections.  The  value  for 
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the  number  of  the  coordinate  system  Is  also  defaulted  since  the  new  elements 
are  In  coordinate  system  number  two,  which  Is  specified  In  the  first  wire 
input  (UR)  command. 


' 
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If  one  follows  through  the  numbering  of  the  subsections.  It  Is  found 
that  the  last  pair  of  newly  generated  subsections  Is  numbered  27  and  28. 

As  will  be  shown  later  In  the  discussion  of  the  actual  Input  deck,  the 
addition  of  five  more  cards  will  completely  generate  the  fine-mesh  wire 
grid,  once  the  coordinate  systems  have  been  defined.  This  shows  the  power 
that  Is  afforded  by  the  macrogenerator  commands. 

E.  Renumbering 

When  the  geometry  Is  generated  the  command  stream  Is  usually  struc- 
tured In  order  to  make  the  job  as  easy  as  possible.  Little  thought  Is 
normally  given  to  what  subsections  have  what  segment  numbers  attached  to  them. 
Consequently,  the  numbering  that  results  may  not  be  such  as  to  make  full 
use  of  the  efficiency  of  the  banded  matrix  Iteration  (BMI)  scheme.  Figure 
33  shows  a simple  example  of  this. 

Suppose  that  there  are  two  antennas,  each  of  which  Is  one-quarter  of 
a wavelength  long.  They  are  parallel  to  each  other  and  separated  by 
1/100  of  a wavelength.  Each  antenna  can  be  generated  by  a single  wire  Input 
(WR)  command,  specifying  that  each  antenna  be  segmented  into  five  subsec- 
tions. Thus,  each  subsection  will  be  one-twentieth  of  a wavelength  long. 

This  Is  shown  In  figure  33(a).  This  means  that  the  center-to-center  separa- 
tions between  parallel  subsections  Is  five  times  smaller  than  the  center-to- 
center  seoaratlon  between  adjacent  subsections  on  the  same  antenna.  The 
interactions  between  subsections  on  separate  antennas  are  greater  than  the 


«« 


5 

4 

3 

2 


10 

9 

8 

7 

6 


(a)  ORIGINAL  NUMBERING 


9 

7 

5 

3 


10 

8 

6 

4 

2 


(fictitious) 


(b)  RENUMBERED 

NOTES:  1)  Antenna  Separation  = 0.01a 
2)  Segment  Lengths  = 0.05* 

Figure  33.  EXAMPLE  OF  RENUMBERING  (Not  to  Scale) 


Interactions  between  subsections  on  the  same  antenna. 

When  the  BMI  technique  Is  used.  Its  efficiency  Is  Increased  by 
having  as  many  of  the  largest  Interactions  as  possible  clustered  near  the 
main  diagonal  of  the  Interaction  matrix.  In  this  example,  the  subsections 
should  be  numbered  as  shown  In  Figure  33(b).  In  the  scheme  shown  In  Figure 
33(a)  the  bandwidth  would  need  to  be  at  least  five  In  order  to  Include  the 
large  Interaction  between  subsections  one  and  six.  The  numbering  In  Figure 
33(b)  would  require  a bandwidth  of  three  to  Include  all  the  major  Inter-  . 
actions;  that  Is,  those  between  subsections  one  and  two,  one  and  three, 
and  one  and  four. 

Figure  34  shows  the  forms  of  the  Renumber  (RN)  command.  The  symbol  RN 
denotes  that  this  Is  the  renumber  command.  The  symbols  II  through  In  are 
the  new  segment  numbers  that  are  to  be  associated  with  the  segment  In 
that  position  In  the  list.  That  Is,  11  located  In  position  one  of  the 
list  Is  the  Integer  Identifier  that  will  be  associated  with  what  Is  origi- 
nally subsection  one.  In  general,  position  k corresponds  to  the  kth 
subsection  generated  since  the  beginning  of  the  geometry  Input  or  since 
the  last  renumber  command.  The  Integer  located  In  this  position  is  the  new 
segnent  number  that  will  be  associated  with  the  kth  subsection  generated. 

Three  things  must  be  kept  In  mind  when  using  this  conmand.  First,  no 
two  nunbers  In  this  list  may  ever  be  the  same.  This  will  give  the  same 
subsection  Identifier  to  two  different  subsections.  Second,  the  list 
may  not  exceed  256  Items  due  to  the  structure  of  the  code.  This  does  not 
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Figure  34.  RENUMBER  COMMAND 
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mean  that  all  the  numbers  Input  be  between  one  and  256  or  that  they  follow 
any  particular  sequence.  However,  If  there  are  N subsections  and  all  of 
them  are  to  be  renumbered,  then  all  of  the  Integers  one  through  N must  ap- 
pear In  a renumber  command.  Finally,  when  a renumber  command  Is  encountered, 
all  of  the  subsections  generated  since  the  last  renumber  command  are  renum- 
bered according  to  the  sequence  specified  by  the  Integers  In  the  list.  If 
this  Is  the  first  renumber  conmand,  then  all  of  the  subsections  generated 
since  the  beginning  of  the  geometry  Input  will  be  resequenced.  Therefore, 

If  the  first  form  of  the  coirmand  Is  used,  no  more  than  256  subsections 
should  be  generated  before  a renumber  command  Is  encountered. 

A shortcut  Is  available  and  Is  shown  as  the  second  form  In  Figure  34. 

When  a negative  Integer  Is  encountered  In  the  list  such  as  -Im  In  Figure 
34,  the  new  sequence  number  of  the  subsection  corresponding  to  that  position 
In  the  list  will  be  the  absolute  value  of  the  negative  Integer;  and  the  new 
sequence  numbers  will  be  Incremented  by  one  until  the  number  of  segments 
less  one  Identified  by  the  next  Integer  In  the  list  (Ip  In  Figure  34)  has 
been  counted.  Figure  35  shows  an  example  of  this.  Suppose  that  subsections 
one  through  ten  have  been  generated  and  that  they  are  to  be  renumbered. 

These  are  represented  by  the  line  labeled  "OLD"  In  Figure  35.  By  employ- 
ing the  conmand  shown  in  the  figure  one  ends  up  with  the  sequence  shown  In 
the  line  labeled  "NEW". 

Item  one  In  the  list  has  the  value  two.  Thus,  what  was  subsection  one 
Is  now  subsection  two.  Four  Is  the  value  of  Item  two  In  the  list,  and  the 
result  Is  that  subsection  two  now  becomes  subsection  four.  Item  three  Is  a 
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Figure  35.  RENUMBER  EXAMPLE 
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negative  five.  Subsection  three  becomes  subsection  five.  The  next  Item  In 
the  list  Is  four.  That  means  that  the  next  three  subsections  are  sequen- 
tially numbered,  starting  with  six  since  the  preceding  subsection  Is  num- 
bered five.  Therefore,  segments  four,  five  and  six  become  subsections  six, 
seven  and  eight,  respectively.  The  last  four  subsections  are  renumbered  In 
the  same  manner  as  the  first  two  subsections. 

The  example  shown  In  Figure  34  is  the  renumber  command  that 
would  result  in  the  segment  numbering  shown  in  Figure  33. 

This  completes  the  discussion  of  the  GEMACS  geometry  language.  It  Is 
seen  to  be  a very  simple  language  and  yet  It  Is  sufficiently  general  so  that 
one  can  develop  any  geometry  to  any  degree  of  accuracy.  For  example,  the  model 
of  the  hut  and  antennas  requires  714  subsections;  but  only  51  commands 
are  needed.  How  this  is  accomplished  is  shown  In  the  next  subsection. 

F.  Sample  Problem  Geometry  Input  Deck 

The  first  thing  that  must  be  done  is  to  set  up  the  global  parameters 
for  the  problem.  This  is  done  by  using  the  general  commands  discussed  In 
section  IV-A.  This  part  of  the  deck  Is  shown  In  Figure  36.  Line  one  sets 
up  all  the  radii  that  will  be  used  In  the  analysis  of  the  coupling.  In  this 
case  all  the  subsections  will  have  the  same  radius,  that  Is,  the  radius  of 
the  physical  antenna  wires. 

Figure  12  shows  the  total  wire  grid  that  will  be  generated.  It  also 
shows  the  two  coordinate  systems  that  will  be  employed.  In  both  cases  the 
X-axis  comes  out  of  the  figure.  Coordinate  system  one  Is  Identical  to  the 
global  system  and  has  axes  Y and  Z.  Coordinate  system  two  Is  the  system  In 
the  darkened  box  and  has  axes  Y*  and  Z'.  Its  origin  is  displaced  from  the 
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global  origin  by  being  located  In  the  global  system  at  Y equals  0.3  and  Z 
equals  0.2  meters.  Lines  two  and  three  In  Figure  36  provide  this  Informa- 
tion to  GEMACS . 

Figure  37  shows  the  cards  needed  to  generate  the  fine-mesh  wire  grid. 

Much  of  this  has  already  been  discussed  In  connection  with  Figure  31. 

Card  one  specifies  an  element  which  Is  one  row  of  the  fine-mesh  wire  grid. 
Cards  two  through  five  generate  segments  one  and  two  In  Figure  32  and  have 
been  discussed  with  reference  to  Figures  20  and  21.  Card  six  generates  the 
darkened  lines  In  Figure  32.  Card  seven  generates  subsection  29  In  Figure 
38. 

Figure  38  shows  the  results  of  the  first  eight  cards  shown  In  Figure 
37.  A basic  building  block  has  been  generated  which  Is  translated  by  the 
coirmand  given  In  card  nine  a total  of  nine  times.  The  result  of  this  trans- 
lation Is  shown  In  Figure  39  as  subsections  30  to  290.  Card  ten  In  Figure 
37  closes  the  grid  by  drawing  a line  across  the  top  of  the  grid  which  Is 
segmented  Into  14  subsections.  This  generates  segments  291  to  304.  This 
completes  the  generation  of  the  fine-mesh  wire  grid  model  of  the  top  of  the 
hut  In  the  vicinity  of  the  two  antennas  of  Interest. 

Figure  40  shows  the  cards  that  are  needed  to  generate  the  models  of  the 
antennas.  The  diagram  In  the  figure  shows  the  results  of  the  operations 
along  with  the  subsection  numbering.  The  first  card  generates  a model  of  the 
transmitting  antenna.  It  has  a length  of  one-tenth  of  a meter  and  Is  divided 
up  Into  two  subsections.  It  has  been  given  a tag  Identifier  of  three  to  flag 
these  segments  In  the  segment  table  output  by  GEMACS.  The  second  card 
generates  a model  of  the  receiving  antenna,  which  has  a length  of  three- 
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tenths  of  a meter,  six  subsections  and  a tag  Identifier  of  four.  It  Is  seen 
from  figure  40  that  subsections  305  and  306  are  on  the  transmitting  antenna 
and  subsections  307  to  312  are  on  the  receiving  antenna. 

The  third  card  shown  In  Figure  40  Is  a renumber  command  that  Is  used 
to  renumber  all  of  the  segments  that  have  thus  far  been  generated.  It  Is 
used  to  renumber  the  antennas  so  that  their  segment  numbers  are  more  In  line 
with  the  subsection  numbers  of  the  nearby  wire  mesh  segments.  The  result 
Is  shown  In  the  figure  In  the  Appendix.  The  first  155  subsections  retain 
their  original  numbers.  Then  the  wire  mesh  segments  numbered  156  to  163 
are  renumbered  as  162  to  169.  The  remaining  141  subsections  originally 
numbered  164  to  304,  are  renumbered  as  172  to  312.  The  two  transmitting 
antenna  subsections  are  then  renumbered  as  170  and  171,  as  shown  In  Figure 
40.  Finally,  the  subsections  on  the  receiving  antenna  are  renumbered  as  shown 
In  the  parentheses  In  Figure  40.  The  final  numbering  scheme  for  the  fine- 
mesh  wire  grid  Is  as  shown  In  the  figure  In  the  Appendix. 

This  takes  care  of  all  the  modeling  that  will  be  done  In  coordinate 
system  two.  The  coarse-mesh  wire  grid  will  now  be  generated,  and  will  be 
defined  in  coordinate  system  number  one.  This  development  will  pretty 
much  follow  the  basic  lines  used  in  the  generation  of  the  fine-mesh  wire 
grid. 

It  Is  necessary  to  first  define  one  of  the  basic  building  blocks  that 
will  be  used  over  and  over.  Figure  41  outlines  these  basic  building  blocks 
by  the  darkened  lines.  The  numbers  are  the  subsection  numbers  that  will 
eventually  be  used  as  Identifiers  for  the  segments  making  up  these  building 
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blocks.  The  block  of  Interest  at  this  point  Is  GRID2  and  Is  made  up  of 
segments  313  and  314. 

Figure  42  shows  the  cards  that  are  needed  to  generate  GRID2  and  one  of 

the  rows  that  will  be  used  to  develop  most  of  the  coarse-mesh  grid.  This 

element  Is  called  R0W2  and  Is  shown  In  Figure  43  as  the  darkened  lines.  In 
Figure  42  cards  one  and  two  introduce  the  definitions  of  R0W2  and  GRID2, 
respectively,  the  latter  being  one  of  the  basic  building  blocks  shown  In 
figure  41.  Cards  three,  four  and  five  complete  the  definition  of  GRID2, 
which  is  made  up  of  subsections  313  and  314.  Card  six  translates  these 
segments  twice,  generating  segments  315  to  318  (Figure  43). 

At  this  point  the  second  of  the  basic  building  blocks,  GRID3,  shown 
as  subsections  319  and  320  In  Figure  41,  must  be  defined.  This  is  done  In 

cards  seven  through  ten  In  a manner  Identical  to  the  generation  of  GRID2. 

This  block  Is  then  translated  by  the  command  In  card  11  thirteen  times, 
generating  subsections  321  to  346,  ending  at  Y equals  1.0.  Cards  12  and 
13  again  translate  GRID2,  each  one  time,  to  generate  segments  347  to  350. 

The  reason  two  commands  are  needed  Is  because  the  displacement  distance 
along  the  Y-axis  changes  for  each  translation. 

The  element  R0W2  Is  completed  by  defining  subsection  351,  which  closes 
off  the  right-hand  side  of  the  element  R0W2.  This  Is  done  In  card  14,  while 
card  15  closes  the  set  of  cards  belonging  to  the  element  R0W2  (Figure  43). 

The  element  R0U2  Is  then  translated  along  the  Z-axIs  one  time  to  gener- 
ate the  second  row  of  elements  In  the  coarse-mesh  wire  grid.  For  this  row 
the  subsection  numbers  extend  from  352  to  390,  as  shown  In  Figure  43. 
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Figure  42.  GENERATION  OF  ELEMENT  R0W2 


4-68 


At  this  point  It  Is  time  to  generate  the  third  of  the  basic  building 
blocks,  GRID4,  shown  as  subsection  391  and  392  In  Figure  41,  and  another 
entire  row  which  will  be  translated  to  generate  the  middle  section  of  the 
coarse  grid.  The  set  of  commands  to  do  this  Is  shown  In  Figure  44.  The 
results  of  the  generation  of  the  element  R0W3  are  shown  In  Figure  45. 

In  Figure  44  cards  one  and  two  set  up  the  definitions  of  the  elements 
R0W3  and  GRID4.  Cards  three,  four  and  five  complete  the  generation  of  GRI04 
resulting  In  segments  391  and  392  In  Figure  41.  Cards  six,  seven  and  eight 
translate  the  element  GRI04  to  generate  subsections  393  to  400.  Finally, 
card  nine  results  In  the  generation  of  subsection  401  which  closes  off  the 
right-hand  side  of  the  element  R0W3. 

Note  that  subsection  29  (Figure  38)  of  the  fine-mesh  grid  and  subsec- 
tion 397  (Figure  45)  of  the  coarse-mesh  grid  are  really  the  same  wire 
segment.  GEMACS  will  recognize  this  fact  and  null  out  the  latter  subsection 
In  all  further  computations. 

Returning  to  Figure  44,  card  ten  completes  the  definition  of  the  ele- 
ment R0W3.  Card  11  translates  this  element  along  the  Z-axIs  to  generate 
the  central  portion  of  the  coarse-mesh  wire  grid.  The  translation  generates 
subsections  402  to  500  as  shown  In  Figure  46.  Overlapping  wires  will 
again  be  generated;  but  as  explained  previously,  these  will  be  nulled  out. 

All  that  remains  to  be  done  Is  to  develop  the  coarse-mesh  wire  grid 
for  values  of  Z greater  than  0.7  meter.  This  Is  accomplished  most  simply 
with  the  deck  of  cards  shown  In  figure  47.  The  result  Is  shown  In  Figure  48. 
First,  define  another  element,  R0W4,  which  Is  Identical  to  the  element  R0W2. 
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DF 

R0W3 

2 

DF 

GRID4 

3 

WR 

0.0  0.0 

0.25  0.0  0.0  0.2 

4 

WR 

0.0  0.0 

0.2  0.0  0.1  0.2 

5 

DE 

SGRID4 

6 

XL 

DF  GRID4 

2 0.0  0.1  0.0 

7 

XL 

DF  GRID4 

1 0.0  1.0  0.0 

8 

XL 

DF  GRID4 

1 0.0  1.1  0.0 

9 

WR 

0.0  1.2 

0.25  0.0  1.2  0.2 

10 

DE 

$R0W3 

11 

XL 

DF  R0W3 

9 0.0  0.0  0.05 

Figure  44.  GENERATION  OF  MIDDLE  SECTION  OF  COARSE  GRID 
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DF 

R0W4 
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XL 
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R0W2 
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0.0 

0.7 
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DE 

$R0W4 
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XL 

DF 

R0W4 

4 0.0 

0.0 

0.1 

5 

WR 

0.0 

0.0 

1.2  0.0 

0.3 

1.2 

3 

6 

WR 

0.0 

0.3 

1.2  0.0 

1.0 

1.2 

14 

7 

WR 

0.0 

1.0 

1.2  0.0 

1.2 

1.2 

2 

8 

END  OF  GEOMETRY  HUTTOP 

Figure  47.  COMPLETION  OF  COARSE-MESH  WIRE  GRID 


FIGURE  48.  FINAL  SUBSECTIONS  OF 
COARSE -MESH  WIRE  GRID 
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This  Is  accomplished  In  cards  one,  two  and  three.  The  element  R0W4  Is  made 
up  of  subsections  501  to  539.  Again,  any  overlapping  subsections  between 
the  fine-mesh  and  coarse-mesh  wire  grids  will  be  nulled  out.  This  element 
then  gets  translated  four  times,  generating  subsections  540  to  695.  The 
final  edge  of  the  hut  Is  modeled  by  cards  five,  six  and  seven  as  a wire 
with  19  subsections,  numbered  696  to  714. 

Subsections  501  to  695  could  have  been  generated  by  just  translating 
the  element  R0W2  Instead  of  generating  another  element.  But  that  would  have 
required  four  translation  conmands.  For  larger  problems  the  technique  used 
here  Is  the  more  efficient  way. 

This  then  completes  the  card  deck  needed  to  model  the  problem  given 
In  section  I.  Card  eight  In  Figure  47  signifies  this  to  GEMACS.  The  final 
numbering  on  each  of  the  subsections  for  the  top  of  the  hut  Is  shown  In  the 
figure  In  the  Appendix.  The  final  numbering  of  the  antennas  Is  shown  by 
the  numbers  In  parenthesis  In  Figure  40. 

Now  it  is  time  to  take  a look  at  the  task  command  language  to  see  what 
the  commands  mean,  how  they  are  constructed  and  how  they  are  put  together 
to  solve  a problem. 


V.  GEMACS  COMMAND  LANGUAGE 


i 


; 


The  GEMACS  command  language  Is  a free-fleld,  keyword-oriented  Input 
stream.  The  order  of  the  Items  Is  not  Important,  and  the  Items  In  each 
command  are  delimited  by  a blank  or  a comma.  An  Item  Is  considered  to 
be  all  of  the  Input  associated  with  a particular  parameter,  such  as 
THETA  = 90.  Each  Item  may  consist  of  several  entries,  each  entry  being 
called  a field.  Thus,  THETA  » 90  has  three  fields,  the  variable  THETA, 
the  equal  sign  and  the  value  90.  Blanks  may  be  Imbedded  between  the 
fields,  but  cannot  appear  within  a field.  If  THETA  were  Input  as  THE  TA, 
the  code  would  set  the  error  flag  since  two  symbols  were  found  when  one 
was  expected.  The  error  flag  would  Inhibit  execution  of  the  analysis. 

All  of  the  commands  are  read  prior  to  the  execution  of  any  command. 
Errors  are  noted  as  they  occur,  and  the  existence  of  an  error  does  not 
prevent  the  reading  of  further  commands.  All  of  the  Input  deck  Is 
scanned  for  errors  during  the  submittal  of  the  code  for  execution. 

There  are  several  conventions  that  should  be  kept  In  mind  when 
any  of  the  commands  are  being  set  up.  These  basically  are  Involved  with 
comment  cards,  comments  on  cards  and  continuation  cards.  Comment  cards 
are  those  cards  that  are  Inserted  within  a deck  to  Identify  the  deck  or 
to  explain  the  purpose  of  commands  or  Items  on  commands.  They  are  non- 
executable statements  and  are  Ignored  by  GEMACS  once  they  have  been 
output  back  to  the  user  In  the  listing  of  Input  deck.  They  are  Identified 
as  comment  statements  by  the  presence  of  a dollar  sign  ($)  as  the  first 
nonblank  character  on  the  card. 
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In  a similar  manner  short  comments  may  be  appended  right  In  the 
conmand  Itself  Instead  of  appearing  as  a separate  card.  This  Is  done 
by  placing  a space,  a dollar  sign  ($)  and  the  conment  after  the 
last  Item  In  the  command.  Again,  GEMACS  does  nothing  with  this 
portion  of  the  command  other  than  to  list  It  In  the  listing  of  the 
Input  deck. 

The  continuation  card  of  a command,  or  more  than  one  line  of  Input,  Is 
signified  by  the  presence  of  a continuation  character  in  the  first  column 
of  each  succeeding  Input  card  for  this  command.  This  continuation  character 
Is  an  asterisk  (*),  and  again  it  must  be  in  the  first  column  of  the  card. 

These  are  the  only  basic,  general  conventions  that  must  be  observed 
in  setting  up  a command  stream  for  GEMACS.  Cautions  regarding  the  occurrence 
of  coirmands  before  the  occurrence  of  other  conmands  do  exist  and  will  be 
pointed  out  as  they  are  encountered. 

The  coimands  will  be  discussed  pretty  much  in  the  order  In  which 
they  would  appear  In  a typical  GEMACS  analysis.  The  general  form  of  the 
command  will  be  given  and  those  items  that  may  be  defaulted  will  be 
Indicated  by  the  use  of  square  brackets  ([  ])  in  the  figures.  The  presence 
of  braces  ({  })  in  the  figures  Indicates  that  a multiple  choice  exists 
for  a particular  item.  The  need  for  a user-defined  symbol  is  Indicated 
by  the  letter  S.  If  a previously  defined  symbol  is  required,  this  will 
be  shown  by  the  letters  DS.  If  the  variable  can  use  either  a new  or 
previously  defined  symbol,  this  will  be  Indicated  by  the  letters  SDS. 

These  codes  along  with  others  used  in  the  commands  are  shown  in  Figure  49. 


The  symbolic  names  that  are  defined  by  the  user  must  in  all  cases 


S - PREVIOUSLY  UNDEFINED  SYMBOL 
DS  - PREVIOUSLY  DEFINED  SYMBOL 
SDS  - EITHER  S OR  DS 


N - NUMERIC  VALUE 
DSN  - EITHER  DS  OR  N 
KW  - KEYWORD 

A - ALPHANUMERIC  WORD  UP  TO  SIX  CHARACTERS 
(A-Z  AND  0-9)  WITH  THE  FIRST  CHARACTER 
BEING  (A-Z) 


FIGURE  49.  COMMAND  CODES 
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consist  of  one  to  six  characters,  the  first  character  being  one  of  the 
letters  of  the  alphabet  (A-Z).  The  other  five  characters  may  be  letters 
of  the  alphabet  or  one  of  the  digits  (0-9).  An  additional  restriction 
is  that  certain  keywords  are  used  by  GEMACS  for  specific  purposes  (such  as 
default  names)  and  may  not  be  used  as  symbolic  names.  These  are  shown 
in  Figure  50. 

There  yet  remains  one  final  restriction  on  the  definition  of  symbolic 
names.  Suppose  the  Interaction  matrix  Is  labelled  XXXXXX,  where  the  X's 
are  legal  characters.  Then  If  the  matrix  Is  to  be  decomposed  Into  lower 
and  upper  triangular  matrices,  the  user  must  not  define  symbols  such  as 
XXXLWR,  XXXUPR  OR  XXXPVT  since  these  names  will  be  generated  internally 
by  the  code  to  Identify  the  files  In  which  the  decomposed  matrices  are 
stored.  However,  as  will  be  shown  later,  these  symbols  may  be  referenced 
in  PRINT,  WRITE  and  PURGE  commands. 

After  the  general  form  of  the  command  there  will  be  shown  a specific 
example  of  the  command  geared  to  the  solution  of  the  example  problem. 

A complete  list  of  the  commands  is  given  in  Table  1 of  the  GEMACS  user's 
manual  [10]. 

• | 

A.  Introductory  Comments 

In  setting  up  a deck  of  Input  commands  it  is  best  to  preface  the  deck 
with  a series  of  comment  cards  that  somewhat  identify  the  problem  being 
worked.  Since  the  input  deck  is  the  first  part  of  the  hard  copy  printout, 
this  serves  to  readily  Identify  this  particular  analysis  in  a stack  of 
similar  outputs.  A sample  that  Is  being  used  for  the  example  is  shown 
In  Figure  51.  Note  that  for  ease  of  Identification  the  dollar  sign  has 


cw 

Cl 

C2 

DM 

DP 

DR 

DT 

DW 

DX 

DY 

DZ 

IS 

LU 

NP 

NR 

ON 

PI 

P2 

R1 

R2 

SC 

SW 

T1 

T2 

VS 

XI 

X2 

Y1 

Y2 

Z1 

Z2 

THREE 

LETTER 

KEYWORDS 

ABS 

CDP 

ECC 

END 

FRQ 

ILP 

INV 

LUD 

OFF 

PHI 

RDP 

SEQ 

SET 

FOUR  LETTER 

KEYWORDS 

AXIS 

BAND 

BNDW 

COND 

EPSR 

ESRC 

LOOP 

PLOT 

PRLC 

READ 

SCDP 

SEGS 

SIZE 

SRDP 

SRLC 

TAGS 

TIME 

TYPE 

VSRC 

■ 1 " 

ZGEN 

ZIMP 

FIVE  LETTER  KEYWORDS 

CONJG 

CPINC 

CPNUM 

DEBUG 

LABEL 

PARTN 

PIVOT 

PRINT 

PULSE 

PURGE 

SOLVE 

THETA 

TRACE  VALUE  WRITE 
SIX  LETTER  KEYWORDS 

BACSUB  CHKPNT  COLPSE  CONVRG  EFIELO  EXPAND  FI LEI D GMDATA 
LINLIN  LINLOG  LINPLR  LOGLIN  LOGLOG  LOGPLR  MAXITR  PCESIN 
REDUCE  REFLCT  REPLCE  RSTART  SINCOS  SVMDEF  TRANSP  WIPOUT 
ZCODES  ZLOADS  ZMATRX 


FIGURE  50.  RESERVED  SYMBOLS 


$••••••••••••••••••••••••••••••••••••••••••••• 

$ 

$ COUPLING  BETWEEN  TWO  MONOPOLES 

$ (0. 25  and  0 .75  WAVELENGTH)  LOCATED  ON 

$ TOP  OF  A COMMUNICATIONS  HUT.  THE  SHORTER 
$ ANTENNA  HAS  A SOURCE  OF  ONE  VOLT.  THE 

$ OTHER  IS  LOADED  BY  50  OHMS.  FREQUENCY 

$ EQUALS  750  MHZ. 

$ 


FIGURE  51.  SAMPLE  INTRODUCTORY  COMMENTS 
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been  indented  to  column  3,  and  the  text  is  set  off  by  the  use  of  lines  of 
asterisks. 

U.  Global  Parameters 

The  first  executable  commands  usually  consist  of  up  to  five  cards 
which  fix  some  of  the  permanent  parameters  of  the  analysis  and  avoid 
the  repetition  of  these  values  on  subsequent  commands.  This 
specifies  the  number  of  files  available  for  use  by  GEMACS,  the  amount 
of  CPU  time  in  minutes  allocated  for  the  job,  the  frequency  of  Interest  In 
MHz  at  which  the  analysis  is  to  be  performed,  and  the  conductivity  and  per- 
mittivity of  the  ground  plane,  if  the  structure  Is  not  perfectly  conducting 
or  if  the  antennas  are  not  located  in  free  space. 

For  this  particular  problem  the  three  cards  shown  in  Figure  52  would 
be  used.  If  the  top  of  the  hut  were  a finite  conductor,  then  the 
conductivity  and  permittivity  would  be  specified  using  the  symbols  COND 
and  EPSR,  respectively.  The  conductivity  is  in  mhos/meter  and  EPSR  Is 
the  dielectric  constant  relative  to  that  of  free  space. 

Note  that  when  NUMFIL  is  set  to  17  this  number  of  files  must  be 
made  available  since  GEMACS  has  no  way  of  checking  on  this.  Two  of 
these  files  are  for  input  and  output  and  their  numbers  are  installation 
dependent.  They  are  presently  set  as  five  for  input  and  six  for  output. 

File  seven  is  the  checkpoint  file.  Files  one  through  four  are  used  for 
internal  scratch  space.  Files  eight  through  17  (in  this  case)  would  be 
used  to  store  data  sets  as  identified  by  symbols  defined  by  the  user. 

Examples  of  such  sets  would  be  the  currents  on  the  subsections,  the  Inter- 
action matrix  and  the  excitation  matrix. 

The  TIME  specification  is  useful  when  one  has  no  Idea  how  long  the 
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NUMFIL  = 17  $ NUMBER  OF  FILES  AVAILABLE 
TIME  = 570  $ NUMBER  OF  CPU  MINUTES  ALLOCATED 
FRQ  = 750.0  $ FREQUENCY  OF  ANALYSIS  IN  MHz 


Figure  52.  SAMPLE  GLOBAL  PARAMETERS 


analysis  will  take.  What  happens  Is  that  when  the  CPU  time  used  by  the 
job  equals  or  exceeds  the  amount  of  time  specified  on  the  TIME  card, 

GEMACS  will  automatically  checkpoint  and  terminate  the  analysis  thus  saving 
the  results  of  the  work  performed  up  to  that  point.  Nothing  will  be  lost 
or  duplicated  when  a restart  Is  performed. 

However,  the  time  specified  on  the  TIME  card  must  be  less  than  that 
specified  on  the  computer  system  job  control  card.  Otherwise,  when 
this  latter  time  Is  reached,  the  system  will  terminate  the  analysis;  and 
the  results  of  the  work  done  since  the  last  checkpoint  will  be  lost. 

It  will  be  possible  to  restart,  but  some  of  the  work  performed  during 
the  previous  execution  will  be  duplicated. 

The  last  of  the  global  parameters  specified  In  Figure  52  Is  the 
frequency  in  megahertz  at  which  the  analysis  Is  to  be  performed.  The 
symbol  for  this  parameter  Is  FRQ.  Its  specification  here  means  that 
this  Item  can  be  defaulted  on  all  commands  requiring  this  parameter. 

C.  Geometry  Generation 

The  next  executable  command  will  read  In  the  geometry  data  set.  This 
command  Is  shown  In  Figure  53.  The  symbol  HUTTOP  has  been  chosen  to 
label  the  geometry  data  set.  In  addition,  the  system  card  reader  Is  the 
file  from  which  the  data  may  be  obtained.  This  Is  the  default  value  for  LU, 
the  logical  unit.  Keep  In  mind  that  If  the  file  name  had  been  defaulted, 
then  the  Item  for  the  geometry  data  set  may  not  appear  In  any  subsequent 
command;  that  Is,  It  must  always  be  a default  parameter  if  It  Is  defaulted 
In  the  GMDATA  command.  The  actual  commands  that  are  used  to  generate  the 
geometry  are  explained  and  Illustrated  In  section  IV  of  this  manual. 


5-9 


si 


: 


GMDATA  [-  S]  [LU  = N] 
GMDATA  = HUTTOP 

Figure  GEOMETRY  GENERATION  COMMAND 


D.  Electrical  Environment 

The  electrical  environment  of  the  structure  to  be  analyzed  should 
be  specified  next.  This  would  Include  the  effects  of  loads,  the  excitation 
(electric  fields  or  voltage  sources  or  a combination  of  both)  and  the 
ground  parameters.  These  last  are  usually  specified  with  the  set  of 
global  parameters,  as  COND  and  EPSR  (Section  V.B).  The  frequency  at 
which  the  analysis  Is  to  be  performed  Is  also  usually  specified  along 
with  the  global  parameters,  although  It  can  also  be  specified  on  a number 
of  other  comnands.  Once  It  Is  specified.  It  may  then  be  defaulted  on  all 
subsequent  comnands.  The  surest  way  to  avoid  error  Is  to  specify  the 
frequency  as  a global  parameter  and  then  default  It  In  all  subsequent 
commands  requiring  the  frequency  as  an  Input.  As  shown  In  Section  IX, 
the  frequency  may  also  be  updated  as  a global  parameter  by  use  of  an 
arithmetic  statement. 

Figure  54  shows  the  general  form  of  the  command  used  to  specify  the 
loading,  as  well  as  the  specific  form  used  In  the  analysis  of  the  sample 
problem.  This  command  allows  one  to  place  loads  on  the  structure  Identified 
In  the  GMDATA  Item  (HUTTOP  In  this  case).  See  Section  V.C.  The  load 
data  will  be  stored  In  the  symbol  SDS,  which  Is  LOAD  In  this  case. 

The  type  of  load  and  the  value  of  the  parameters  are  specified  by  choosing 
one  of  the  possibilities  In  the  large  braces.  COND  Is  used  to  specify  the 
conductivity  of  the  subsections  In  units  of  mhos  per  meter.  ZIMP  Is  used 
to  specify  a lumped  load  Impedance  In  ohms.  Two  values  are  required  for 
this  Item,  the  resistance  and  the  reactance.  In  the  example  problem  there 
Is  a purely  resistive  load  of  50  ohms.  Hence  the  Input  Is  50.0,  0.0. 

Note  that  the  yalues  given  for  COND  or  ZIMP  are  Independent  of  frequency. 
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Individual  circuits,  which  are  functions  of  frequency,  may  also  be 
specified  with  equivalent  overall  values  of  resistance  (R),  Inductance  (L) 
and  capacitance  (C),  in  that  order  of  specification.  Either  parallel  or 
series  circuits  may  be  used,  as  shown  In  Figure  55.  The  units  are  ohms 
for  R,  millihenries  for  L and  microfarads  for  C.  A parallel  or  series 
circuit  is  Input  as  PRLC  or  SRLC,  respectively,  as  shown  In  the  Inner 
braces  in  Figure  54.  Note  that  all  three  values  must  be  Input  even 
if  one  or  more  of  the  quantities  does  not  exist  In  the  equivalent  circuit. 
The  Inductance  and  capacitance  are  functions  of  frequency  in  this  case. 

The  TAGS/SEGS  Item  Identifies  those  subsections  which  are  to  be  loaded 
with  the  values  specified  in  this  command.  In  this  case  segment  156  is  to 
be  loaded  with  50  ohms.  As  can  be  seen  from  Figure  40,  this  is  the  bottom 
subsection  of  the  receiving  antenna. 

If  other  subsections  were  to  have  this  same  load,  then  their  segment 
identifiers  would  also  be  listed,  separated  by  coronas  or  blanks.  If 
consecutive  segments  are  to  be  loaded  with  the  same  value,  then  only  the 
first  and  last  segment  numbers  need  be  listed,  separated  by  a hyphen.  Thus, 
the  TAGS/SEGS  Item  could  be  Input  as  SEGS  equals  156-161,  since  there  are  6 
subsections  on  this  antenna,  If  all  the  subsections  were  to  be  loaded. 

The  alternative  would  be  to  list  all  the  egment  ni  ibers  separated  by 
coronas  or  blanks  (e.g.,  156,  157,  ....  160,  161). 

As  mentioned  In  the  beginning  of  Section  IV,  the  tag  Is  simply  another 
way  of  identifying  a group  of  subsections.  Unlike  the  subsection  number 
It  need  not  be  a number  unique  to  each  segment.  As  a matter  of  fact,  all 
714  subsections  In  HUTTOP  may  have  the  sane  tag,  and  that  tag  may  be  zero 
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FIGURE  55.  EQUIVALENT  CIRCUITS  FOR 

LOAD  SPECIFICATION  COMMAND 
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or  any  other  number.  The  tag  may  thus  be  used  as  a convenient  way  to 
reference  a group  of  subsections,  all  of  which  have  the  same  load,  or 
all  of  which  have  the  same  excitation,  or  all  of  which  belong  to  the  same 
antenna,  or  any  other  division  devised  by  the  user.  It  can  be  thought  of  as 
a flag  identifying  a group  of  subsections  all  of  which  have  a common 
characteristic.  Therefore,  when  the  geometry  was  modeled,  segment  156 
was  given  a tag  of  four;  and  the  identifier  In  the  command  could  have 
been  TAGS  = 4 instead  of  SEGS  ■ 156,  but  only  If  segments  157  to  161  were 
given  some  other  tag  since  they  are  not  loaded  In  this  problem. 

Note  that  more  than  one  ZLOADS  command  may  appear  In  the  command  stream, 
but  each  cortmand  must  reference  the  same  load  data  set  (LOAD  In  this  case). 
Also  more  than  one  type  of  load  may  be  placed  on  the  same  subsection. 

Finally,  and  very  importantly,  all  the  ZLOADS  commands  must  precede  the 
coimand  which  causes  the  Interaction  matrix  to  be  generated.  This  is 
necessary  since  the  symbol  for  the  load  data  set  Is  a required  item  for 
Input  on  the  Interaction  matrix  generation  card  (Section  V.E). 

Figure  56  shows  the  general  form  of  the  command  used  to  specify  the 
antenna  excitation,  as  well  as  the  specific  form  used  in  the  analysis  of  the 
example  problem.  As  can  be  seen,  It  is  quite  similar  to  the  load  speci- 
fication cormand  shown  in  Figure  54.  The  file  for  the  excitation  data  Is 
labeled  as  SOURCE.  The  keyword  VSRC  (voltage  source)  is  a flag  to  the  code 
that  an  excitation  is  being  placed  directly  on  one  or  more  subsections,  as 
opposed  to  illuminating  the  entire  structure  with  an  Incident  electric  field. 
This  latter  can  also  be  done  with  GEMACS , and  Its  command  Is  explained  In 
Section  VIII. 
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SDS  = VSRC  [(DS)]  [FRO  - DSN]  V = DSN,  DSN  ^ags|  = NyN N 


SOURCE  = VSRC  (HUTTPP)  V = 1. 0,0,0  SEGS  = 170 


FIGURE  56.  EXCITATION  SPECIFICATION  COMMAND 
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The  geometry  data  set  Is  referenced  In  the  parentheses  following  the 
keyword.  Recall  that  this  would  be  defaulted  If  the  default  name  were  used 
In  the  geometry  generation  command  (Section  V.C). 

The  frequency  of  the  analysis  may  also  be  defaulted  If  It  has  been 
specified  In  the  list  of  global  parameters  or  on  a previous  conmand. 

Since  It  Is  a global  parameter  In  the  example.  It  Is  shown  defaulted  In 
the  example  In  Figure  56. 

The  value  of  the  voltage  to  be  placed  across  the  terminals  Is  specified 
as  a complex  number,  giving  the  real  and  Imaginary  parts  of  the  excitation. 

In  almost  all  cases  the  voltage  Is  given  as  1.0,  0.0  since  the  electromag- 
netic fields  problem  Is  a linear  one.  In  other  words.  If  the  excitation  were 
Increased  by  some  factor,  say  two  or  ten,  then  the  resultant  field  at  some 
point  In  space  would  be  Increased  by  the  same  factor,  for  example  two  or 
ten.  The  same  Is  true  for  the  calculation  of  coupling.  If  the  Input  voltage 
Is  Increased  by  a factor  of  five,  the  received  voltage  would  be  Increased  by 
. a factor  of  five;  but  the  coupling  (the  ratio  of  received  to  excitation 

voltages)  would  remain  the  same  for  either  value  of  excitation.  Hence,  there 
Is  no  loss  In  generality  by  always  specifying  an  Input  voltage  of  1.0  + jO.O 
volts. 

The  last  Item  In  this  command  Is  the  specification  of  which  subsections 
are  excited.  This  Is  done  In  exactly  the  same  manner  as  for  the  load 
specification  conmand.  The  excitation  In  this  case  Is  placed  across 
subsection  170  (Figure  40).  Since  this  subsection  has  been  given  a tag  of 
3,  the  Item  could  be  TAGS  ■ 3 Instead  of  SEGS  ■ 170,  If  subsection  171 
had  some  other  tag.  The  same  comments  apply  regarding  multiple  cards.  All 
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excitation  specification  commands  must  refer  to  the  same  data  set  (SOURCE); 
and  all  subsections  with  the  same  excitation  can  be  listed  on  the  same 
command,  the  subsection  or  tag  numbers  being  separated  by  commas.  If 
consecutively  numbered  subsections  have  the  same  excitation,  only  the 
first  and  last  subsection  numbers  need  be  given,  separated  by  a hyphen. 

Refer  to  the  discussion  of  the  load  specification  command  for  examples  of 
these  Input  conventions. 

There  Is  one  point  to  keep  In  mind  when  using  this  conrnnd  or  the 
electric  field  excitation  command  (ESRC,  Section  VIII).  If  the  frequency 
at  which  the  analysis  Is  being  performed  Is  changed  (as,  for  example. 

In  a parametric  study.  Section  IX),  then  the  excitation  matrix  Is  reini- 
tialized to  zero.  Consequently,  a new  VSRC  or  ESRC  command  must  be 
input  even  though  the  excitation  on  each  of  the  subsections  remains  the 
same. 

E.  Interaction  Matrix  Generation 

The  general  and  specific  forms  for  this  command  are  shown  In  Figure  57. 
The  keyword  SINCOS  Is  at  this  point  in  time  an  Inflexible  quantity  Indica- 
ting the  expansion  set  that  Is  used  to  represent  the  current,  much  like  a 
Taylor  series  Is  used  to  approximate  some  mathematical  function  such  as 
the  exponential  or  logarithmic  functions.  Other  expansion  functions  will  be 
Incorporated  Into  future  versions  of  GEMACS  and  different  keywords  will  be 
used,  but  for  the  present  time  only  SINCOS  may  be  entered  for  this  Item. 

Once  again  the  name  of  the  geometry  data  set  must  be  entered  If  the 
default  name  has  not  been  used  In  the  geometry  specification  comnand. 

The  frequency  Is  also  an  Item  In  this  command  that  may  be  defaulted. 
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ZGEN  SINCOS  (GNDATA  = DS3  C™  ■ DSlO  ZMATRX  - S 
[ZLOADS  - DS)  [COND  - DSN]  [EPSR  - DSN] 

ZGEN  SINCOS  GNDATA  = HUTTOP  ZMATRX  = ZIJ  ZLOADS  - LOAD 

FIGURE  57.  INTERACTION  MATRIX  GENERATION  COMMAND 


Since  this  was  Included  as  one  of  the  global  parameters,  It  is  not  Included 
on  the  card  for  the  sample  problem  input  data  deck. 


The  next  item  Is  a required  Item  and  identifies  the  file  In  which  the 
elements  of  the  interaction  matrix  will  be  stored.  Although  up  to  six 
characters  may  be  used  for  a symbolic  name,  ZIJ  Is  being  used  since  the 
entries  in  the  Interaction  matrix  are  referred  to  as  the  "zlj's",  which  Is 
short  for  the  Ijth  element  of  the  Z-matrlx.  Reviewing  the  MOM  matrix 
equation,  (Figure  5),  it  will  be  seen  that  this  is  the  matrix  which  is 
multiplied  by  the  I-matrix  (current  matrix)  to  yield  the  excitation  matrix 
(E-matrix). 

The  next  Item  of  this  command  is  the  Identification  of  the  load 
matrix.  Remember  that  this  symbol  was  defined  on  the  load  specification 
cormand  (Figure  54).  This  item  is  shown  as  one  which  may  be  defaulted,  but 
this  is  true  only  if  no  loads  have  been  placed  on  any  of  the  subsections.  > 

If  a load  matrix  has  been  defined,  its  symbol  must  appear  in  this  command. 

The  next  two  items  of  this  command  specify  the  conductivity  and 
permittivity  of  the  ground  plane,  if  one  exists.  Again,  it  is  recommended 
that  these  quantities  be  specified  in  the  set  of  commands  which  define  the 
global  parameters  and  should  be  defaulted  here.  The  quantities  then  need  to 
be  specified  only  once;  and,  since  they  are  located  in  one  well-defined  area 
of  the  input  deck,  they  can  be  quickly  found  and  modified  as  needed. 

At  this  point  the  global  parameters  and  the  data  for  the  geometry,  for 
the  load  and  the  excitation  have  been  input  into  GEMACS.  From  these  data 
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GEMACS  has  generated  the  elements  of  the  Interaction  matrix  and  the 
excitation  matrix.  It  Is  now  ready  to  solve  the  basic  MOM  matrix  equation 
(Figures  7 and  8). 

F.  Matrix  Equation  Solution 

Two  solution  techniques  presently  exist  within  GEMACS,  and  It  Is 
anticipated  that  more  will  be  added  to  the  code.  These  are  a full-matrix 
solution  technique,  In  which  a Gauss-Jordan  algorithm  Is  used;  and  the 
banded  matrix  Iteration  (BMI)  technique.  In  which  the  Interaction  matrix  Is 
banded,  decomposed  and  used  In  an  Iterative  scheme  to  solve  for  the  currents 
on  the  geometry.  The  former  technique  Is  discussed  In  Section  VII.  The  BMI 
technique  (Figures  7 and  8)  is  discussed  in  this  section. 

The  first  step  is  to  develop  a banded  matrix  from  the  Interaction 
matrix  [Z].  The  object  is  to  get  all  of  the  larger  elements  Into  this 
banded  matrix.  In  addition,  one  must  also  choose  those  elements  that  signi- 
ficantly contribute  to  the  electromagnetic  phenomena  associated  with  the 
problem,  such  as  edges  and  other  nearby  perturbing  Influences.  This  Is 
not  an  easy  task;  and  about  the  only  guide  available  is  past  experience, 
which  at  this  point  Is  next  to  nonexistent.  A bandwidth  (M  In  Figure  7) 
is  selected.  It  should  be  large  enough  to  assure  convergence  to  a 
reasonably  accurate  solution  and  yet  small  enough  to  obtain  that  solution  In 
a practical  amount  of  time. 

It  is  better  to  choose  a large  bandwidth  In  order  to  have  GEMACS 
yield  faster  convergence  than  to  have  a smaller  bandwidth  requiring  a 
large  number  of  Iterations  [5,pp.  3-5], 

In  the  example  all  of  the  Interactions  Involving  those  elements  In 
coordinate  system  number  two  are  placed  within  the  banded  portion  of  the 


5-21 


MBMrtta 


I 


n 

i 


Interaction  matrix.  This  includes  subsections  one  to  312,  which  represent 
the  antennas  and  the  top  of  the  hut.  To  these  are  added  those  subsections 
which  border  the  edge  of  the  hut  in  the  vicinity  of  the  antennas.  Looking 
at  the  figure  in  the  Appendix,  this  Includes  the  first  500  subsections. 

The  bandwidth  (that  is,  the  number  of  diagonals  on  either  side  of  the  main 
diagonal)  is  therefore  499,  which  is  rounded  off  to  500. 

These  are  the  data  that  are  needed  to  construct  the  BAND  command,  shown 
in  Figure  58.  The  symbol  DS  represents  the  name  that  has  been  assigned  to 
the  interaction  matrix.  This  is  "ZIJ"  in  the  example.  The  value  for  BNDW 
Is  the  number  of  diagonals  on  one  side  or  the  other  of  the  main  diagonal 
that  will  appear  In  the  banded  matrix,  whose  symbol  is  given  by  SDS  in  the 
general  form  of  the  command.  Therefore,  in  the  example  the  value  of  N is 
500. 

A restriction  needs  to  be  kept  In  mind  when  naming  the  banded  matrix. 

The  same  symbolic  name  may  not  be  used  for  both  the  banded  matrix  and  the 
original  Interaction  matrix.  The  reason  is  that  both  matrices  are  referenced 
in  the  comnand  that  implements  the  BMI  technique  and  some  means  must  be 
available  to  differentiate  between  them. 

Therefore,  for  the  example  the  symbol  BNDZIJ  has  been  chosen  since  it 
is  almost  an  acronym  for  "banded  ZIJ  matrix."  The  final  form  of  the  command 
that  will  be  used  in  the  solution  of  the  problem  is  also  shown  in  Figure  58. 

Once  the  Interaction  matrix  has  been  banded,  the  banded  portion  must  be 
decomposed  into  lower  and  upper  triangular  matrices  as  shown  in  the 
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SDS  = BAND  (DS)  BNDW  = N 


BNDZIJ  = BAND  (ZIJ)  BNDW  = 500 


Fisure  58.  THE  BAND  COMMAND 
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solution  process  Indicated  in  Figure  8.  This  Is  accomplished  by  using  the 
matrix  decomposition  (LUD)  conmand  shown  In  Figure  59. 

The  symbol  DS  Is  the  name  of  the  banded  matrix  that  has  been  developed. 

In  this  case  BNDZIO.  The  symbol  for  the  decomposed  matrix  may  be  either  a 
new  name  or  some  previously  defined  name.  Also,  In  the  case  the  name  of 
the  decomposed  matrix  may  be  the  same  as  the  name  for  the  banded  matrix. 

This  Is  the  method  usually  chosen,  and  that  Is  what  is  in  Figure  59. 

Once  the  banded  matrix  has  been  decomposed,  the  BMI  comnand  Is  used  to 
solve  the  MOM  matrix  equation  for  the  current  on  the  wires  in  the  grid  model. 
The  form  of  this  command  and  the  specific  command  to  be  used  to  solve  the 
example  are  shewn  In  Figure  60.  In  this  command  DS1  represents  the  name 
of  the  decomposed  banded  matrix,  BNDZIJ  in  this  example.  The  currents  that 
are  to  be  calculated  are  stored  under  the  name  represented  by  SDS1 , which 
In  the  example  Is  CURENT.  The  excitation  that  has  been  set  up  Is  stored  in 
DS2,  or  SOURCE  in  the  example.  The  symbol  DS3  is  the  name  of  the  interaction 
matrix  that  was  calculated.  Recall  that  It  Is  ZIJ.  The  final  symbol  In  the 
equation  is  SDS2,  the  set  of  currents  from  the  previous  iteration.  It  may 
be  the  same  as  SDS1.  An  initial  value  may  be  pre-loaded  in  SDS2  (Section 
V.K),  but  it  has  been  found  that  leaving  the  initial  currents  on  the 
subsections  equal  to  zero  does  not  significantly  decrease  the  solution  time. 
Thus,  the  array  for  the  currents  is  usually  not  preloaded  and  SDS2  is  nor- 
mally the  same  as  SDS1. 

The  next  item,  CONVRG,  is  the  convergence  criterion  that  will  be  used  to 
test  the  currents  that  have  been  derived  at  each  iteration.  At  this  tine 
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SDS  = LUD  CDS) 
BNDZIJ  = LUD  (BNDZIJ) 


Figure  59.  MATRIX  DECOMPOSITION  COMMAND 


GEI'ACS  has  three  convergence  tests  built  Into  It.  The  pros  and  cons  of  each 
of  these  tests  have  been  discussed  in  the  BDM  quarterly  reports  [6,8]  and 


are  touched  on  briefly  in  the  engineering  manual  [9].  Basically  they  are  the 
Boundary  Condition  Relative  Error  (BCRE) , the  Iterative  Relative  Error  (IRE) 
and  the  Predicted  Relative  Error  (PRE) . The  BCRE  compares  the  boundary 
condition  obtained  using  the  calculated  currents  to  the  original  excitation. 
The  IRE  compares  the  currents  at  one  iteration  to  those  obtained  from  the 
previous  iteration,  convergence  being  attained  when  there  Is  a very 
small  change  in  the  value  of  the  new  set  of  currents.  The  PRE  Is  a function 
of  the  IRE  and  Is  based  on  the  fact  that  the  latter  can  be  represented 
In  an  approximate  fashion  by  a exponential  function,  especially  when  the 
convergence  rate  is  fairly  rapid.  The  PRE  was  found  to  yield  more  useful 
data  relating  to  convergence;  and,  as  such,  was  felt  to  be  the  best  of 
the  three  available  convergence  tests.  As  a consequence,  the  default 
,alue  for  this  item  is  the  PRE.  This  is  the  one  used  in  the  solution  of  the 
example  problem. 

The  next  item  is  the  value  of  the  convergence  criterion  in  percent  at 
which  the  iteration  process  will  terminate.  The  values  of  each  of  the 
convergence  criteria  are  computed  at  each  step  of  the  iteration  process 
according  to  equations  73,  74  or  76  of  the  engineering  manual  [9].  For 
the  predicted  relative  error  (PRE)  it  is  basically  a function  of  the  ratio 
of  the  current  vector  at  one  iteration  to  the  current  vector  at  the 
preceding  iteration. 

The  choice  of  VALUE  is  dependent  on  the  problem  and  what  quantity  is 
being  investigated.  If  it  is  the  electric  field  far  from  the  antenna,  a 
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value  of  ten  percent  will  yield  good  results.  When  working  close  to  the 
antennas,  as  In  the  example,  the  VALUE  Item  should  be  given  Its  default 
value,  which  Is  one  percent.  This  requires  more  Iterations,  but  the  nature 
of  the  problem  necessitates  the  acceptance  of  the  Increased  cost. 

The  final  Item  in  this  command  Is  the  maximum  number  of  Iterations  that 
will  be  allowed.  It  may  happen  that  the  way  the  grid  Is  set  up  or  the 
choice  of  the  bandwidth  casues  the  solution  to  converge  very  slowly.  In 
that  case  an  answer  will  be  obtained  eventually,  but  at  a high  cost.  The 
MAXITR  Item  provides  a check  which  will  not  allow  GEMACS  to  grind  away 
forever.  Suppose,  for  example,  that  the  value  of  MAXITR  Is  ten.  If 
convergence  Is  not  reached  after  the  tenth  Iteration,  GEMACS  will  terminate 
the  analysis  and  print  out  the  values  of  current  at  that  point.  It  will  also 
output  the  value  of  the  convergence  criterion  In  percent  and  predict  how 
many  more  Iterations  are  required.  Based  on  this  Information  the  value 
of  MAXITR  can  be  changed  If  not  too  many  more  Iterations  are  required  or 
the  bandwidth  may  be  Increased  In  an  attempt  to  obtain  a faster  rate  of 
convergence;  or,  as  a last  resort,  a new  wire  grid  model  may  be  set  up. 

The  analysis  Is  then  restarted  from  the  last  checkpoint  Section  V.H). 

Note  that  this  Item  cannot  be  defaulted.  A nominal  value  to  use  as 
a starting  point  Is  ten.  The  form  of  this  command  that  is  used  In  the 
analysis  of  the  problem  Is  shown  in  Figure  60. 


At  this  point  GEMACS  has  now  calculated  the  current  on  each  of  the 
subsections.  These  data  are  all  that  are  needed  to  solve  the  coupling 
problem.  Therefore,  the  next  operation  Is  to  have  GEMACS  print  out  these 


current  values.  The  commands  that  are  used  for  control  of  the  output  are 
discussed  In  Section  VI. 

G.  Checkpoint  and  Example  Problem  Input  L i s ting 


It  has  been  mentioned  several  times  that  GEMACS  has  a checkpoint  and 
restart  capability.  This  feature  allows  the  user  to  restart  his  analysis 
at  any  point  In  the  stream  if  It  has  been  Interrupted  for  any  reason 
whatsoever.  The  importance  is  that  the  user  may  restart  from  very  nearly 
the  point  at  which  the  analysis  had  left  off,  thus  eliminating  a lot  of 
costly  duplication. 

In  order  to  make  use  of  the  restart  option  one  or  more  checkpoints 
must  have  been  taken  during  the  analysis  by  using  the  checkpoint  command. 
When  a checkpoint  is  taken,  a snapshot  of  the  computer  status  is  stored 
away  on  some  peripheral  device  or  file,  such  as  a magnetic  tape.  All  of 
the  data  input  to  or  generated  by  GEMACS  are  stored  in  an  orderly  manner 
on  this  file.  In  this  way  the  status  of  the  analysis  at  that  point  in 
time  is  saved;  and  it  is  available  for  later  use  as  a starting  point  from 
which  a modified  analysis  stream  may  be  begun. 

The  checkpoint  command  is  shown  in  Figure  61.  Since  most  of  the  items 
can  be  defaulted,  it  is  generally  a very  simple  input  card  to  generate. 

The  symbol  CHKPNT  signifies  the  type  of  command.  The  first  item  (LU) 
designates  what  file  number  the  checkpoint  tane  or  disc  is  loaded  on. 

This  is  one  of  those  files  that  must  be  accounted  for  when  specifying 
the  number  of  files  as  a global  parameter  (NUMFIL  in  Figure  52).  As  a rule 
one  uses  the  default  value  which  puts  the  checkpoint  storage  device  on 
logical  unit  seven. 
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The  FILEID  Item  can  be  used  to  label  the  checkpoint  file  with  an 
Identifier  significant  to  the  job.  This  is  useful  when  several  checkpoint 
files  are  floating  around.  GEMACS  compares  this  name  with  the  one  that 
Is  given  In  the  restart  command  (Section  V.H);  and  If  the  two  differ,  a 
warning  Is  printed  out  advising  the  user  to  double  check  to  make  certain 
that  the  correct  checkpoint  data  Is  being  used.  The  default  name  for  the 
file  In  this  and  the  restart  command  Is  CHKPNT. 

The  next  Item,  CPINC,  specifies  how  often  the  checkpoint  Is  to  be 
taken,  counting  In  minutes  of  computer  processing  time  (not  wall  clock 
time).  This  Is  the  actual  time  that  the  computer  is  working  on  the  problem. 
For  longer  problems  It  is  reconmended  that  a checkpoint  be  taken  every 
60  minutes.  Taking  a checkpoint  on  a large  problem  may  require  five 
minutes;  and  If  taken  too  often,  a large  percentage  of  the  time  could 
be  spent  just  transferring  data  to  the  checkpoint  file.  Separating  the 
checkpoints  by  too  great  a time  runs  the  risk  of  termination  a long  time 
after  the  last  checkpoint,  which  may  cause  expensive  duplication  of  all 
the  operations  since  that  last  checkpoint  was  taken. 

If  this  Item  Is  defaulted,  an  immediate  checkpoint  Is  taken  regardless 
of  what  task  has  just  been  completed.  For  reasons  to  be  explained  shortly, 
a checkpoint  should  be  taken  just  after  the  Interaction  matrix  has  been 
completely  generated.  In  this  way  a record  of  the  completed  matrix  will 
be  available  If  a change  In  the  bandwidth  is  desired.  Remember,  the  genera- 
tion of  the  Interaction  matrix  Is  an  expensive  process. 

The  final  Item,  NR,  provides  control  over  whether  or  not  an  historical 
record  of  the  checkpoints  taken  Is  kept.  When  the  symbol  "NR"  Is  typed 
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on  the  card,  GEMACS  will  not  rewind  the  checkpoint  file.  The  next  check- 
point will  start  right  after  the  last  record  of  the  previous  checkpoint. 

This  Is  true  whether  the  checkpoints  are  taken  at  specified  Intervals  of 
time  or  taken  Immediately.  When  restarting,  any  checkpoint  taken  during 
the  course  of  the  analysis  will  be  available  for  access. 

If  this  Item  Is  defaulted,  then  after  each  checkpoint  is  taken  the 
file  Is  rewound  and  the  subsequent  checkpoint  will  then  start  at  the  be- 
ginning of  the  file.  In  effect  one  Is  wiping  out  checkpoints  as  new 
data  are  stored  on  the  file.  The  disadvantage  of  this  course  of  action 

* I 

Is  that  when  restarting  one  can  access  only  the  last  checkpoint  written. 

Rewinding  and  overwriting  save  storage  space,  but  they  limit  the  flexibility 
when  restarting. 

Sample  checkpoint  conmands  are  also  shown  In  Figure  61.  The  examples 
are  Identical  to  the  conmands  used  In  the  analysis  of  the  sample  problem. 

In  the  first  example,  it  is  specified  that  a checkpoint  will  be  taken 
after  every  60  minutes  of  computer  processing  time.  Note  that  the  file 
is  to  be  rewound  before  each  new  checkpoint  Is  taken.  The  historical 
record  will  be  lost,  but  storage  space  will  be  saved;  this  will  he 
significant  since  there  are  so  many  elements  In  the  Interaction  matrix. 

Such  a timed  checkpoint  command  should  appear  early  In  the  task  command 
Input  stream.  In  the  example  It  Is  just  after  the  read  geometry  (GMDATA) 
command,  as  shown  In  Figure  62. 

The  second  example  in  Figure  61  shows  what  an  Immediate  checkpoint 
command  would  look  like.  The  default  values  for  the  file  name  and  device 
number  have  again  been  used,  as  v?ell  as  a default  for  the  "MR"  parameter. 
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!♦»♦ **»*♦*****♦**♦**»♦*»*»*»»**♦»* *******♦»*****♦****♦»*»***»»*********** 
I************************************************************************ 
$ 

$ COUPLING  BETWEEN  TWO  MONOPOLES  (0.25  and  0.75)  LOCATED  ON 
$ TOP  OF  A COMMUNICATIONS  HUT.  THE  SHORTER  ANTENNA  HAS  A 

$ SOURCE  OF  ONE  VOLT.  THE  OTHER  IS  LOADED  BY  50.0  OHMS.  THE 

$ FREQUENCY  EQUALS  750.0  MHZ. 

$ 

f ************************************************************* *********** 
$ ************************************************************************ 
NUMFIL-17 
TIME-570 
FRQ- 750.0 
DEBUG  ON  ILP 
GMDATA-HUTTOP 
DEBUG  OFF 
CHKPNT  CPINC-60 

ZLOADS-LOAD  GMDATA-HUTTOP  ZIMP-50. 0,0.0  SEGS-156 
SOURCE-VSRC (HUTTOP ) V-l. 0,0.0  SEGS-170 
ZGEN  SINCOS  GMDATA-HUTTOP  ZMATRX-ZIJ  ZLOADS-LOAD 
CHKPNT 

BNDZIJ-BAND(ZIJ)  BNDW-500 
BNDZI J-LUD(BNDZIJ) 

CHKPNT 

BNDZI J*CURENT-SOURCE-ZIJ*CURENT  MXITR-10 
PRINT  SOURCE  CURENT 
END  OF  COMMAND  STREAM 

FIGURE  62.  EXECUTION  LANGUAGE  INPUT  STREAM 
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Finally,  no  input  Is  specified  for  the  increment  parameter.  Therefore, 
the  checkpoint  will  be  taken  as  soon  as  this  command  is  encountered  in 
the  execution  sequence.  As  can  be  seen  from  the  sample  problem  input  in 
Figure  62,  checkpoints  are  to  be  taken  as  soon  as  the  interaction  matrix 
has  been  generated  and  just  after  the  banded  matrix  has  been  decomposed. 

Whenever  a checkpoint  Is  taken  during  the  analysis,  a standard  message 
is  printed  on  the  GEMACS  output.  An  example  of  this  is  to  be  seen  on 
page  20  of  the  Appendix.  The  first  line  gives  the  unique  numerical  iden- 
tifier of  the  checkpoint  and  the  CPU  time  at  which  the  checkpoint  Is 
being  started.  There  then  follows  a listing  of  the  names  of  the  common 
areas  being  stored  on  the  file  as  well  as  the  number,  name  and  size  of  all 
peripheral  files  that  have  been  generated  during  the  analysis.  The  last 
part  of  the  printout  gives  the  CPU  time  at  which  the  checkpoint  was  com- 
pleted, the  length  of  time  it  took  to  take  the  checkpoint  and  the  size 
of  the  checkpoint  file  in  number  of  words  written.  All  times  are 
given  in  minutes  of  computer  processor  time  (CPU). 

With  the  exception  of  the  output  commands  (Section  VI)  a discussion 
of  all  the  execution  coirmands  that  are  Included  in  the  original  Input 
to  perform  the  coupling  analysis  using  the  BMI  technique  has  been  completed. 
The  complete  execution  connand  input  stream  is  shown  In  Figure  62.  This 
is  pretty  nuch  the  way  that  the  commands  would  be  sequenced  In  order  to 
obtain  the  currents  that  are  flowing  on  any  structure.  Still  other  com- 
mands are  available,  and  they  are  discussed  In  the  following  four  sub- 
sections and  sections  VI  through  IX. 
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H.  Restarting 


When  analyzing  a structure  which  requires  a large  number  of  subsections. 
It  Is  necessary  to  divide  up  the  process  Into  several  sittings  or  "runs". 
This  means  that  the  job  must  be  submitted  more  than  once  possibly  be- 
cause other  jobs  require  the  computer  resources  being  used.  The  restart 
command  Initiates  a run  someplace  In  the  middle  of  the  analysis  Instead 
of  starting  over  from  the  beginning  each  time.  Th^s  is  made  possible  by 
using  the  data  that  have  been  stored  away  on  the  checkpoint  tape. 

The  format  for  the  restart  command  Is  shown  In  Figure  63.  The  symbol 
Identifying  the  type  of  the  conmand  Is  RSTART.  The  next  two  items, 

LU  and  FILE  ID,  are  the  same  as  those  In  the  checkpoint  command  (Section 
V.G).  If  not,  then  the  entries  for  these  two  Items  must  be  identical  In 
both  types  of  commands. 

The  last  Item  In  this  conmand,  CPNUM,  specifies  the  checkpoint  number 
from  which  the  restart  Is  to  begin.  As  explained  In  Section  V.G.,  each 
checkpoint  Is  given  an  unique  numerical  Identifier.  This  Is  true  whether 
the  file  is  rewound  between  checkpoints  or  not  and  whether  It  Is  a timed 
checkpoint  or  not.  It  Is  this  Identifier  that  Is  Input  for  the  CPNUM 
item.  As  can  be  seen  In  the  sample  problem  output  on  page  22  of  the 
Appendix,  the  checkpoint  number  immediately  precedes  the  description  of 
the  data  that  are  being  stored  on  the  file.  The  rest  of  the  description 
includes  a list  of  the  common  blocks  read  from  the  checkpoint  file  together 
with  their  sizes,  followed  by  a list  of  the  names  and  sizes  of  all  the 
peripheral  files.  If  a name  were  given  to  the  checkpoint  file.  It  would 
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appear  before  the  number  of  this  file  and  the  number  of  the  checkpoint. 

The  example  shown  In  Figure  63  shows  what  the  total  execute  command 
list  would  look  like  If  one  were  simply  restarting  from  some  point  In  the 
analysis  after  an  execution  was  prematurely  terminated.  It  would  contain 
only  the  restart  and  end  (Section  V.J)  commands. 

It  should  be  pointed  out  that  one  logical  unit  can  be  used  to  restart 
from  and  a second  unit  used  to  take  subsequent  checkpoints  on.  In  this 
way  the  Integrity  of  the  data  on  the  original  checkpoint  file  can  be 
maintained.  To  do  this  some  logical  unit  (LU)  other  than  seven  must  be 
specified  in  the  restart  command  and  a clean  storage  file  loaded  on 
seven  to  accumulate  later  checkpoint  data.  A good  alternate  to  seven  would 
be  three  or  four,  since  at  present  these  are  unused. 

In  addition  to  simply  restarting  one  can  modify  the  original  Input 
stream;  for  example  to  Increase  or  decrease  the  bandwidth  with  which  the 
BAUD  and  BMI  conmands  are  to  work.  This  Is  done  by  using  the  WIPOUT 
comnand,  whose  format  Is  shown  In  Figure  64.  The  symbol  WIPOUT  Indicates 
the  type  of  command,  and  the  N's  are  the  command  sequence  numbers  to  be 
wiped  out.  Note  that  there  Is  a difference  between  card  numbers  and  com- 
mand sequence  numbers.  The  comment  cards  shown  In  Figure  62  are  cards  one 
through  ten,  but  they  are  not  commands.  Card  11  contains  command  one,  and 
card  25  Is  command  15.  Also  If  a command  requires  one  or  more  continuation 
cards,  the  continuation  cards  do  not  Increase  the  command  sequence  number. 

The  first  example  In  Figure  64  Illustrates  the  process  of  eliminating 
particular  commands  from  the  original  task  list.  Any  number  of  Individual 
tasks  may  be  deleted  In  this  manner. 
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WIPOUT  N N ...  N 


— OR — 


W I POUT  N N 99999 


WIPOUT  3 6 12 


(fictitious) 


WIPOUT  12  99999 


Figure  64.  THE  WIPOUT  COMMAND 
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There  Is  a shortcut  for  the  WIPOUT  command,  and  Is  also  shown  in 
Figure  64.  The  number  99999  has  the  significance  that  all  commands, 
starting  with  the  command  whose  sequence  number  Is  equal  to  the  value  of 
the  N preceding  the  number  99999  down  to,  and  Including,  the  WIPOUT 
command  are  all  deleted.  To  see  this,  think  of  the  list  of  commands  as 
generating  a task  list.  Such  a list  can  be  seen  on  pages  A-10  and  A-ll  of 
the  Appendix.  An  explanation  of  this  list  Is  to  be  found  In  Section  VI. E. 
The  commands  In  the  restart  also  generate  a list  which  Is  appended  to  the 
bottom  of  the  original  list.  If  the  N were  12  then  all  the  tasks  from  12 
through  the  wipe  opt  command  would  be  eliminated.  The  second  example  In 
Figure  64  shows  such  a command. 

The  use  of  this  second  example  In  the  Input  stream  for  the  sample  pro- 
blem (Figure  62)  would  eliminate  all  the  commands  below  the  first  Immediate 
checkpoint;  that  Is,  from  the  BAND  command  down.  The  WIPOUT  command  could 
then  be  followed  by  a new  BAND  command,  specifying  a different  bandwidth. 
This  new  command  would  then  need  to  be  followed  by  a sequence  of  commands 
Identical  to  the  original  sequence  since  that  original  sequence  was  deleted. 
One  command  cannot  be  replaced  by  a new  command  and  not  disturb  the  se- 
quence following  the  command  that  was  replaced.  The  original  command 
would  be  deleted  and  the  replacing  command  put  out  of  place  at  the  end  of 
the  modified  task  list.  This  Is  further  explained  In  Section  VI. E. 

In  the  second  part  of  the  sample  problem  shown  In  the  Appendix,  the 
analysis  Is  restarted  after  the  Interaction  matrix  has  been  generated.  This 
Is  shown  on  page  A-48  of  the  Appendix.  A full  matrix  solution  process 
(Section  VII)  Is  then  Implemented  to  solve  for  the  subsection  currents. 
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There  Is  one  more  point  to  keep  In  mind.  If  the  operation  has  already 
begun  and  not  been  completed.  Its  command  cannot  be  wiped  out.  As  can 
be  seen  from  the  sample  output,  the  first  checkpoint  occurred  while  the 
Interaction  matrix  generation  command  (ZGEN)  was  being  executed.  It  could 
not  have  been  deleted  In  the  restart.  This  Is  one  advantage  of  the  Immediate 
checkpoint  as  used  In  the  sample  Input  stream.  The  full  matrix  has  been 
generated  and  stored  , and  the  actual  solution  process  has  not  started.  Thus, 
a solution  process  different  from  the  original  can  be  used  on  a restart. 


As  a matter  of  fact,  this  Is  what  was  done  starting  on  page  49  of  the 
Appendix.  Also  the  bandwidth  can  be  easily  changed  on  a restart  since  the 
original  BAND  command  has  not  yet  begun  to  be  executed. 

I.  PURGE 

In  order  to  save  as  much  storage  space  as  possible,  all  files  and  data 
that  will  no  longer  be  used  In  the  analysis  should  be  eliminated.  This  Is 
accomplished  by  the  use  of  the  PURGE  command,  shown  In  Figure  65.  This 
also  has  the  advantage  of  making  the  checkpoint  file  shorter  and  reducing 
the  time  It  takes  to  transfer  the  data  to  this  file. 

The  command  itself  Is  quite  simple.  It  consists  of  the  mnemonic  PURGE, 
followed  by  the  list  of  symbols  whose  files  are  to  be  purged.  One  example 
that  Is  often  used  Is  to  purge  the  Interaction  natrlx  after  it  has  been 
decomposed.  This  Is  done  only  when  using  the  full  matrix  solution. 

It  cannot  be  done  when  using  the  Bill  solution  process  until  after  the 
currents  have  been  obtained,  since  the  original  Interaction  matrix  is 
referenced  by  the  Bf!I  command. 

t 
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PURGE  SDS  SDS  ...  SDS 
PURGE  BNDZIJ  BNDUPR  BNDLWR 


Figure  65.  THE  PURGE  COMMAND 


The  example  shown  In  Figure  65  is  taken  from  the  command  language 
stream  used  in  the  second  part  of  the  sample  problem.  In  this  part  a full 
matrix  solution  is  used  to  find  the  currents  on  the  wire  grid  model. 
Therefore,  the  banded  matrix  and  its  upper  and  lower  derivatives  are  no 
longer  needed.  The  command  shown  In  Figure  65  will  eliminate  these  files 
from  storage  and  from  all  future  checkpoints. 

J.  END 

The  final  command  that  will  appear  in  the  task  execution  conmand  stream 
is  the  END  command.  As  shown  In  Figure  66,  it  consists  simply  of  the 
three-letter  symbol  "END".  It  Is  used  to  indicate  to  GEMACS  that  all  the 
commands  ha  been  input  and  processed. 

As  is  the  case  with  the  END  command  in  the  geometry  input  (Section  IV. A), 
there  is  no  need  to  separate  the  comment  by  a dollar  sign  ($).  All  that 
is  needed  is  a blank  space  or  a comma.  The  message  shown  in  the  example 
in  Figure  66  is  a typical  one  used  in  all  input  decks. 

K.  The  SET  Command 

This  command,  shown  in  Figure  67,  is  used  to  initialize  or  change  the 
data  associated  with  some  set.  One  use  of  this  command  was  mentioned  in 
connection  with  the  banded  matrix  Iteration  command  (Section  V.F).  There 
It  was  stated  that  an  initial  value  could  be  used  for  the  currents  to  possi- 
bly speed  up  the  Iteration  process.  Use  of  this  command  could  also  be  an 
alternate  way  to  develop  the  excitation  matrix.  Instead  of  using  the  vol- 
tage source  connand  (VSRC)  or  to  develop  the  load  matrix  (ZLOADS  command), 
both  discussed  in  section  V.D. 
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SET  SDS  = N,  im  [Rl=Nl  [R2=N ] CC1=N ] tC2=NJ 


SET  SOURCE  = 0.0,0, 0 R2  = 100 
SET  SOURCE  = 1. 0,0.0  R1  = 33 


(fictitious) 


Figure  67.  THE  SET  COMMAND 


As  can  be  seen  In  Figure  67,  the  mnemonic  for  this  command  is  SET.  The 
symbol  SDS  represents  the  name  of  the  file  being  generated.  The  data  are 
input  in  either  real  or  complex  form,  indicated  by  “N,[N]“.  If  the  data 
are  real,  the  second  value  is  defaulted.  If  the  data  are  complex,  then  tv/o 
values  must  be  input,  even  if  one  or  the  other  is  zero.  Note  that  only  one 
real  or  one  complex  value  may  be  input  per  command.  This  value  will  then 
be  used  for  all  the  elements  that  exist  within  the  row  and  column  limits 
specified  in  the  command. 

The  last  four  items  shown  in  the  command  in  Figure  67  specify  these  row  - 
and  column  limits.  Any  or  all  of  these  items  may  be  defaulted.  R1  and  R2 
are  the  row  limits  of  the  elements  to  which  the  values  will  be  applied.  Cl 
and  C2  are  the  lowest  and  highest  column  limits,  respectively.  The  default 
value  is  one  for  R1  and  Cl;  and  R2  and  C2  are  defaulted  to  R1  and  Cl,  respec- 
tively. This  latter  holds  true  even  if  R1  and  Cl  have  been  defaulted. 

An  example  of  the  use  of  this  command  is  also  shown  in  Figure  67.  It 
assumes  that  an  excitation  matrix,  labeled  SOURCE,  already  exists  and  that 
it  consists  of  100  elements.  It  is  desired  to  change  the  boundary  condi- 
tion such  that  only  one  subsection  is  excited  by  one  volt.  The  first 
command  will  set  all  of  the  elements  to  zero.  Since  the  excitation  is  a 
complex  number,  two  values  are  input,  one  real  and  one  imaginary.  The  default 
value  for  R1  has  been  used,  while  R2  has  been  set  to  the  number  of  rows  in 
the  matrix.  Since  there  is  only  one  column  in  the  matrix.  Cl  has  been 
defaulted  to  one;  and  C2  has  been  defaulted  to  Cl. 

The  second  cormand  then  sets  element  33  to  one  volt.  In  this  command 
R2  has  been  defaulted  to  Rl.  Again,  Cl  and  C2  have  each  assumed  their 
default  values. 
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Note  that  If  element  66  was  also  to  have  been  set  to  one  volt,  then  a 
separate  command  would  have  been  needed.  Setting  R2  to  66  In  the  second 
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VI.  OUTPUT 


There  are  three  types  of  output  provided  by  GEMACS:  the  standard 
boiler-plate,  those  data  specifically  requested  by  the  user;  and  error 
messages  and  debug  Information  needed  by  the  user  when  he  has  a problem 
running  the  computer  code.  Examples  of  the  first  two  may  be  seen  In  the 
output  listing  of  the  sample  problem  contained  In  the  Appendix. 
i . A.  Standard  Output 

Figure  68  shows  a list  of  those  types  of  data  that  are  always  printed 
out  by  GEMACS  regardless  of  the  problem  or  the  specific  data  requested  by 
the  user. 

The  first  output  provided  by  GEMACS  Is  the  list  of  keywords  that  may 
not  be  used  as  symbol  names.  This  Is  Identical  to  the  list  that  Is  shown 
In  Figure  50.  This  list  Is  printed  out  strictly  as  a convenience 
In  case  there  is  an  error  message  stating  that  one  of  the  names 
that  Is  used  Is  not  acceptable.  However,  other  reasons  for  such  a message 
would  be  a name  that  Is  greater  than  seven  characters  In  length  or  one 
whose  Initial  character  Is  not  one  of  the  letters  of  the  alphabet. 

The  date  and  time  are  next  provided  by  GEMACS  to  the  user.  This  Infor- 
mation really  comes  In  handy  when  one  starts  piling  up  stacks  of  computer 
paper  and  loses  track  of  the  order  In  which  things  were  done. 

A list  of  the  task  commands  Input  by  the  user  Is  next  printed  out. 
Normally  It  Is  just  a straight  listing  of  this  part  of  the  Input  deck.  How- 
ever, If  any  errors  were  detected  on  any  of  the  commands,  the  error  will  be 
printed  out  Immediately  after  the  erroneous  card.  The  subsequent  commands 
will  also  be  processed  and  errors  searched  for,  but  the  analysis  will  be 
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RESERVED  KEYWORDS 
DATE  AND  TIME 
INPUT  TASK  COMMANDS 
GLOBAL  PARAMETERS  (if  any) 
INPUT  GEOMETRY  COMMANDS 
SEGMENT  DATA  TABLE 
LOADS  (if  any) 

EXCITATION 
TASK  PROCESSING 
ANTENNA/LOAD  PARAMETERS 


Figure  68.  STANDARD  DATA  OUTPUT 


terminated  after  the  last  card  of  the  geometry  Input  has  been  scanned. 

The  global  parameters  will  be  printed  out  next.  If  any  have  been  spec- 
ified. These  data  Include  the  number  of  files  available,  the  time  allowed 
for  execution,  the  frequency  In  megahertz  at  which  the  analysis  Is  being 
performed,  and  the  conductivity  and  permittivity  of  any  Imperfect  ground. 

The  global  parameters  are  followed  by  a listing  of  the  geometry  Input 
deck.  As  In  the  case  of  the  command  language  deck  this  will  be  a straight 
listing  unless  an  error  was  detected  In  the  format.  Typographical  errors 
that  result  In  the  erroneous  location  of  one  or  more  subsections  will  not 
cause  an  error  message.  These  will  have  to  be  ferreted  out  by  the  user. 

These  first  five  Items  are  printed  out  mainly  as  a check  for  the  user 
to  satisfy  himself  that  what  he  thought  he  put  Into  GEMACS  is  actually  what 
did  get  Into  GEMACS.  The  first  analyzed  data  that  are  printed  out  are  the 
specifics  relating  to  each  of  the  subsections.  The  data  are  contained 
In  the  segment  data  table,  whose  headings  are  listed  In  Figure  69.  Each 
line  In  the  table  contains  the  data  for  one  particular  subsection. 

The  first  column  represents  the  tag  number  assigned  to  each  subsection 
as  It  Is  generated.  Looking  at  the  output,  It  Is  seen  that  most  of  the 
first  312  subsections  have  a tag  of  two.  This  Indicates  subsections  of 
the  fine-mesh  grid  which  was  defined  In  coordinate  system  number  two. 

Two  subsections  each  have  a tag  of  three.  Indicating  the  transmitting 
antenna.  Six  subsections  have  a tag  of  four  Indicating  the  receiving 
antenna.  The  rest  of  the  subsections  belong  to  the  coarse-mesh  grid  and 
have  a tag  of  one. 
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TAG 

XN 

XC 

XP 

YN 

YC 

YP 

ZN 

ZC 

ZP 

RADIUS 

LENGTH 

END1 

ISEG 

END2 


Figure  69.  SEGMENT  DATA  TABLE  HEADINGS 
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The  next  nine  columns  list  the  (x,y,z)  coordinates  of  the  ends  and  cen- 
ter points  of  the  subsection.  The  letter  "N"  stands  for  the  first  coordi- 
nate point  on  the  wire  Input  (WR)  card;  the  letter  "P"  stands  for  the  second 
coordinate  point  on  the  wire  card.  The  letter  "CH  Indicates  the  center 
coordinate  of  the  subsection.  Therefore,  reading  the  table  from  left  to 
right,  the  start,  center  and  end  points,  respectively,  for  each  coordinate 
are  given.  By  reading  this  table  carefully  one  can  then  tell  whether  or  not 
each  subsection  Is  properly  positioned  In  the  wire  model.  Note  that  these 
values  are  In  meters  In  the  global  coordinate  system  only. 

The  next  two  columns  of  the  segment  table  list  the  radius  and  length 
In  meters,  respectively,  of  each  subsection  of  the  wire  model.  By  glancing 
through  the  table  It  Is  seen  that  each  subsection  has  a radius  of  0.000776 
meters  and  that  the  length  of  each  of  the  first  312  subsections  Is  0.05 
meters.  For  the  coarse-mesh  part  of  the  grid  the  length  of  each  subsection 
Is  either  0.05  or  0.1  meters.  These  data  are  to  be  expected  based  on 
the  figures  In  Section  IV  showing  the  generation  of  the  wire  grid. 

The  next  column,  labeled  "END1",  specifies  what  Is  the  lowest  numbered 
segment  attached  to  the  negative  end  of  the  subsection.  The  word  "negative" 
Is  used  In  the  sense  that  one  travels  along  the  subsection  from  this  end 
(END1)  or  point  to  the  other  end  (END2)  or  point.  There  Is  also  a signi- 
ficance attached  to  the  sign  associated  with  segment  numbers  listed  In  this 
column.  The  number  In  the  first  row  of  the  sample  output  has  no  sign,  while 
the  number  In  the  second  row  has  a negative  sign.  The  absence  of  a sign 
Indicates  a connection  with  the  positive  end  of  that  subsection,  and  the 
presence  of  the  negative  sign  shows  a connection  to  the  negative  end  of  the 
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subsection  Indicated.  Therefore,  the  negative  end  of  segment  one  Is 
connected  to  the  positive  end  of  subsection  30;  and  the  negative  end  of  sub- 
section two  Is  connected  to  the  negative  end  of  segment  358.  This  Is  shown 
In  Figure  70,  which  Is  an  expanded  view  of  part  of  the  wire  grid  model. 

The  arrows  Indicate  the  reference  direction  of  current  flow.  The  reason 
why  plus  one  Is  not  In  row  two  Instead  of  -358  will  be  explained  shortly. 

The  next  to  last  column,  ISEG,  simply  contains  the  subsection  number 
for  which  the  data  In  that  row  apply.  With  one  exception, the  numbers  In 
this  column  appear  only  once.  That  case  Is  Illustrated  In  the  row  just 
below  ISEG  equals  396.  The  segment  number  as  well  as  the  tag  for  what 
should  be  subsection  397  Is  zeroed  out.  This  Indicates  that  subsection  397 
Is  a duplicate  of  a preceding  segment.  That  this  duplicate  subsection 
Is  numbered  29  can  be  determined  by  comparing  Figures  38  and  45.  This  can 
also  be  verified  by  checking  the  list  of  duplicate  subsections  which  occurs 
just  before  the  segment  data  table.  In  the  case  of  the  sample  problem  there 
are  24  such  duplicate  pairs  of  subsections. 

The  final  column  In  the  segment  data  table  gives  the  lowest  numbered 
subsection  and  Its  end  which  Is  connected  to  the  positive  end  of  the  sub- 
section whose  row  Is  being  studied. 

GEMACS,  In  determining  the  data  for  the  columns  labeled  "ENDl"  and  "END2", 
always  looks  forward  from  the  subsection  number  under  consideration.  The 
first  larger  number  It  finds  Is  the  number  that  would  be  entered  In  the 
appropriate  column.  When  the  count  reaches  the  highest  segment  number, 

GEMACS  then  cycles  back  to  segment  one  and  starts  counting  upward  looking 
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for  a subsection  connected  to  the  point.  When  It  reaches  the  number  of  the 
subsection  listed  under  ISEG,  then  the  searching  ceases.  In  this  manner 
a circular  linked  list  Is  established.  This  will  aid  In  determining 
what  subsections  are  connected  to  a particular  point. 

As  an  example  of  Its  use,  one  of  the  things  to  be  accomplished  when 
setting  up  the  wire  grid  model  Is  to  have  the  base  of  the  antennas 
connected  to  an  Intersection  of  the  wire  grid.  This  can  be  determined  by 
use  of  the  circular  linked  list.  The  procedure  would  be  as  follows. 

First  of  all  locate  the  subsections  that  have  tags  of  three  and 
four.  Then  find  the  one  subsection  In  each  tag  group  that  has  one  end 
located  at  X equals  0.0.  The  search  narrows  down  to  segments  170,  171  and 

. 

156  to  161.  The  second  search  pinpoints  subsections  156  and  170.  The 

, 

(x,y,z)  coordinates  of  170  and  156  are  (0.0,  0.75,  0.45)  and  (0.0,  0.55, 
0.45)  , respectively.  Now  segment  170  Is  connected  to  the  positive  end  of 
172,  which  Is  connected  to  the  negative  end  of  subsection  173,  which  Is 
connected  to  the  minus  end  of  segment  135,  which  Is  finally  connected  to 
the  plus  end  of  169,  which  then  points  back  to  the  negative  end  of  170. 
Thus,  subsections  170,  172,  173,  135  and  169  all  meet  at  one  Intersection. 
Similarly,  It  can  be  determined  that  156,  162,  163,  127  and  155  all  meet 
at  one  point.  Thus,  the  two  antennas  are  each  located  at  an  Intersection 
of  the  wire  grid. 

Returning  back  to  Figure  68,  It  Is  seen  that  the  load  data  are  next 
printed  out,  If  any  such  data  exist.  The  sample  problem  output  shows 
the  name  of  the  data  set  (LOAD),  the  data  showing  the  segment  num- 
bers on  which  loads  are  placed  and  the  values  of  those  loads. 
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The  load  type  Is  "ZIMP",  Indicating  a lumped  load  resistance  and  reactance 
Is  also  shown. 


These  data  are  Immediately  followed  by  the  data  on  the  excitation.  The 
geometry  HUTTOP  Is  excited  by  the  values  contained  In  SOURCE.  The  excita- 
tion type,  location  and  strength  are  also  shown. 

The  output  that  appears  between  the  data  for  the  excitation  and  the  ta- 
ble of  values  of  current  traces  the  actual  task  execution.  The  tasks  are 
accomplished  In  the  order  In  which  the  task  commands  are  sequenced.  So, 
the  first  operation  that  occurs,  as  seen  In  the  sample  output.  Is  that  the 
Impedance  matrix  (ZIJ)  Is  being  filled.  Other  auxiliary  data  such  as  the 
names  of  the  geometry  and  load  data  sets  are  printed  out,  as  well  as  the 
electrical  parameters  of  the  problem. 

Then  It  Is  seen  that  GEMACS  extracts  the  banded  matrix  from  the  original 
Interaction  matrix  ZIJ  with  a bandwidth  of  500.  Also  printed  out  are  some 
statistics  relating  to  the  Interaction  matrix  for  the  center  column  of  that 
matrix.  The  band  dominance  factor  Is  the  ratio  of  the  band  norm  to  the  dif- 
ference between  the  column  norm  and  the  band  norm,  where  the  norm  Is  the 
square  root  of  the  sum  of  the  squares  of  each  element  In  that  column.  The 
column  norm  Is  obtained  by  using  all  of  the  elements  In  the  column,  while 
the  band  norm  Is  derived  using  only  those  elements  that  appear  In  the 
banded  matrix.  The  higher  this  ratio  the  better,  since  a higher  ratio  says 
that  more  of  the  larger  Interactions  are  Included  In  the  banded  matrix. 

The  next  operation  Is  the  decomposition  of  the  banded  matrix.  It  Is 
shown  that  diagonal  pivoting  Is  being  used  by  the  N after  the  "pivot  equals" 
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Identifier.  This  Is  the  only  pivoting  allowed  In  the  present  version  of 
GEMACS  because,  as  was  reported  In  one  of  the  BDM  quarterly  reports  [6], 
pivoting  Increases  the  bandwidth  required  and  this  would  decrease  the  effi- 
ciency. Also,  for  the  matrices  being  dealt  with,  pivoting  has  not  been  found 
necessary.  As  a matter  of  fact,  this  parameter  Is  not  even  an  Input  In  the 
decomposition  command  (LUD). 

Also  printed  out  are  the  maximum  and  minimum  diagonal  values  of  the  banded 
matrix,  as  well  as  the  ratio  of  these  two  values.  This  ratio  can  be  used 
as  a measure  of  the  111-condltloness  of  the  banded  matrix.  The  higher  the 
ratio,  the  more  unstable  Is  the  matrix  and  the  greater  the  effects  of  compu- 
ter round-off  error  on  the  accuracy  of  the  solution. 

The  next  data  printed  out  relate  to  the  BMI  solution  technique.  The 
equation  as  Input  In  the  command  is  printed  out  along  with  the  maximum 
nunber  of  Iterations  allowed,  the  convergence  criterion  chosen  and  the 
value  of  that  criterion  to  be  reached  before  Iteration  is  stopped.  The 
iteration  process  can  then  be  followed.  The  predicted  number  of  Iterations 
that  will  be  needed  and  the  values  of  each  of  the  convergence  criteria 
calculated  at  that  iteration  are  also  printed  out. 

Note  that  the  first  four  Iterations  show  that  zero  required  Iterations 
are  necessary.  GEMACS  does  not  make  such  a prediction  until  It  sets  up 
a history  of  convergence,  which  requires  four  Iterations  to  accomplish. 

After  this,  the  figure  for  the  number  of  predicted  Iterations  Is  based  on 
the  algorithm  derived  In  Section  F of  RADC-TR-75-272  [6]. 
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When  convergence  or  divergence  is  reached,  a message  Is  printed  Indicating 
this  fact.  The  final  values  of  the  various  convergence  criteria  are  given, 
and  It  can  be  seen  that  the  final  value  of  the  chosen  criterion,  the  pre- 
dicted relative  error  (PRE),  Is  less  than  the  value  specified  In  the  BMI 
, command. 

The  final  data  that  are  considered  part  of  the  standard  output,  as  listed 
In  Figure  68,  are  those  data  relating  to  antenna  Input  and  output  power  and 
transmitting  antenna  Input  Impedance.  These  are  listed  under  the  heading 
"Antenna/Load  Parameters".  For  each  subsection  that  is  excited  (170  In  this 
case)  the  magnitude  and  phase  of  the  input  Impedance  Is  printed  out  as  well 
as  the  power  In  watts  delivered  to  the  antenna.  If  there  Is  also  a load  at- 
tached to  this  subsection,  then  the  power  loss  In  this  load  In  watts  Is  also 
given. 

* If  the  subsection  Is  just  simply  loaded,  then  the  first  three  items  have 

a value  of  zero.  The  power  loss  In  the  load  Is  the  only  non-zero  entry 
for  this  subsection.  This  Is  done  for  subsection  156  In  the  output. 

One  final  point  to  note  Is  that  interspersed  throughout  the  execution 
trail  are  messages  Indicating  that  a checkpoint  Is  being  taken.  The  message 
and  Its  explanation  are  discussed  In  section  V.G. 

B.  Calculation  of  Coupling 

The  antenna  and  load  parameters  that  have  just  been  discussed  allow  the 
calculation  of  the  coupling  that  exists  between  any  pair  of  antennas  In  the 
system.  The  data  supplied  for  subsection  156  specify  the  power  delivered  by 
the  receiving  antenna  to  the  load.  As  can  be  seen,  this  value  Is  0.146 
milliwatts.  The  power  delivered  to  the  transmitting  antenna  (subsection  170) 


Is  given  as  6.22  milliwatts.  Calculation  of  ten  times  the  log  to  the  base 
ten  of  the  ratio  of  these  two  powers  gives  the  coupling  between  the  two 
antennas  as  -16.5  dB. 

Measurements  were  made  on  a similar  structure  [11].  Two  monopoles  were 
mounted  on  flat  plates,  one  of  which  measured  approximately  5.33  wavelengths 
on  a side.  The  other  plate  measured  approximately  11.4  wavelengths  on  a 
side.  The  transmitting  antenna  was  located  approximately  1.67  wavelengths 
from  the  nearest  edge  as  shown  In  Figure  9A  of  reference  [11].  The  recei- 
ving antenna  was  located  1.67  and  3.17  wavelengths  from  the  nearest  edges 
and  0.5  wavelength  from  the  transmitting  antenna.  The  measured  coupling 
was  about  -21  dB. 

As  discussed  In  the  report  several  reasons  can  be  cited  to  explain  why 
a discrepancy  can  exist  between  the  measured  and  predicted  data.  First, 
the  measurements  were  not  performed  In  a reflectionless  area;  and  stray 
radiation  could  cause  an  Interference  effect  to  exist  such  that  scattered 
fields  from  nearby  surfaces  could  cause  cancellation  of  the  direct  field 
from  the  transmitting  to  the  receiving  antenna. 

Secondly,  the  plate  In  the  example  Is  considerably  smaller  than  the  one 
on  which  the  measurements  were  made.  Also,  the  antennas  In  this  problem 
are  much  closer  to  the  edges  than  they  were  in  the  experimental  setup.  The 
edge  effects  can  cause  unaccounted-for  reflections. 

Third,  for  the  larger  antennas  It  was  difficult  to  keep  the  two  mono- 
poles perfectly  parallel  and  both  orthogonal  to  the  ground  plane.  The 
mathematical  model  has  both  of  these  characteristics  built  In.  Thus,  a 
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polarization  loss  factor  could  cause  a lower  value  of  coupling  to  be 
measured  In  the  experiment. 

It  should  finally  be  noted  that  the  measured  data  that  are  presented  In 
the  report  are  an  average  of  several  measurements  taken  on  different  days 
by  different  personnel.  The  measured  data  varied  by  one  to  two  dB  over 
the  course  of  the  experiment,  especially  when  the  longer  antennas  were  In- 
volved. 

Pending  more  careful  experiments  on  a setup  that  more  nearly  approximates 
the  mathematical  model.  It  can  be  said  that  GEMACS  provides  a good  engi- 
neering approximation  to  a very  tricky  problem  In  the  near  fields  of  collo- 
cated antennas. 

C.  The  PRINT  and  WRITE  Commands 

In  addition  to  the  data  that  are  automatically  output  by  GEMACS,  the 
data  associated  with  any  symbolic  name  can  also  be  obtained.  There  are  two 
commands  available  to  obtain  this  Information.  The  first  of  these  Is 
the  "PRINT"  command,  shown  In  Figure  71.  This  command  will  result  In  the 
entire  contents  of  the  file  being  printed  out.  This  command  Is  most  usually 
used  to  print  out  the  current  In  amperes  and  the  excitation  In  volts  per  meter 
on  each  of  the  subsections. 

As  can  be  seen  In  Figure  71  the  mnemonic  for  this  command  Is  the 
word  "PRINT".  This  Is  then  followed  by  a list  of  the  symbols,  whose  files 
are  to  be  printed  out.  In  the  sample  problem  the  voltage  (excitation)  and 
the  current  are  to  be  printed  out.  Therefore,  the  symbols  "SOURCE"  and 
"CURRENT"  are  specified  In  the  command. 
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There  are  two  formats  In  which  the  data  are  printed  out.  If  the  data 
are  real,  then  ten  values  per  line  will  be  printed  out.  If  the  data  are 
complex,  only  two  values  per  line  will  be  printed.  In  this  case  each  value 
will  have  Its  real  and  Imaginary  components  listed,  as  well  as  the  magni- 
tude and  phase.  Examples  of  the  complex  format  can  be  seen  In  the  output 
for  the  example  problem. 

Also  Included  with  the  data  Is  a heading  that  gives  the  type  of  data 
(real  or  complex)  and  the  symbolic  name  of  the  file  whose  contents  are 
being  listed.  The  lineage  that  Is  given  shows  what  other  symbolic  names 
are  related  to  the  one  being  considered.  Thus,  the  current  Is  derived  from 
the  banded  matrix  (BNDZIJ),  which  In  turn  Is  obtained  from  the  Interaction 
matrix  (ZIJ),  which  Is  finally  dependent  on  the  geometry  data  set  (HUTTOP). 
The  excitation  set,  SOURCE,  Is  only  derived  from  the  geometry  (HUTTOP). 

There  are  times,  however,  when  the  data  set  Is  too  large  to  print  all 
of  the  values;  and  yet  some  portion  of  the  data  needs  to  be  displayed  In 
order  to  check  on  the  consistency  of  the  data.  One  example  of  such  a data 
set  would  be  the  elements  of  the  Interaction  matrix.  For  the  sample  problem 
over  500,000  lines  of  data  would  be  printed  on  10,000  sheets  of  paper.  The 
WRITE  command,  shown  In  Figure  72,  can  be  used  In  this  situation. 

Using  this  command  part  of  the  data  associated  with  the  symbol  DS  Is 
output  on  the  file  specified  In  the  LU  Item  In  a binary  format,  which  can 
be  used  as  Input  by  other  programs.  The  default  value  of  the  LU  Item  Is  the 
computer  system  high-speed  printer.  If  the  printer  Is  used,  then  the  data 
are  output  In  the  standard  ASCII  character  set. 

Also,  If  the  output  Is  being  written  on  the  system  printer,  then  the 
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FILE  ID  Item  Is  Ignored,  even  If  some  value  Is  Input.  If  some  storage  file 
Is  used,  such  as  a card  or  paper  tape  punch,  then  a FILEID , If  specified, 
will  be  the  first  word  output  to  this  device.  If  no  FILEID  Is  specified, 
then  no  Identifying  word  will  be  written  on  the  output  file.  This  word 
can,  therefore,  be  used  to  Insure  that  the  proper  data  set  has  been  accessed. 

Since  the  data  can  be  thought  of  as  being  stored  In  a matrix  format.  It 
can  be  said  that  the  four  remaining  Items  specify  the  row  and  column 
limits  of  the  data  that  are  to  be  printed.  Any  combination  of  these  Items 
may  be  defaulted,  the  default  values  being  dependent  on  what  parameters 
have  been  defaulted.  The  values  of  R1  and  Cl  specify  the  first  row  and  first 
col  min,  respectively,  of  the  data  to  be  output.  Their  default  values  are 
each  one. 

The  values  of  R2  and  C2  are  the  last  row  and  column,  respectively, 
of  the  data  to  be  printed  out.  Their  default  values  are  the  values  Input 
for  R1  and  Cl,  respectively.  If  these  latter  values  have  been  specified. 

If  they  have  not  been  specified,  then  R2  and  C2  will  default  to  the  number 
of  rows  and  columns,  respectively.  In  the  matrix. 

It  then  follows  that  If  all  four  Items  have  been  defaulted,  the 
entire  matrix  will  be  printed  out,  just  as  If  the  PRINT  command  had  been 
used. 

The  example  shown  In  Figure  72  will  Illustrate  these  statements.  In 
this  example  It  Is  requested  that  part  of  the  Interaction  matrix  be  printed 
out  on  the  system  printer.  The  FILEID  then  Is  not  needed,  and  It  has  been 
defaulted.  The  first  row  to  be  printed  Is  row  ten.  Since  the  parameter  R2 
Is  defaulted,  then  only  the  data  In  the  specified  columns  of  row  ten  will 
be  printed  out  since  R2  will  default  to  Rl  In  this  case.  The  column  limits 
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are  one  and  650,  determined  by  the  fact  that  Cl  has  been  defaulted  (and 
hence  Is  one)  and  C2  has  been  given  as  650.  If  both  Cl  and  C2  had  been 
defaulted,  then  the  entire  row  would  have  been  printed  out. 

0.  Electric  Field  Output 

In  addition  to  determining  the  currents  on  the  structure,  GEMACS  also 
calculates  the  electric  field  due  to  those  currents.  This  can  be  done  at 
any  point  In  the  near-  or  far-fleld  region  of  the  radiating  structure. 

These  data  are  then  printed  out  In  a tabular  form  and  In  an  optional  graphic 
form.  The  command  that  accomplishes  these  functions  Is  shown  In  Figure 
73. 

The  symbol  represented  by  SDS  Is  the  name  of  the  data  file  In  which  the 
electric  field  data  are  to  be  stored.  This  Item  may  be  defaulted.  In 
which  case  the  data  are  not  stored  and  cannot  be  later  referenced.  Use 
of  this  option  In  no  way  affects  the  printing  out  of  the  tabular  or  graphic 
forms  of  the  electric  field.  The  advantage  of  the  default  option  Is  that 
It  reduces  the  number  of  data  storage  files  required  by  GEMACS. 

The  defined  symbol  In  parenthesis  following  the  command  mnemonic  EFIELD 
Is  the  name  of  the  data  set  containing  the  subsection  currents. 

The  U,  V and  W In  the  next  nine  Items  represent  the  geometric  variables 
In  the  three  conmonly  used  coordinate  systems:  X,Y  and  Z In  the  rectangular 
system;  R,e  and  Z In  the  cylindrical  system;  and  R,e  and  * In  the  spherical 
system.  The  symbols  Ul,  DU  and  U2  represent  the  Initial  value,  the  Incre- 
ment and  the  final  value  of  the  variable  U.  The  same  Is  true  of  the  three 
Items  relating  to  V and  W. 

Each  of  these  three  variables  may  represent  any  one  of  the  geometric 

I variables  In  any  one  of  the  three  coordinate  systems.  That  Is,  Ul  may 
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represent  either  XI,  Y1  or  Z1  In  the  rectangular  system;  R1 , el  or  Z1 
In  the  cylindrical  system;  or  Rl,  el  or  $1  In  the  spherical  system. 

Then,  VI  would  represent  the  Initial  value  of  one  of  the  two  remaining 
variables  In  that  coordinate  system;  and  W1  would  be  the  Initial  value 
of  the  third  coordinate  system  variable. 

This  freedom  of  representation  allows  the  user  to  Implicitly  specify 
how  rapidly  each  variable  Increments  with  respect  to  the  other  two 
coordinates.  U varies  the  least  rapidly,  V varies  more  rapidly;  and  W 
Is  the  Independent  variable  against  which  the  ratio  of  the  electric  field 
vector  to  the  maximum  field  Is  measured.  That  Is,  the  coordinate  W 
Is  stepped  through  the  range  Wl  to  W2  In  Increments  of  DW  before  VI 
Is  Incremented  by  the  value  of  DV.  Similarly,  the  coordinate  V is  stepped 
through  Its  range  before  the  value  of  U Is  changed. 

The  example  shown  in  Figure  74  will  make  this  clearer.  This  command 

will  print  and  plot  the  electric  field  behavior  In  the  near  field  of  the 
top  of  the  hut.  The  complete  output  data  are  shown  in  the  Appendix.  It 
can  be  seen  that  the  electric  field  vector  data  are  plotted  as  a function 
of  Y first  for  X * 0.0  and  Z * 1.4.  The  next  tables  and  plots  are  for 
X » 0.0  and  Z a 1.7  and  X * 0.0  and  Z * 2.0.  At  this  point  the  value  of  X Is 

incremented  to  ten,  and  the  value  for  Z reverts  back  to  the  original  value. 

There  then  follow  three  more  sets  of  data  corresponding  to  Z * 1.4,  1.7  and 
2.0,  respectively. 

The  optional  graphic  display  Is  controlled  by  the  following  six-choice 
Item.  If  this  Item  Is  defaulted,  then  only  a tabular  listing  of  the 
data  will  be  output  by  GEMACS.  If  one  of  the  six  choices  Is  present, 
then  the  plot  will  be  In  either  a rectangular  or  polar  form  with  axes 
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EFIELD  (CURENT) 


X1=0.0  DX=10 . 0 X2=10 , 0 
Zl=1.4  DZ=0 . 3 Z2=2 . 0 
Y1=0.0  DY=0 , 1 Y2=1.2 


LINLIN 


Figure  74.  THE  EFIELD  COMMAND  EXAMPLE 
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In  either  a linear  or  logarithmic  progression.  Each  of  the  six  choices 
Is  made  up  of  two  merged  mnemonics  as  follows: 

1.  LINLIN  - Both  the  dependent  and  Independent  variables  are  plotted 

linearly  on  a rectangular  graph 

2.  LINLOG  - The  Independent  variable  Is  plotted  linearly  and  the 

dependent  variable  Is  plotted  logarithmically  on  a 
rectangular  graph 

3.  LOGLIN  - The  Independent  variable  Is  plotted  logarithmically 

and  the  dependent  variable  is  plotted  linearly  on  a 
rectangular  graph 

4.  LOGLOG  - Both  the  Independent  and  dependent  variables  are  plotted 

logarithmically  on  a rectangular  graph 

5.  LINPLR  - The  Independent  variable  Is  plotted  as  a function  of 

angle  from  the  reference,  and  the  dependent  variable 
is  plotted  linearly  on  a polar  graph 

6.  LOG PL R - The  Independent  variable  Is  plotted  as  a function  of 

angle  from  the  reference,  and  the  dependent  variable 
Is  plotted  logarithmically  on  a polar  graph. 

In  each  of  these  plots  the  Independent  variable  will  be  one  of  the 
geometric  variables:  X,Y,Z  In  the  Cartesian  coordinate  system;  R,e,  Z 
in  the  cylindrical  coordinate  system;  and  R,e,$  in  the  spherical  coordi- 
nate system.  The  normalized  magnitude  of  the  electric  field  vector  will 
be  plotted  as  the  dependent  variable.  This  magnitude  Is  not  the  abso- 
lute value  of  the  electric  field  at  the  point  being  analyzed.  It  Is  the 
ratio  of  the  field  at  that  point  to  the  maximum  field  calculated  for  all 
■■■'*!■  Wwmbimrallcw  points  generated  by  the  cm— and.  Therefore,  1%  is .2* ... 
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normalized  value,  the  normalization  factor  being  printed  out  just  before 
the  tabulated  values  generated  by  the  EFIELD  command. 

Note  that  If  the  LOG  option  Is  chosen,  then  the  graphical  and  tabu- 
lated value  Is  20  times  the  common  logarithm  of  the  ratio  of  the  point 
field  strength  to  the  maximum  field  strength.  This  represents  the  power 
In  dB  that  the  field  strength  at  that  observation  point  Is  down  from  the 
maximum  value  of  field  strength  at  all  the  observation  points. 

The  simplest  case  Is  that  of  the  LINLIN  rectangular  plot  shown  In  Figure 
75  In  which  both  axes  vary  linearly.  In  this  case  the  Independent 
variable  Is  plotted  across  the  top  of  the  page,  and  the  dependent 
variable  Is  plotted  down  the  page.  The  origin  of  the  plot  Is  located  In 
the  upper  left-hand  comer  of  the  plot.  Figure  76  shows  the  same  data 
plotted  In  a LINLOG  format.  In  which  the  Independent  variable  Is  plotted 
linearly  across  the  top  of  the  plot;  and  the  dependent  variable  Is  plotted 
logarithmically  down  the  side  of  the  page. 

The  data  for  these  plots  are  shown  In  Figure  77.  The  top  line  shows 
the  normalization  constant  that  Is  used  to  obtain  the  absolute  value  of 
the  electric  field  vector  at  the  observation  points.  The  next  line 
gives  the  values  for  the  geometric  coordinates  which  remain  constant 
for  these  data,  In  this  case  X and  Z.  The  left-hand  column  shows  the 
variation  of  the  Independent  variable,  while  the  next  six  columns  give 
the  magnitude  and  phase  of  the  three  components  of  the  electric  field. 

The  right-hand  column  gives  the  power  ratio  In  dB  between  the  total  electric 
field  at  that  observation  point  and  the  normalization  factor. 

For  example,  at  Y ■ 2.0  meters,  the  total  electric  field  vector  Is 
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approximately  0.4  mV/m.  Dividing  this  by  the  normalization  factor 
yields  approximately  0.39,  which  Is  the  value  plotted  on  Figure  75.  Taking 
the  common  logarithm  of  this  value  and  multiplying  by  20  yields  -8.25  dB, 
the  value  that  Is  plotted  In  Figure  76. 

The  same  procedures  for  Interpreting  the  plots  and  tabulated  data 
can  be  used  when  studying  the  polar  plots  In  the  cylindrical  and  spherical 
coordinate  systems.  Figures  78,  79  and  80  show  the  LINPLR,  LOGPIR  and 
the  tabulated  values  respectively,  of  a system  being  analyzed  In  the  cylin- 
drical coordinate  system. 

Again  the  first  line  of  Figure  80  lists  the  normalization  factor,  and 
the  second  line  gives  the  values  of  the  constant  geometrical  coordinates. 
The  rest  of  the  table  Is  read  Identically  to  that  shown  In  Figure  77, 
except  that  the  Independent  variable  Is  now  theta  Instead  of  Y. 

In  the  plots  In  Figures  78  and  79,  the  origin  Is  located  In  the  center 
of  the  figure.  The  reference  axis  (the  X-axis)  Is  shown  In  each  figure 
along  with  the  direction  of  Increasing  theta,  where  theta  Is  the  angle 
around  the  Z-axIs  In  the  X-Y  plane  In  the  cylindrical  system. 

As  an  example,  Figure  80  Indicates  that  for  theta  equal  to  30.0  degrees 
the  magnitude  of  th$  electric  field  Is  approximately  0.35  mV/m,  which  Is 
approximately  0.32  of  the  maximum  electric  field  at  all  observation 
points.  Multiplying  this  by  the  cosine  of  30  degrees  results  In  an 
X-axis  value  of  0.275.  Multiplying  0.32  by  the  sine  of  30  degrees  gives 

a Y-axis  value  of  0.16.  These  results  are  shown  by  the  dotted  lines  In 

I * 

figure  78. 

Taking  20  times  the  comnon  logarithm  of  0.32  results  In  the  normalized 
dB  qaln  shown  In  the  right-hand  column  of  Figure  80.  In  order  to 


Figure  79.  THE  CYLINDRICAL  LOGPLR  PLOT 
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correlate  this  number  with  the  plotted  values  of  Figure  79  this 
normalized  value  must  be  algebraically  added  to  100  dB,  the  dynamic 
range  of  the  polar  plots.  The  result  for  theta  equal  to  30  degrees  Is 
approxalmately  90  dB.  Multiplying  this  first  by  the  cosine  of  30  degrees 
and  then  by  the  sine  of  30  degrees  yields  77.95  dB  and  45  dB  for  the  X- 
and  Y-axis  values,  respectively. 

There  are  two  points  to  be  noted  about  the  polar  plots.  First,  they 
are  not  square.  Therefore,  curves  of  constant  power  down  are  ellipti- 
cal rather  than  circular.  Second  the  sign  on  the  axis  values  should  be 
Ignored.  Some  values  are  preceded  by  a negative  sign  only  because  of  the 
limitations  Implicit  In  ANSI  standard  FORTRAN  coding. 

The  analysis  of  the  plots  and  tabulated  data  in  the  spherical  coordinate 
system  Is  Identical  to  that  just  presented  for  the  cylindrical  coordinate 
system.  The  only  point  to  note  is  that  In  the  spherical  coordinate 
system  the  Independent  variable  may  be  either  phi  or  theta.  The  angle 
phi  Is  measured  In  the  X-Y  plane  and  Is  positive  when  rotating  counter- 
clockwise around  the  Z-axIs.  The  reference  (or  zero-degree  direction) 
is  the  X-axis.  The  theta  angle  is  measured  from  the  Z-axIs  and  Is  posi- 
tive when  traveling  from  this  axis  to  the  X-Y  plane.  It  has  the  same  X- 
axls  reference  In  the  plots  as  was  used  for  all  the  angular  coordinates  dis- 
cussed. 

Finally,  when  angular  coordinates  are  used  as  the  independent  variable 
In  a rectangular  plot,  they  are  treated  the  same  as  any  of  the  linear 
Independent  variables.  They  are  plotted  on  the  X-axis  and  may  vary  lin- 
early or  logarithmically. 
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Error  Messages  and  Debug  Output 

There  Is  an  extensive  set  of  messages  available  to  the  user  that 
Is  printed  when  GEMACS  encounters  an  error  during  Input  processing  or 
during  the  performance  of  an  analysis.  These  are  automatically  printed 
out  without  the  need  for  a user  request.  In  addition,  a dictionary  will 
soon  be  available  as  a supplement  to  the  user's  manual  which  will  alpha- 
betically list  the  messages  as  they  appear  and  provide  Information  re- 
garding the  subroutine  that  generated  the  message,  the  cause  of  the  error 
and  the  steps  to  be  taken  to  eliminate  the  error. 

As  mentioned  before.  If  an  error  Is  found  during  the  reading  of  the 
Input  deck,  GEMACS  will  print  the  appropriate  error  message  and  continue 
processing  the  rest  of  the  Input  cards.  If  subsequent  errors  are  found, 
further  error  messages  are  printed.  However,  execution  of  the  analysis 
process  will  not  be  Initiated.  GEMACS  will  terminate  after  the  Input  pro- 
cessing has  been  completed,  and  It  will  print  out  the  contents  of  the  Input 
deck  with  the  error  message  Immediately  following  an  Improper  coitmand. 

If  the  error  should  occur  during  the  execution  of  the  analysis, 

GEMACS  will  terminate  the  analysis  at  that  point,  print  out  an  appropriate 
error  message  and  take  a checkpoint.  If  a checkpoint  conmand  (CHKPNT) 
has  appeared  In  the  command  stream  prior  to  the  command  that  Initiated 
the  operation  In  which  the  error  occurred.  Since  there  Is  a checkpoint 
file  available.  It  Is  possible  to  restart  from  the  command  at  which  the 
error  occurred  once  the  source  of  the  error  has  been  located  and  corrected. 
However,  keep  In  mind  that  If  a conmand  Is  changed,  the  erroneous  conmand 
and  all  subsequent  commands  must  be  wiped  out  when  restarting  and  replaced 
by  a new  command  stream  (refer  to  the  discussion  In  connection  with  Figure 


64), 
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In  addition  to  the  error  messages  a wealth  of  Information  Is  availa- 
ble regarding  the  actual  processing  that  goes  on  during  the  execution 
of  an  analysis.  These  data  can  be  obtained  through  the  use  of  the  DEBUG 
command,  shown  In  Figure  81.  The  mnemonic  for  this  command  Is  the  word 
DEBUG.  If  the  option  chosen  Is  "ON",  then  the  data  generated  during  the 
execution  of  the  following  command  Is  printed  out.  Included  In  this  print- 
out will  be  the  occurrence  of  any  major  operation,  such  as  the  storage  or 
retrieval  of  major  blocks  of  data.  Examples  of  such  output  can  be  seen 
In  example  two  of  the  GEMACS  user's  manual  [10]. 

If  the  option  "OFF"  Is  chosen,  the  program  and  the  output  return  to 
the  normal  mode  of  operation.  All  commands  between  DEBUG  ON  and  DEBUG 
OFF  will  be  affected  and  will  have  the  diagnostic  Information  printed 
out  during  the  execution  of  their  operations. 

If  the  "TRACE"  option  Is  chosen,  then  the  debug  operation  Is  turned 
on.  However,  In  addition  to  the  Information  normally  printed  out  GEMACS 
will  also  provide  the  user  with  Information  regarding  the  entry  Into  or 
the  exit  from  various  subroutines  during  the  execution  of  an  operation. 

This  will  aid  In  following  the  program  flow  and  can  be  of  use  In  tracking 
down  "bugs"  within  the  code. 

If  the  "STATS"  option  if  chosen,  then  statistics  are  collected  and  output, 
which  describe  what  subroutines  were  accessed,  how  often,  the  amount 
of  CPU  time  expended  In  each  subroutine  and  the  percentage  of  the  total 
CPU  time  spent  In  each  subroutine.  These  data  are  useful  only  when 
modifying  GEMACS  to  Increase  efficiency.  Also,  the  statistics  are  meaning- 
ful only  when  the  entire  analysis  has  been  accomplished  without  the  need 
for  a restart,  since  In  this  case  all  statistics  are  reinitialized  to  zero. 
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Figure  81.  THE  DEBUG  COMMAND 
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In  addition,  the  use  of  this  option  significantly  increases  the  cost  of 
an  analysis  since  a large  number  of  subroutine  calls  Is  Involved  In  the 
collection  of  these  timing  statistics. 

In  contrast  to  the  other  options  available,  once  the  STATS  option  Is 
turned  on  it  cannot  be  turned  off.  However,  statistics  are  automatically 
collected  whenever  the  TRACE  option  is  In  effect.  In  this  case,  when  the 
TRACE  option  is  turned  off,  the  collection  of  statistics  automatically 
ceases. 

If  option  "ILP"  Is  chosen,  then  the  debug  operation  will  be  on  during 
the  reading  of  the  Input  deck.  All  of  the  task  tables,  symbolic  names  and 
other  useful  data  will  be  output,  which  are  helpful  in  finding  "bugs" 

In  that  part  of  the  code  which  reads  the  commands  and  sets  up  the  execu- 
tion tables  based  on  the  data  contained  In  the  commands.  This  option  is 
Illustrated  on  pages  A-3  to  A-ll  of  the  Appendix. 

The  example  shown  In  Figure  81  is  taken  from  example  two  of  the 
GEMACS  user's  manual  [10],  The  only  command  that  Is  affected  Is  the  one 
that  Initiates  the  BMI  solution  process.  The  first  command  would  Initiate 
the  debug  operation,  while  the  last  command  would  terminate  the  operation. 
If  this  latter  card  were  not  present,  debug  output  would  be  printed  out  for 
all  commands  following  the  BMI  command.  Caution  must  be  exercised  when 
using  the  DEBUG  command,  since  It  may  cause  the  printing  of  large  volumes 
of  data.  This  is  especially  true  when  the  debug  Is  turned  on  during 
the  execution  o^  the  Interaction  matrix  generation  command  (ZGEN). 
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In  addition  to  the  BMI  matrix  equation  solution  technique  (Section  V) 
GEMACS  has  the  capability  to  solve  the  MOM  matrix  equation  (Figure  5)  using 
the  entire  interaction  matrix  by  a Gauss-Jordan  elimination  process.  The 
sequence  of  commands  shown  In  Figure  62  Is  very  nearly  the  same  as  would  be 
used  for  the  latter  solution  technique.  The  only  difference  Is  that  the 
BAND,  LUD  and  BMI  commands  would  be  replaced  by  the  commands  discussed  In 
this  section. 

There  are  essentially  two  equivalent  ways  to  Implement  the  full  matrix 
solution  technique.  The  simpler  of  the  two  is  to  use  the  SOLVE  command, 
shown  In  Figure  82.  The  three  Input  quantities  for  this  command  correspond 
one-for-one  to  the  quantities  In  the  MOM  matrix  equation,  shown  In  Figure  5. 
That  is,  DS1  represents  the  name  of  the  Interaction  matrix,  SDS  Is  the  sym- 
bolic name  of  the  set  of  currents;  and  DS2  Is  the  name  of  the  excitation 
matrix.  This  can  be  seen  directly  by  studying  the  example  shown  In  the 
figure.  Here  the  same  symbols  are  used  In  the  SOLVE  command  as  were  used 
In  the  BMI  command  (Figure  60). 

Use  of  this  conmand  results  In  the  lower/upper  decomposition  of  the  Inter- 
action matrix  followed  by  the  forward  el  1ml nation/backward  substitution 
process.  As  such  It  combines  the  LUD  and  BACSUB  commands,  a combination  of 
which  Is  the  second  method  of  Implementing  the  full  matrix  solution  process. 

The  decomposition  command  (LUD)  Is  the  same  as  Is  used  In  conjunction 
with  the  BMI  solution  process  and  Is  discussed  In  Section  V (Figure  59). 
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SOLVE  DS1  * SDS  - DS2 


SOLVE  ZIJ  • CURENT  = SOURCE 


Figure  82.  THE  SOLVE  COMMAND 
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The  forward  ellmlnatlon/backward  substitution  command  (BACSUB),  shown 
In  Figure  83  finds  the  solution  for  an  excitation,  given  a previously 
decomposed  Interaction  matrix.  Therefore,  In  Figure  83,  DS1  represents 
the  name  of  the  Interaction  matrix  after  It  has  been  decomposed,  SOS  Is 
the  name  of  the  current  vector;  and  DS2  Is  the  symbolic  name  of  the  exci- 
tation matrix. 

The  example  shown  In  Figure  83  lists  the  two  commands  as  they  would 
appear  In  a command  stream  after  the  Interaction  matrix  generation  command. 
The  first  command  would  decompose  the  Interaction  matrix,  ZIJ,  and  store 
It  In  the  file  labeled  LUDZIJ.  Then,  the  second  comnand  would  perform 
the  ellmlnatlon/substltutlon  process  using  LUDZIJ  and  the  excitation  matrix, 
SOURCE,  to  calculate  the  subsection  currents  and  store  them  In  the  matrix 
whose  name  Is  CURENT. 

Although  the  commands  shown  In  Figures  82  and  83  result  in  Identical  out- 
puts, they  each  have  their  own  advantages.  When  using  the  SOLVE  command, 
the  analyst  needs  to  generate  only  one  comnand  to  calculate  the  currents. 
Furthermore,  the  generation  of  the  decomposed  matrix  set  Is  automatically 
accomplished  by  GEMAC3  and  need  not  be  a concern  of  the  user. 

However,  If  the  problem  at  hand  requires  the  derivation  of  the  subsection 
currents  for  a number  of  nonllnearly  related  excitation  sets,  then  the 
commands  shown  In  the  example  In  Figure  83  should  be  used.  The  reason  is 
that  the  Interaction  matrix  and  Its  decomposed  triangular  matrices  are  not 
dependent  on  the  excitation  specified.  As  a result,  the  Interaction  matrix 
needs  to  be  decomposed  only  once  using  the  LUD  comnand,  the  resultant 
matrix  then  being  used  In  a series  of  BACSUB  commands,  one  such  comnand 
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